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PREFACE 


The favorable reception of the First Edition of this volume 
appears to have sustained the authors' belief in the need of a 
book on mathematics beyond the calculus, written from the 
point of view of the student of applied science. The chief 
purpose of the book is to help to bridge the gap which separates 
many engineers from mathematics by giving them a bird's-eye 
view of those mathematical topics which are indispensable in 
the study of the physical sciences. 

It has been a common complaint of engineers and physicists 
that the usual courses in advanced calculus and differential 
equations place insufficient emphasis on the art of formulating 
physical probletns in mathematical terms. There may also be a 
measure of. truth in the criticism that many students with pro- 
nounced utilitarian leanings are obliged to depend on books 
that are more distinguished for rigor than for robust uses of 
mathematics. 

This book is an outgrowth of a course of lectures offered by 
one of the authors to students having a working knowledge of the 
elementary calculus. The keynote of the course is the practical 
utility of mathematics, and considerable effort has been made to 
select those topics which are of most frequent and immediate use 
in applied sciences and which can be given in a course of one 
hundred lectures. The illustrative material has been chosen for 
its value in emphasizing the underlying principles rather than 
for its direct application to specific problems that may confront 
a practicing engineer. 

In preparing the revision the authors have been greatly aided 
by the reactions and suggestions of the users of this book in both 
academic and engineering circles. A considerable portion 
of the material contained in the First Edition has been rear- 
ranged and supplemented by further illustrative examples, proofs, 
and problems. The number of problems has been more than 
doubled. It was decided to omit the discussion of improper 
integrals and to absorb the chapter on Elliptic Integrals into 
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much enlarged chapters on Infinite Series and Differential 
Equations. A new chapter on Complex Variable incorporates 
some of the material that was formerly contained in the chapter 
on' Conformal Representation. The original plan of making 
each chapter as nearly as possible an independent unit, in order 
to provide some flexibility and to enhance the availability of the 
book for reference purposes, has been retained. 

I. S. S. 

E. S. S. 

Madison, 'Wisconsin, 

September, ld41. 
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CHAPTER I 

INFINITE SERIES 

It is difficult to conceive of a single mathematical topic that 
occupies a more prominent place in applied mathematics than 
the subject of infinite series. Students of applied sciences meet 
infinite series in most of the formulas they use, and it is quite 
essential that they acquire an intelligent understanding of the 
concepts underlying the subject. 

The first section of this chapter is intended to bring into 
sharper focus some of the basic (and hence more difficult) notions 
mth which the reader became acquainted in the first course in 
calculus. It is followed by ten sections that are devoted to a 
treatment of the algebra and calculus of series and that represent 
the minimum theoretical background necessary for an intelligent 
use of series. Some of the practical uses of infinite series are 
indicated briefly in the remainder of the chapter and more fully 
in Chaps. II, VII, and VIII. 

1. Fundamental Concepts. Familiarity with the concepts 
discussed in this section is essential to an understanding of the 
contents of this chapter. 

Function. The variable y is said to be a function of the variable 
X if to every value of x under consideration there corresponds at least 
one value of y. 

If X is the variable to which values are assigned at will, then 
it is called the independent variable. If the values of the variable 
p are determined by the assignment of values to the independent 

1 
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variable Xj then y is called the dependent variable. The functional 
dependence of y upon x is usually denoted by the equation* 

y = 

Unless a statement to the contrary is made, it will be supposed 
in this book that the variable x is permitted to assume real 
values only and that the corresponding values of y are also real. 
In this event the function f{x) is called a real function of the real 
variable x. It will be observed that 

(1-1) y = y/x' 

does not represent a real function of x for all real values of Xy for 
the values of y become imaginary if x is negative. In order that 
the symbol f{x) define a real function of Xy it may be necessary to 
restrict the range of values that x may assume. Thus, (1-1) 
defines a real function of x only if a; ^ 0. On the other hand, 
y = \/x^ — I defines a real function of x only if |a:| ^ 1. 

Sequences and Limits. Let some process of construction 
yield a succession of values 


•Cl, X2j Xzy * * * ) * * * > 

where it is assumed that every Xi is followed by other terms. 
Such a succession of terms is called an infinite sequence. Exam- 
ples of sequences are 


(o) 1, 2, 3, • • • , n, 

1 _ 1 1 _ 1 

W 2’ 4' 8’ 16' 


(c) 0, 2, 0, 2, 


,1 + (-!)”, • • • . 




Sequences will be considered here only in connection with the 
theorems on infinite series, t and for this purpose it is necessary 
to have a definition of the limit of a sequence. 

Definition. The sequence Xiy 3:2, • • • , • • • is said to 

converge to the constant L as a limit if for any preassigned positive 
number e, however smally one can find a positive integer p such that 

\xn — I/| < € for all n > p, 

* Other letters are often used. In particular, if more than one function 
enters into the discussion, the functions may be denoted by fi(x)y / 2 (a?), etc.; 
hy/(x), g{x)y etc.; by F(x)y G(x), etc. 

t For a somewhat more extensive treatment, see I. S. Sokolnikoff, 
Advanced Calculus, pp. 3-21. 
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For convenience, this definition is frequently written in the 
compact form 


lim Xn = L, 

n— ► «o 


and L is called the limit of the sequence. If a variable x takes 
on these successive values Xi^ x^, * Xn, ••• y then x is said 
to approach L as a limit. It follows from this definition that, 
of the sequences given above, (h) converges to the limit 0, whereas 
Ca) and (c) are not convergent. 

As an illustration, let the variable x assume the set of values 

Xi = 0.1, X2 = 0.11, Xs = 0.111, • • • . ^ 

It is easily seen that 

lim Xn = 

n— ♦ « U 


that is, corresponding to any « > 0, one can find a positive 
integer p such that 

\L - Xn\ = - X,| < 6 


for all values of n greater than p. Observe that 

1 _ 1 1 _ 1 1 _ _ 1 
g xi 90» g Xa 900 » • • • , 9 9 • lO"’ 


Hence, for any € that is chosen, it is necessary to demand that n 
be large enough so that 

1 _ 1 . 

9 “ 9 • 10" • 

The inequality is equivalent to ♦ 

9 • lO* > 

€ 

and, taking logarithms to the base 10,* 
log 9 + n > log i 


or 

n > — (log 9 + log «) = - log 9€. 

* From the definition of the logarithm, it follows that, if > B, then 
log A > log B. 
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Thus, if p is chosen as any ihteger greater than 1— log 9c|, the 
inequality 


1 

9 


Xn\ 


< € 


will be satisfied for all values of n greater than p. 

Infinite Series. Let ui, ws, * • * be an infinite sequence 

of real functions of a real variable x. Then the symbol 
00 

(1-2) V Un{x) = Ui{x) + U 2 {X) + • • • + Unix) + * * ’ 

n-1 

is called an infinite series, 

if, in (1-2), X is assigned some fixed value, say x = a:o, there 
results the series of constants 


(1-3) 5^ Unix^), 

n — l 

Denote by Sn(a;o) the nth partial sum, that is, the sum of the 
first n terms, of the series (1-3) so that 

Bnix^) = ni(aJo) + W2(a:o) + * • * + Unix^)^ 

As n increases indefinitely, the sequence of constants 

«i(a:o), 52(a;o), * • • , Sn(xo), • • • 

either will converge to a finite limit S or it will not converge to 
such a limit. If 

lim Sn(xo) = S, 

n— ♦ « 

the series (1-22 is said to converge to the value S for x = Xo.* 
If the series (1-2) converges for every value of x in some intervalf 
(a, 6), then the series is said to be convergent in^he interval (a, 6). 
As an example, consider the series 

(1-4) 1 -t- X “b X® -1- • • • + x'*"^ H- • • • . 

If X = (1-4) becomes 

* 2 ~ 4 ~ ~ 2"“"^ ^ ' 

* This limit S is usually called the sum of the eenes (1-3). 
t This means that x can assume any real value between a and h and that 
a and h can be thought of as the end points of an interval of the x-axis. 
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which is convergent to the value 2. In order to establish this 
fact, note that 

1 _i_ 1 . 1 .1 

8» = 1 + 2 + 4+ -- -+^i 

is a geometric progression of ratio 3^, so that 



Hence, the absolute value of the difference between 2 and Sn 
is which can be made arbitrarily small by choosing n 

sufficiently large. 

On the other hand, if x = — 1, the series (1-4) becomes 

which does not converge; for S 2 n = 0 and S 2 n-i = 1 for any choice 
of n and, therefore, lim Sn does not exist. Moreover, if a; = 2, 

the series (1-4) becomes 

1 + 2 + 4 + • • • + 2--1 + • • • , 

so that Sn increases indefinitely with n and lim Sn does not exist. 

n—* a® 

If an infinite series does not converge for a certain value of x, 
it is said to diverge or he divergent for that value of x. It will 
be shown later that the series (1-4) is convergent for — 1 < a; < 1 
and divergent for all other values of x. 

The definition of the limit, as given above, assumes that the 
value of the limit S is known. Frequently it is possible to infer 
the existence of S without actually knowing its value. The 
following example will serve to illustrate this point. 

Example. Consider the series 

+ • 

and compare the sum of its first n terms 

0-1 

s» - 1 + 21 + 31 -t" nl 

with the sum of the geometrical progression 
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« .1 
^n=l+2 + ^+ ‘‘’+^ 


The corresponding terms of Sn are never less than those of but, no 
matter how large n be taken, Sn is less than 2. Consequently, Sn < 2; 
and since the successive values of s« form an increasing sequence of 
numbers, the sum of the first series must be greater than 1 and less than 
or equal to 2. A geometrical interpretation of this statement may help 
to fix the idea. If the successive values of Sn, 

8i = 1 , 

«2 = 1 + ^, = 1.5, 

+ + + = 

«. = H-^ + |i + ^+ ^ = 1.717, 

are plotted as points on a straight line (Fig. 1), the points representing 
the sequence si, S 2 , * * * , Sn, * * * always move to the right but never 


progress as far as the point 2. It is intuitively clear that there must be 
some point s, either lying to the left of 2 or else coinciding with it, which 
the numbers Sn approach as a limit. In this case the numerical value 
of the limit has not been ascertained, but its existence was established 
with the aid of what is known as the fundamental principle. 

Stated in precise form the principle reads as follows: If an infinite 
set of numhers Si, Sa, • • • , «n, • • • forme an increasing sequence (that is, 
Sn > Sn, when N > n) and is such that every Sn is less than some fixed 
number M (that is, Sn< M for all values of n), then Sn approaches a limit 
s that is not greater than M (that is, lim Sn = s ^ M). The formulation 

n—* «o 

of the principle for a decreasing sequence of numbers Si, S 2 , * , 

Sn, * , which are always greater than a certain fixed number m, will 

be left to the reader. 

2. Series of Constants. The definition of the convergence of 
a series of functions evidently depends on a study of the behavior 
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of , series of constants. The reader has had some acquaintance 
with such series in his earlier study of mathematics, but it seems 
desirable to provide a summary of some essential theorems that 
will be needed later in this chapter. The following important 
theorem gives the necessary and sufficient condition for the 
convergence of an infinite series of constants: 

qo 

Theorem. The infinite series of constants S Un converges if 

n=» 1 

and only if there exists a positive integer n such that for all positive 
integral values of p 

|Sn+p — Sn| = |Wn+l + i^n+2 + * * ' + Un^J[ < €, 

where t is any preassigned positive constant. 

The necessity of the condition can be proved immediately by recalling 
the definition of convergence. Thus, assume that the series converges, 
and let its sum be /S, so that 

lim Sn — S 

n— > 00 

and also, for any fixed value of p, 

lim Sn+p = S. 

n—* 00 

Hence, 

lim (Sn+, — Sn) = lim (m»+i + u.+t + • • • + «.+,) = 0, 
n— » 00 n--> 00 

which is another way of saying that 

|Wn+l + Un+2 + * * ' + Un^p\ < € 

for a sufficiently large value of n. 

The proof of the sufficiency of the condition requires a fair degree 
of mathematical maturity and will not be given here.* 

This theorem is of great theoretical importance in a variety of 
investigations, but it is seldom used in any practical problem 
requiring the testing of a given series. A number of tests for 
convergence, applicable to special types of series, will be given in 
the following sections. 

It may be remarked that a sufficient condition that a series 
diverge is that the terms Un do not approach zero as a limit when 
n increases indefinitely. Thus the necessary condition for con- 
vergence of a series is that lim Un = 0, but this condition is not 

n— > w 

* See SoKOLNiKOFF, I. S., Advanced Calculus, pp. 11-13. 
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sufficient; that is, there are series for which lim Un = Obut wh^ch 

n— > 00 

are not convergent. A classical exainple illustrating this case 
is the harmonic series 

+ ' ■ ■ +1+ ■ ■ ■ ’ 


in which Sn increases without limi t as n increases. 

Despite the fact that a proof of the divergence of the harmonic 
series is given in every good course in elementary calculus, it will 
be recalled here because of its importance in subsequent con- 
siderations. Since 


1 


+ ■ 


71 - 1-1 71 - 1-2 


-t- 


-I- 


71-1-71 


> n 


2n 2’ 


it is possible, beginning with any term of the series, to add a 
definite number of terms and obtain a sum greater than 
If 71 = 2, 




71 » 4, 


1 + 1 - j - 1 _(- 1 



71 = 8, 
71 = 16, 


1+JL + 

9 ^ 10 ^ 





17 ^ 18 ^ 


-I- A 

^ 32 



Thus it is possible to group the terms of the harmonic series 

in such a way that the sum of the terms in each parenthesis 
exceeds and, since the series 

1 + I + I + 5 ' 

is obviously divergent, the harmonic series is divergent also. 
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3. Series of Positive Terms. Thils section is concerned with 
series of the type 

flO 

On = Ol + 02 + * * * + fln + • • ' , 

n-1 

where the On are positive constants. It is evident from the 
definition of convergence and from the fundamental principle 
(see Sec. 1) that the convergence of a series of positive constants 
will be established if it is possible to demonstrate that the partial 
sums Sn remain bounded. This means that there exists some 
positive number M such that Sn < M for all values of n. 
The proof of the following important test is based on such a 
demonstration. 

00 

Comparison Test. Let ^ an he a series of positive terms^ 

n-1 

00 

and let S bn he a series of positive terms that is known to converge, 

n-1 

00 

Then the series S <^n is convergent if there exists an integer p such 

n-1 

00 

thaiy for n ^ p, an ^ bn- On the other hand, if S Cn is a series of 

n-1 

positive terms that is known to be divergent and if an ^ for 
00 

n ^ p, then S ctn is divergent also, 

n-1 

Since the convergence or divergence of a series evidently is not 
affected by the addition or subtraction of a finite number of terms, the 

proof will be given on the assumption that p = 1. Let «» = ai + ai 

00 

+ • * • + ttn, and let B denote the sum of the series S bn and Bn its 

n — 1 

nth partial sum. Then, since a« ^ hn for all values of n, it follows that 
Sn ^ Bn for all values of n. Hence, the Sn remain bounded, and the 

oo 

series S an is convergent. On the other hand, if a« ^ Cn for all values of 

n-1 

00 00 

n and if the series 2) Cn diverges, then the series 2! an will diverge also, 

n— 1 n— 1 

There are two series that are frequently used as series for 
comparison, 
o. The geometric series 

« 

(3-1) a + ar + ar* + • • • + ar~ + * • • , 
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which the reader will recall* is convergent to ^ i'’! ^ 

and is divergent if |r| ^ 1. 

6. The p series 

(3-2) i + + 

which converges if p > 1 and diverges if p ^ 1. 

Consider first the case when p > 1, and write (3-2) in the form 

(3-3) 1 + (?; + p) + (5 + p + p + ^) + • • • 

+ ( 

where the nth term of (3-3) contains terms of the series 
(3-2). Each term, after the first, of (3-3) is less than the corre- 
sponding term of the series 

i4-2*-i-4*4-~4- • • • 4- 2^“^ - — - I- • • • 

or 

(3-4) 1 -f •+• ^2 p-i) 2 + ’ • * + ^2p~i)n-i + * * ' • 

Since the geometric series (3-4) has a ratio 1/2^^ (which is less 
than unity for p > 1), it is convergent and, by the comparison 
test, (3-2) will converge also. 

If p = 1, (3-2) becomes the harmonic series which has been 
shown to be divergent. 

If p < 1, 1/n^ > 1/n for n > 1, so that each term of (3-2), 
after the first, is greater than the corresponding term of the 
harmonic series; hence, the series (3-2) is divergent alsc. 

Example 1. Test the series 

i+s+p+-'-+p+---- 

T^e geometric series 

1 + ^, + ^, + •••+!:+••• 

* Since the sum of the geometric progression of ^ terms a + ar + ar* 

4. ... -I- jg equal to ^ (1 — r*). 

1 r 1 — r 


^ + 
( 2 ““*)** ' 


1 ) 
2 » - l)pj 
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is known to be convergent, and the terms of the geometric series are 
never less than the corresponding terms of the given series. Hence, the 
given series is convergent. 

Example 2. Test the series 

^ ^ log 2 ^ log 3 ^ log 4 ^ + log n + 

Compare the terms of this series with the terms of the p series for 

p = 1, 

^ + i + l + ■ ■ ■ +^+ ■ * ■ • 

The given series is divergent, for its terms (after the first) are greater 
than the corresponding terms of the p series, which diverges when 

p = 1. 

00 

Ratio Test. i/The series S- cin of positive terms is convergent if 


and divergent if 


lim = r < 1 

n— ♦ «o On 

lim ^ > 1. 

n~» 00 




lim the series may converge or diverge. 

n—* 00 On 


Consider first the case when r < 1, and let q denote some constant 
between r and 1. Then there will be some positive integer N such that 


Hence. 


On+l ^ 
On 


3 


for all N, 


oy+i < aifqf 

OjV+2 < on+iQ < ONq^i 

as+s < < ONq^i 


and 

Oat+i + Oy+a + Oiv+3 + ’ * * as{q + 3* + 3* + ' ’ • )• 

Since ? < 1, the series in the right-hand member is convergent; there- 
fore, the series in the left-hand member converges, also. It follows that 


the series 2) a» is convergent. 

n — l 

If the limit of the ratio is greater than 1, then a„+i > On for every 

00 

n ^ JV so that lim On 5 *^ 0, and hence the series S o» is divergent. 

• n - 1 
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It is important to observe that this theorem makes no reference 
to the magnitude of the ratio of an+i/un but deals solely with the 
limit of the ratio. Thus, in the case of the harmonic series the 
ratio is On+i/un = n/{n + 1), which remains less than 1 for all 
finite values of n, but the limit of the ratio is precisely equal 
to 1. Hence the test gives no information in this case. 

Example 1. For the series 


1 



+ ~ + — + 


+ 



lim 

n-4 «0 ttn n— ► «0 


n + 1 , 2"-^ 
2^ n 


lim 

n--> «o 



1 

2 


and, therefore, the series converges. 
Example 2. The series 




i+A 

10 ^ 10 * 


+ — + 
^ 10*^ 


+ 



is divergent, for 


Um “=±1 = Um = lim ^ = «. 

n— ► CO An n— ♦ CO 10"'*'*' 71 1 n— ♦ "o 10 


Example 3. Test the series 
1 

1-2 

Here 

Qn-f 1 1 


JL 4- + 


1 


(2n — l)2n 


lim ^ 

n-» 00 An 


lim 


lim 


^ (2n — l)2n 

(2a + l)(2n + 2) ’ 1 


4n^ — 2» 


= lim 




n--» « 4n* + 6n + 2 n--> • 1 I § I 1 

2n 2n2 


== 1 . 


Hence, the test fails; but if the given series be compared with the 

p series for p = 2, it is seen to be convergent. 

00 

Cauchy’s Integral Test. Let S be a series of positive 

n-cl 

terms swh that a»+i < a». If there exists a positive decreasing 
function f{x), for x 1, such that f{n) = On, then the given series 
converges if the integral 

fix) dx 

exists; the series diverges if the integral does not exist. 
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The proof of this test is deduced eajsily from the following 
graphical considerations. Each term an of the series may be 
thought of as representing the area of a rectangle of base unity 
and height /(n) (see Fig. 2). The sum of the areas of the first 
n inscribed rectangles is less than so that 

8n+l — ai ^ 02 + Us + * • * + On+1 ^ ^ /(^) 

But f(x) IS positive, and hence 

fix) tte < " fix) dx. 

If the integral on the right exists, it follows that the partial sums 
are bounded and, therefore, the series converges (see Sec. 1). 



The sum of the areas of the circumscribed rectangles, Oi + 02 
+ • • • + On, is greater than f{x) dx; hence, the series will 
diverge if the integral does not exist. 

Example 1. Test the harmonic series 


1+5 + 1 + 


In this case, f{x) = - and 

X 



+ • • • 



and the series is divergent. 


lim log n 

n-4 « 


‘‘O » 
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« 1 

Exam^ple 2. Apply Cauchy's test to the p series 2/ where p > 0. 

n*l ' 

Taking /(x) = — » observe that 

jpp 

J ’" dx 1 , " ^ 

, “ 1 — * ' > if p 1, 

1 j;p 1 — p l’ 

= log *1", if p = 1. 

00 

^ ~ exists if p > 1 and does not exist if p ^ 1. 

PROBLEMS 

1. Test for convergence 

. . 1 . 1 . 1 , 1 . _ . 

+ ■ ■ ■ +4+ ■ ' ■ ’ 

W 1 + 11 + 1 + • • • ; 

1 . 1 • 2 . 1 • 2 • 3 , 

2^ 3^ + ■ ■ ■ » 

«' + l + S + 7+-'-’ 

/ \ 1 I 2 . 3 I 
(?) 2 ^ ^ » 


ih) 


1 +_J_ + 

I o 1 o ' 


1 


+ 


2 log 2 ' 3 log 3 ' 4 log 4 
2. Use Cauchy’s integral test to investigate the convergence of 




1 


1 -hn* 


+ 


(*) f + + 2T1* + 


3. Show that the series Z an of positive terms is divergent if nan 

n — l 

has a limit L which is different from 0. Hint: Let lim nan = L so that 

n— ► 00 

nan > L — € forn large enough. Hence, a» > 
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§4 


4. Test for 

convergence 

(a) 

•0 

1 

y/n* 


W 

eo 

n ^ 

2n> 


(c) 

ee 

1 

y/n(n 

+ 1)’ 

id) 

00 

n\ 

n*’ 


(e) 

00 

1 


(f) 

00 

n2 

nV 


is) 

1. 

nl ^ 
10^' 


ih) 

.1 

1 

n y/n* 


ii) 

00 

1 

(2n + 

T)*’ 

ii) 

00 

n 

(2n + 

w 


4. Alternating Series. A series whose terms are alternately 
positive and negative is called an alternating series. There is a 
simple test, due to Leibnitz, that establishes the convergence of 
many of these series. 

Test for an Alternating Series. If the alternating series 
Oi — a2 + Gs — G4 + • • • > where the ai are positive^ is such 
that ttn+i < an and lim an = 0, then the series is convergent, 

n— > « 

Moreover^ if S is the sum of the series, the numerical value of the 
difference between S and the nth partial sum is less than Gn+i. 
Since 

S2n == (Gi — Ga) + (G3 — G 4 ) + * • * + (tt2n-l a2n) 

— dl (d2 ds) • • • (d2n—2 d2n— 1 ) d2n> 

it is evident that S2n is positive and also that S2n < di for ah 
values of n. Also, S 2 < S4 < Se < • • * ; so that these partial 
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sums tend to a limit S (by the fundamental principle). Since 
s*»+i = «2n + osn+i and lim 02«+i = 0 , it follows that the partial 

n-+ «o 

sums of odd order tend to this same limit. Therefore, the series 
converges. The proof of the second statement of the test will 
be left as an exercise for the reader. 



is convergent since lim ^ = 0 and — < i- Moreover, S 4 = 1 

n-+ 00 W Ti \ U 

““ H H K differs from the sum S by less than 
Example 2. The series 

i-i+i-i+i-i+i 

2 ^ 3‘ 4^3* 6^3’ 

is divergent. Why? 

6 . Series of Positive and Negative Terms. The alternating 
series and the series of positive constants are special types of the 
general series of constants in which the terms can be either posi- 
tive or negative. 

Definition. If ui + U2 + • • • + Wn + • * • is an infinite 
series of terms such that the series of the absolute values of its terms, 
\ui\ + IW 2 I + • • • + \un\ , is convergent, then the series 

Wi + W 2 + * • • + Wn + • • • is said to be absolutely convergent. 
If the series of absolute values is not convergent, but the given series 
is convergent, then the given series is said to be conditionally 
convergent. 

Thus, • 



is convergent, but the series of absolute values, 

is not, so that the original series is conditionally convergent. 

If a series is absolutely convergent, it can be shown that the 
series formed by changing the signs of any of the terms is also a 
convergent series. This is an immediate result of the following 
theorem: 



$8 


INFINITE SERIES 


17 


Theorem. If the series of absolute values 2 |Mn| is corwergerU, 

n- 1 

00 

then the series S Wn is necessarily convergent. 

n 1 

Let 


S» = + W2 + * • • + 

and 

tn = |t^l| + ItA2l+ * * * +|t^n|. 

If p« denotes the sum of the positive terms occurring in Sn and 
denotes the sum of the negative terms, then 




(5-1) Sn = Pn - qn 

and 

<• = Pn + ^n. 

00 

The series 2) \un\ is assumed to be convergent, so that 

n»»l 

(5-2) lim tn = lim (p, + g«) = L. 

n— > 00 n-^ « 


But pn and g« are positive and increasing with n and, since (5-2) shows 
that both remain less than L, it follows from the fundamental principle 
that both the pn and qn sequences converge. If 

lim Pn = P and lim qn = Q, 

n-+ « n—* 80 

then (5-1) gives 

lim 8n = lim (p» — gn) = P — 0, 

n— ♦ 89 n-* 80 

00 

which establishes the convergence of S Un. 

n«l 


Moreover, it can be shown that changing the order of the 
terms in an absolutely convergent series gives a series which is 
convergent to the same value as the original series.* However, 
conditionally convergent series do not possess this property. In 
fact, by suitably rearranging the order of the terms of a condi- 
tionally convergent series, the resulting series can be made to 
converge to any desired value. . For example, it is knownf that 
the sum of the series 


1 



• • + 


. . . 


* See SoKOLNiKOFF, I. S., Advanced Calculus, pp. 240-241. 
t See Example 1, Sec. 13. 
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is loge 2. The fact that the sum of this series is less than 1 and 
greater than 3^ can be made evident by writing the series as 

(^ ■ 0 G " i) G " 0 ’ 

which shows that the value of for n > 2; whereas, by 

writing it as 

^ “ G ■ 0 ■ (5 ~ 0 ■ ’ 

it is clear that Sn < 1 for n > 2. Some questions might be 
raised concerning the legitimacy of introducing parentheses in a 
convergent infinite series. The fact that the associative law 
holds unrestrictedly for convergent infinite series can be estab- 
lished easily directly from the definition of the sum of the infinite 
series. It will be shown* next that it is possible to rearrange 
the series 



so as to obtain a new series whose sum is equal to 1. The positive 
terms of this series in their original order are 


, 1 

1 1 1 

3’ 

5’ 7’ 9' ■ 

The negative terms are 

1 1 

1 1 

2’ 4’ 

6’ 8' 


In order to form a series that converges to 1, first pick out, in 
order, as many posit^e terms as are needed to make their sum 
equal to or just greater than 1, then pick out just enough negative 
terms so that the sum of all terms so far chosen will be just less 
than 1, then more positive terms until the sum is just greater 
than 1, etc. Thus, the partial sums will be 


Bi = 1 , 

1 1 1 

^*^2 2 ’ 

1 1 ^ 1 1 

»4=l-2 + 3 + g 


30’ 


* General proof can be constructed along the lines of this example. 
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8*=1- 

37 ^ 1 *“ 

38 * 1 ~ 


1 , 1 , 1 _ 1 
2 ■'■3 ■'■5 4 60 ’ 
1 + 1 + 1 _ 1 + 1 + 1 = 
2 ^ 3^5 4 ^ 7^9 

2 ^ 3^6 4 ^ 7^9 


1307 

1260’ 

1 ^ 1093 
6 1260’ 


It is clear that the series formed by this method will have a sum 
equal to 1. 

As another example, consider the conditionally convergent 
series 


(5-3) 


1 


J_ + _L_J_ 

V2 \/3 Vi 


Let the order of the terms in (5-3) be rearranged to give the 
series 

The nth term of (6-4) is 

1 ■ 1 1_ 

V4n — 3 V4n — 1 V^^ 

which is greater than 



But the series 2 &» is divergent, and it follows that the series 

n — l 

(6-4) must diverge. 

00 

Inasmuch as the series S |wn| is a series of positive terms, the 

n-l 

tests that were developed in Sec. 3 can be applied in establishing 

00 

the absolute convergence of the series S Wn. In particular, the 

n — 1 

ratio test can be restated in the following form: 
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Ratio Tbst. The aeries ^ Unia abaolvtely convergent if 

n-1 


and is divergent if 




lim ^ 

n-* 09 Un 


If the limit is unity, the test gives no information. 
Example 1. In the case of the series 

x^ x^ 

l+a: + j, + 3j + 4i+-'-, 


Um “i±l = Um 
n-*oo Un n->«o n! X^-^ 


lim ± = 0 

n— ► 00 n 


for all values of x. Hence, the series is convergent for all values of 
X and, in particular, the series 

^ n . 2* 28, 

1-2 + ^-3I+ *•* 
is absolutely convergent. 

Example 2. Consider the series 


_i_._j__.__i . 

1 - * ^ 2(1 - a:)* ^ 3(1 - as)* ^ 


« |(n + 1)(1 - *)“+» 


= Um 


l(n + 1)(1 - 


= lim 


-(l + i)(l-x) 


Therefore, the series will converge if 




< 1 or 1 < |1 — »1, 


which is true for « < 0 and for x > 2. 

For X s* 0 and for x * 2 the limit is unity, but if ^ » 0 the series 
becomes the divergent harmonic series 


1 +1 +1 + 




and if X » 2 there results the convergent alternating series 
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je 


'^+1-1+ ‘ + • • • • 

It follows that the original series converges for < 0 and for x ^ 2 
and diverges for 0 x < 2. 

6. Algebra of Series. The following important theorems are 
stated without proof:* 

Theorem 1. Any two convergent series 

C7 = + W2 + • • • + Un + • • • 

F = i;i + 1^2 + • • • + Vn + • • • 
can be added or subtracted term by term to give 
t/ + y = (wi + vi) + (i/2 + V 2 ) + • • * + {Un + t^n) + ' * * 
or 

U - V - {ui — Vl) + (U2 — V2) + • • * + (Un- Vn) + * * ‘ . 

If the original series are both absolutely convergent^ then the resulting 
series will be absolutely convergent also. 

Theorem 2. If 

= i/1 + i/2 + • • • + Wn + • • • 

and 

V ^ Vl + V 2 + • • • + Vn + * • • 

are two absolutely convergent series, then they can be multiplied like 
finite sums and the product series will converge to UV. Moreover, 
the product series will be absolutely convergent. Thus, 

UV = UiVi + U1V2 + U2V1 + UlVz + U2V2 + UsVi + • • • . 

Theorem 3. In an absolutely convergent series the positive 
terms by themselves form a convergerU series and also the negative 
terms by themselves form a convergent series. If in a convergent 
series the positive terms form a divergent series, then the series of 
negative terms is also divergent and the original series is conditionally 
convergent. 

Theorem 4. 7/ i/i + i /2 + • • • + i/n + • • • /« aw abso- 
lutely convergent series and if Mi, ikf* 2 , • * • , M^ ••• is any 
sequence of quantities whose numerical values are all less than some 
positive number N, then the series 

UiMi + U2M2 + • • • + UnMn + • • • 
is'^ absolutely convergent. 

♦ See SoKOLNiKOFF, I. S., Advanced Calculus, pp. 212-213, 241-245. 
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Example. Consider the series 
sin X sin 2x 


sin 3x 


V 2 \ ^33 

This series is absolutely convergent for all values of x, for the series 


1 


1 


p ~ ^8 + 33 ““ 
is absolutely convergent and |sin nx\ ^ 1. 

PROBLEMS 

1 . Show that the following series are divergent: 

, , 5 7,9 11 , , , , 2n + 3 , 

(«)9-Z + R- -« + -'- + + • • 


6 8 

1-1 4-1-1 4-^ -i -4- 
(6) 2 2* "*■ 3 32 4 42 


(c) 


1 


2^3 4^ 


2. Test for convergence, and if the series is convergent determine 
whether it is absolutely convergent. 




^"^3 3-6''"3-6-9 

,, 2 3.4 5 , 


3-5 1-3-5-7 


3 • 6 • 9 • 12 


+ 


1 2^34 

3. For what values of x are the following series convergent? 


/ % I / ,\ 1 1 

(o) ® - 2 + 3 ; 

w 1 21 41 61 + ■ ' ■ ’ 

(c) 1 - a: + a:* -»*+••• ; 

' 

4 . Determine the intervals of convergence of the following series- 
2x . 1/ 2x V . 1/ 2x 


(a) 


•*-5(jT4) ■‘"Ks 




* + 4 2 \» + 4/ ' 3 Va: + 4. 

(6) X + 21a:» + 31»* + 4b« + • • • ; 

/ s , . . »»(«» - 1) , . Mm — l)(m — 2) 

(c) 1 + MX 4 2 i ® +. H 

where m is not a positive integer. 


X* + 
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7. Continuity of Functions. Uniform Convergence. Before 
proceeding with a discussion of infinite series of functions, it is 
necessary to have a clear understanding of the concept of con- 
tinuity of functions. The reader will recall that a function 
f{x) is said to be continuous at a point a: = a:o if lim f(x) = f(xo) 

X—*XO 

regardless of how x approaches Xq, From the discussion of 
the limit in Sec. 1, it appears that this concept can be defined in 
the following way: 

Definition. The function f{x) is continuous at the point 
X = Xo if, corresponding to any preassigned positive number e, 
it is possible to find a positive number 8 such that 

(7-1) |/(x) — f{xo)\ < e whenever \x — Xo\ < 5. 

The foregoing analytical definition of continuity is merely a 
formulation in exact mathematical language of the intuitive 



concept of continuity. If the function f(x) is i^epresented by a 
graph and if it is continuous at the point x = Xo, then it is 
possible to find a strip bounded by the two* parallel lines x — xq 
+ 8 and a; = Xo — 5, such that the graph of the function will lie 
entirely within the strip bounded by the parallel lines y = f{xo) 
+ € and y = /(xo) — € (Fig. 3). But if the function is discon- 
tinuous at some point (such as x = xi), then no interval about 
such a point can be found such that the graph of the function will 
lie entirely within the strip of width 2t, where € is arbitrarily 
small. 

Definition. A function is said to be continuous in an interval 
(a, b) if it is continuous at each point of the interval. 
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If a finite number of functions that are all continuous in an 
interval (a, b) are added together, the sum also will be a continuous 
function in (a, 6). The question arises as to whether this 
property will be retained in the case of an infinite series of con- 
tinuous functions. Moreover, it is frequently desirable to 
obtain the derivative (or integi^al) of a function f{x) by means of 
term-by-term differentiation (or integration) of an infinite series 
that defines f(x). Unfortunately, such operations are not always 
valid, and many important investigations have led to erroneous 
results solely because of the improper handling of infinite series. 
A discussion of such questions requires an introduction of the 
property of uniform convergence of a series. 

It was stated in Sec. 1 that the series 

(7-2) ui(x) + U 2 (x) + • • • + Un(x) + • • • 

is convergent to the value S, when x = Xo, provided that 

(7-3) lim 8n(xo) = S. 

n— » • 

The statement embodied in (7-3) means that for any preassigned 
positive number «, however small, one can find a positive number 
N such that 

|Sn(a:o) — /S| < e for all n ^ N. 

If the series (7-2) is convergent for every value of z in the interval 
(a, b), then the series (7-2) defines a function S{x). Let ajo be 
some value of z in (a, h),'so that 

|s»(a:o) — /S(a:o)l < « whenever n ^ N. 

It is important to note that, in general, the magnitude of N 
depends not only on the choice of «, but also on the value of zo. 
This last remark may be clarified by considering the series 

(7-4) X + (x — !)» + (x — !)»* H- • • • 

-f- (aj — 1)*"“* -!-•••. 

Since 

*»(*) ® (® 1)* (a: — l)a;® -|- • • • -j- (a; — l)a:"“^ ^ 

it is evident that 

lim 8n(z) as lim x" = 0, if 0 x < 1. 
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Thus, S(x) = 0 for all values of x in the interval 0 ^ ® < 1, and 
therefore 

l8n(*) - <S(x)| = lx" - 01 = |X"1. 

Hence, the requirement that |Sn(x) — <S(x)| < c, for an arbitrary 
e, will be satisfied only if x” < c. This inequality leads to the 
condition 

n log X < log «. 


Since log x is negative for x between 0 and 1, it follows that it is 
necessary to have 


n > 


log*, 

log X 


which clearly shows the dependence of N on both e and x. In fact, 
if € = 0.01 and x = 0.1, n must be greater than log 0.01/log 0.1 
= — 2/— 1 == 2, so that N can be chosen as any number greater 
than 2. If € = 0.01 and x = 0.5, N must be chosen larger than 
log 0.01/log 0.5, which is greater than 6. Since the values of 
log X approach zero as x approaches unity, it appears that the 
ratio log c/log x will increase indefinitely and that it will be 
impossible to find a single value of iV which will serve for € = 0.01 
and for all values of x in 0 x < 1. 

It should be noted that the discussion applies to the interval 
(0, 1) and that it might be possible to find an iNT, depending on t 
only, if some other interval were chosen. If the series and the 
interval are such that it is possible to find an iST, for any pre- 
assigned €, which will serve Jor all values of x in the interval, then 
the series is said to converge uniformly in the interval. 

00 

Definition op Uniform Convergence. The series S Un{x) 

n — l 

is uniformly convergerU in the interval (a, h) if, for any e > 0, there 
exists a positive number N, independent of the value of x in (a, h), 
such that 

|)S(x) — Sn(a:)| < € for all n ^ N. 


The distinction between uniform convergence and the type of 
convergence exemplified by the discussion of the series (7-4) 
will become apparent in the discussion of the series 


(7 *6) 1 -f- X + X* + • • • + “b * • * > 

where — ^ ^ x 
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• 1 • 

Since 8n(x) = -z , it follows that 

X X 

S(x) = Ito ,.(x) - Um (j-L - 
Then, 

\S(x) - s„(x)l = |j^|, 
which will be less than an arbitrary € > 0 if 



Again, it appears that the choice of N will depend on both 
X and €, but in this case it is possible to choose an N that will 
serve for all values of a; in (—H, H). Observing that the ratio 
log €(1 — a;)/log \x\ assumes its maximum value, for a fixed e, 
when jc = — 3^, it is evident that if iV is chosen so that 


iog €(^) ^ 
log 3^ 


log Sc 
log 2' 


then the inequality (7-6) will be satisfied for all n N. 

Upon recalling the conditions for uniform convergence, it is 
seen that the series (7-5) converges uniformly for —3^ ^ x ^ 
However, it should be noted that (7-5) does not converge uni- 
formly in the interval ( — 1, 1). For, in this interval, the ratio 
appearing in (7-6) will increase indefinitely as x approaches the 
values ±1. The discussion given above shows that the series 
(7-5) is uniformly convergent in any interval (—a, a), where 
a < 1. 

It may be remarked that the series (7-5) does not even con- 
verge for a; = ± 1. For x = 1, it is obviously divergent, and 
when aj =5 — 1 the series becomes 


1 — — 

If —1 < a: < 1, (7-5) defines the function — , which takes 

X X 

the value 3^ when a; = — 1. 



INFINITE SERIES 


27 


§7 


As is often the case with definitions, the definition of uniform 
convergence is usually difficult to apply when the behavior of a 
particular series is to be investigated. There are available 
several tests for the uniform convergence of series, the simplest of 
which is associated with the name of the German mathematician 
Weierstrass. 

Theorem. (Weierstrass M Test). Let 

(7-7) Ui{x) H- U 2 {x) + • • • + Un(x) + • • • 

be a series of functions of x defined in the interval (a, 6). If there 
exists a convergent series of positive constants, 

Ml + M 2 + • • • + Mn+ * • • , 

such that lwi(x)| S Mi for all values of x in (a, b), then the series 
(7-7) is uniformly and absolutely convergent in (a, b). 

Since, by hypothesis, the series of is convergent, it follows 
that for any prescribed e > 0 there exists an N such that 

Mn+i + Mn +2 -!-•*•<€ for all n ^ iV. 
But \ui(x)\ ^ Mi for all values of x in (a, b), so that 

|Wn+l(iC) -t" Un+ 2 {x) + • • • I ^ Mn+l + -^71+2 -|- * • * <6 

for all n ^ A and for all values of x in (a, 6). Therefore, the 
series (7-7) is uniformly and absolutely convergent in (a, 6). 

The fact that the Weierstrass test establishes the • absolute 
convergence, as well as the uniform convergence, of a series means 
that it is applicable only to series which converge absolutely. 
There are other tests that are not so restricted, but these tests are 
more complex. It should be emphasized that a series may con- 
verge uniformly but not absolutely, and vice versa. 

Example 1. Consider the series 

sin x , sin 2x , , sin nx , 

—+—^+ • • • • • • • 

Since |sin nx\ 1 for all values of x, the convergent series 

i + i+. 

12 ^ 2* ^ -r ^2 ^ 

will serve as an M series. It follows that the given series is uniformly 
and absolutely convergent in any interval, no matter how large. 
Example 2. As noted earlier in this section, the series 

I + X + x^ + ••• + x^ + ••• 
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converges uniformly in any interval (—a, a), where o < 1. The 
series of positive constants 

l+a + a*+ ••• + + ••• 

could be used as an Af series in this case, since this series converges for 
a < 1 and \x*\ ^ a* for x in (—a, a). 


PROBLEMS 


1. Show that the series (7-4) is uniformly convergent in the interval 

( 0 , HI 

2. By using the definition of uniform convergence, show that 

_1 1 1 

a: 4- 1 (x + l)(a; + 2) (a; -f- n — !)(« + n) 

is uniformly convergent in the interval 0 ^ x 1. 

Hird: Rewrite the series to show that Sn{x) = — and therefore 

a; “I- n 

|S(x) - s.(a!)| = 

3. Test the following series for uniform convergence: 


(a) 


cos X , cos 3a; , cos 5x , 

-p- + -35-+-^+ • • • 


f 


sin 2a; , sin 4a; , sin 6a; , 

{b) -yTy + -gry + yy + • * * ; 

(c) 1 4- a; cos 6 x^ cos 2^ 4- a;* cos 30 4- • * • , |a;| ^ a;i < 1 ; 
. ^ cos 2a; cos 3a; cos 4a; 

W 2 ® ” 3 “ ' 4 ® * * ’ » 

(fi) 10a; 4- 10*a;* 4- 10»a;® 4- • • • . 


8. Properties of Uniformly Convergent Series. As remarked 
in the preceding section, the concept of uniform convergence was 
introduced in order to allow the discussion of certain properties 
of infinite series. This section contains the statements* of three 
important theorems concerning uniformly convergent series. 

Theorem 1. Let 


ui{x) 4 - Uiix) 4 - • • • 4 - Unix) 4 - • ‘ • 

be a series such that each Uiiz) is a corUinuous function of x in 
the interval (a, b). If the series is uniformly convergent in (a, 6), 
then the sum of the series is also a continuous function of x in (a, 6). 
* For proofs, see I. S. Sokolnikoff, Advanced Calculus, pp. 256-"262. 
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Corollary. A discontinuous function cannot be represented 
by a uniformly convergent series of continuous functions in the 
neighborhood of the point of discontinuity. 

Theorem 2. If a series of continuous functions j 

ui{x) + U 2 {x) + • • * + Un{x) + • • * , 

converges uniformly to S(x) in (a, b), then 

S(x) dx = ui(x) dx + f^U2(x) do; + • • • H- f^Un(x) dx + • • • . 

Ja Ja Ja 

where a < a < b and a < < b. 

Theorem 3. Let 


Ui(x) + U2{x) + • * • + Un{x) + • • • 

be a series of differentiable functions that converges to S{x) in 
(a, b). If the series 

^i(^) + + * * • + Wn(^) + • * • 

converges uniformly in (a, h), then it converges to S'(x). 

These theorems provide sufficient conditions only. It may be 
that the sum of the series is a continuous function when the series 
is not uniformly convergent. It is impossible to discuss neces- 
sary conditions in this brief introduction to uniform convergence. 
It may happen also that the series is differentiable or integrable 
term by term when it does not converge uniformly. In the 
chapter on Fourier series it will be shown that a discontinuous 
function can be represented by an infinite series of continuous 
functions. In that chapter, it is established that the series 

/ . sin 2x , sin Zx \ 

2(^sm*--^ + -3 ) 

represents the function a; for — tt < a: < t. But, if this series 
be differentiated term by term, the resulting series is 

2(cos X — cos 2a; + cos 3a; — * * • )> 


which does not converge in (— ir, tt); for the necessary condition 
for convergence, namely, that lim \un\ = 0, does not hold for any 

n— > « 

value of X. 

The series used in the first example of Sec. 7, 


,o , ,, sin a: , sin 2a: , sin 3x , 

(o-l; 12 02 o« I 


, Sin nx , 

I 2 ' 

n* 


2 * 


3* 
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is uniformly convergent in any interval (a, b) and as such defines 
a “continuous function S{x), Moreover, the series can be inte- 
grated term by term to produce the integral of S(a;). The 
term-by-term derivative of (8-1) is 

(8-2) c js X + }4 cos 2x + 3^^ cos 3a; -f * • • 

which is convergent in (0, tt), but the M series for (8-2) cannot 
be found since 1 + 3^^ + )^ + * * • is divergent. This merely 
suggests that (8-2) may not converge to the derivative of S(a;), 
but it does not say that it will not. 


PROBLEMS 

1 . Test for uniform convergence the series obtained by term-by-term 
differentiation of the five series given in Prob. 3 of Sec. 7. 

2. Test for uniform convergence the series obtained by term-by- ■ 
term integration of the five series given in Prob. 3, Sec. 7. 

9, Power Series. One of the most important types of infinite 
series of functions is the power series 


(9-1) 2^ ^ uo + aix -f a 2 X^ + . . . 4. QnX^ + • * • , 

n—O 

in which the are independent of x. Some of the reasons for 
the usefulness of power series will become apparent in the dis- 
cussion that follows. 

Whenever a series of functions is used, the first question which 
arises is that of determining the values of the variable for which 
the series is convergent. The ratio test was applied for this 
purpose in the examples discussed in Sec. 5. In general, for a 
power series, 

lim 

n— ♦ 00 

SO that the series converges if 


Un+1 

= lim 


Un 

n— ► • 

Un— 1 


and diverges if 


lim 


lim 


an-i 


Un-l 


< 1 


> 1 . 


Therefore, the series will converge for those values of x for which 


|x| < .lim 


an— 1 


an 
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If lim 


= r, it follows that the series will converge when x 


lies inside the interval ( — r, r), which is called the interval of con- 
vergence, the number r being called the radius of convergence- 

This discussion establishes the following theorem: 

00 

Theorem. If the series S is such that 

n—O 

lim l^i = r. 


then the series converges in the interval —r<x<r and diverges 
outside this interval. The series may or may not converge at the 
end points of the interval. 


Example 1. Consider the series 

/]f.2 /mS 

Since lim = lim I — ?1-~| 


+ ^ + 
n 


Since lim = lim -\ = 1, the series converges for 

— 1 < aj < 1 and diverges for |x| > 1. At the end point a? = — 1 the 
series becomes 



^2 3^4 

which is convergent. At the end point a; = 1 the divergent series 

i+i+i+i+i+ ■■■ 


is obtained. Hence, this power series is convergent for —1 ^ a; < 1. 
Example 2. The series 

1 + a; + 2!a;2 + • • • + n\x^ + * • • 

will serve to demonstrate the fact that there are power series which 
converge only for the value a; = 0. For 

lim 1^' = Um ~ -i’ = lim 1 = 0. 

Obviously, the series converges for a; = 0, as does every power series, 
but it diverges for every other value of x. 

Power series in a; — fc are frequently more useful than the 
special case in which the value of h is zero. A series of this type 
has the form 


ao + ai(a: — A) + a^ix — hy + • • * + dn(x — h)^ + • * • . 
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In this case the test ratio yields 


lim 

n—^ 00 


Wn+l 

= lim 

an 

Un 

n—* oo 

On-l 


\x — A|. 


If this limit is less than 1, the series is convergent; if greater than 
1, the series is divergent; and if the limit is equal to 1, the test 
fails and the values of which make the limit equal to 1, must be 
investigated. Thus, if the series is 




then 

lim 

n—* «o 


(x_^ 

^ ^ 




Un+1 

y 

(x - 1)* (n - 1)* 

Un 

— lim 

n— ► 

(x — 1)”“^ 


■ ir. (* ~ s - s) I* - ii - 1» - II- 


Therefore the series converges if — 1| < 1, or 0 < x < 2, and 
diverges for jx -- 1| > 1, or a; < 0, a; > 2. For a; — 1 = 1, or 
a; = 2, the series becomes 


1+1+^+^+ 


+ ^ + 


which is the p series for p = 2 and is therefore convergent. For 
a: — 1 = — 1, or x = 0, the series becomes 

1-1+ ^2 - p+ ••• + (-l)”p+ • • * . 
which is an alternating series of decreasing terms with Km Wn = 0 

n—> « 

and is therefore convergent. Thus the series is convergent for 
0 ^ a; ^ 2. 

, PROBLEM 

Find the interval of convergence for each of the following series, and 
determine its behavior at the end points of the interval: 

(а) 1 + a; + X? + «*+••• ; 

(б) l«a: + |-g+ • • • ; 

(c) 1 + (® + 1) + + • ' • ; 
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(d) 1 - 2a! + 3x’ - 4x» + • • • ; 

/»* /)(»3 

^ + 2 + ^ ^ + • • • ; 

(/) (x-2) -i(x-2)* + i(x-2)» - • • • ; 

/y*3 ^5 

10. Properties of Power Series. The importance of power 
series in applied mathematics is due to the properties given in 
the theorems of this section, as will be evident from the applica- 
tions discussed in succeeding sections. 

Theorem 1. If r > 0 is the radius of convergence of a power 

oe 

series S anX**, then the series converges absolutely and uniformly for 

n *0 

every value of x in any interval a ^ x < b that is interior to 
(-r, r). 

Since the interval (a, h) lies entirely within the interval (— r, r), it 
is possible to choose a positive number c that is less than r but greater 
than a and 6. The interval (a, 6) will then lie entirely within the 
interval (— c, c); and it follows that, for a ^ ^ 6, 

lanX"! < \anC^\, 

oe 

The series of positive constants S |a»c"| is convergent, for c < r, and, 

n»0 

accordingly, can be used as a Weierstrass M series establishing the 

00 

absolute and uniform convergence of S OnX" in (a, b), 

n — O 


Theorem 2. A power series 2 anX'^ defines a continuous func- 

n-'O 

tion for all values of x in any closed interval (a, b) that is interior to 
the interval of convergence of the series. 

This statement is a direct consequence of the preceding theorem 
and of Theorem 1, Sec, 8. 

Theorem 3. If the radius of convergence of the power series^ 
00 00 

5 CLnX^ is r, then the radii of convergence of the series S nanX'^'^^ 

n»0 n*0 


and 

n =0 


an 

n + 1 


obtained by term-by-term differentiaiion and 


integration of the given series^ are also r. 
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If the radius of convergence can be determined from the ratio 
test, then the proof follows immediately from the fact that if 


lim 

n--» 00 


fln-l 

an 


= r, 


then 


lim 

00 


(n — l)an-.i 
nan 

lim 


= r and 

(n + l)an-i 
nan 


r. 


Since the series obtained by term-by-term differentiation and 
integration are also power series, these processes can be repeated 
as many times as desired and the resulting series will be power 
series that converge in the interval ( — r, r). It follows from 
Theorem 1 that all these series are uniformly and absolutely 
convergent in any interval which is interior to (— r, r). How- 
ever, the behavior of these series at the end points x = —r and 
X — r must be investigated in each case. 

For example, the series 


/|•2 

l+^+|+|+ 


+ ¥ + 


has unity for its radius of convergence. The series converges 
for a; = — 1 but is divergent for x = 1. The series obtained by 
term-by-term differentiation is 


which has the same radius of convergence but diverges at both 
a; = 1 and a; = —1. On the other hand, the series obtained by 
term-by-term integration is 


. x^ . x^ . x^ . , , 

® 2-3 ^ n(n + 1) 

which converges for both a; = 1 and a; = — 1. 

This discussion leads to the conclusion stated in the following 
theorem: 

oe 

Theorem 4. A power aeries S anX^ may he differentiated and 

n *0 

integrated term by term as many times as desired in any closed 
interval (a, 6) that lies entirely within the interval of convergence of 
the given series, 

00 

Theorem 5. If a power series S anX^ vanishes for all values 

n — O 

of X lying in a certain interval about the point a; = 0, then the 
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coefficient of each power of x vanishes, that is, 

Cto “ Oj Cti ~ Oj flt2 ~ Oj * * * > * * * . 


The reader may attempt to construct the proof of this theorem 
with the aid of Theorem 2 of this section. 

11. Expansion of Functions in Power Series. It was stated in 
Theorem 2, Sec. 10, that a power series defines a continuous 
function of x in any interval which lies within the interval of 
convergence. This theorem suggests at once the possibUity of 
using such a power series for the purpose of computation. For 
example, the values of sin x might be obtained by means of a 
power series. Accordingly, it becomes necessary to develop 
some method of obtaining such a power series, and this section 
is devoted to a derivation of Taylor’s formula and a discussion of 
Taylor’s series. 

One of the simplest proofs of Taylor’s formula will be given 
here.* It assumes that the given function f(x) has a continuous 
nth derivative throughout the interval (a, 6). Taylor’s formula 
is obtained by integrating this nth derivative n times in suc- 
cession between the limits a and x, where x is any point in (a, 6). 
Thus, 


\f^^ix) dx = (a), 

Ja |a 

I I {dxy = j f^^~^^{x) dx — j dx 

Ja Ja Ja Ja 

= - {x - a)/(’‘-«(o), 

(dxy = (o) — (x — o)/<’*“*'(a) 


m x 

■ 




f ■ r 

Ja Ja 


/(”>(x) (dx)” = fix) - fia) - (x - o)/'(o) 
(x _ . . . _ 


2 ! 


-/"(a) - 


(n - 1)! 


/(>•-« (a). 


For other proofs, see I. S. Sokolnikoff, Advanced Calculus, pp. 291-295. 
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Solving for f(x) gives 

(ll-l) fix) = /(o) + (* - a)fia) + f'ia) 

+ • • • + + *- 

where 

(11-2) = X* ‘ 

The formula given by (11-1) is known as Taylor's formula and 
the particular form of Rn given in (11-2) is called the integral 
form of the remainder after n terms. The foregoing can be 
stated in the form of a theorem. 

Taylor's Theorem. Any function f{x) that possesses a con- 
tinuous derivative (x) in the interval (a, b) can he expanded in 
the form (11-1) for all values of x in (a, 6). 

The term Rn, which represents the difference between f{x) 
and the polynomial of degree n — 1 in a; — a, is frequently more 
useful when expressed in a different form. Since* 

J^®/(n)(x) dx ^ (x — a)f^^^{^), where a < ^ < x, 

repeated integration gives 

(11-3) • • £ /(»>(x) idx)- = /<»>({). 

The right-hand member of (11-3) is the Lagrangian form of the 
remainder after n terms. 

The special form of Taylor's formula that is obtained by 
setting a = 0 is known as the Maclaurin formula. In this case 

(11-4) fix) = m +mx +rm r, + • ■ ■ 

where 

♦ The student will recall from elementary calculus that 
f v(.x) dx = (h - aMi), 


where o < { < 6. 
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Taylor^s formula with the Lagrangian form of the remainder 
is often encountered in a somewhat different form, which results 
from setting a; — a = A. Since a < ^ < x, { can be written 
in the form a + dh, where 0 < d < 1, Hence, (11-1) becomes 

(11-5) /(a + h)= m +r(a)h+r{a) • • • 

+ 0fi) where 0 < d < 1. 

In this derivation of Taylor^s formula, it was assumed that 
f(x) possesses a continuous nth derivative, and as a result it 
appeared that then f(x) could be expressed as a polynomial of 
degree n in x — a. It should be noted, however, that only the 
first n coefficients of this polynomial are constants, for the 
coefficient of (x — a)" is a function of f and the value of f is 
dependent upon the choice of x. It may happen that f{x) 
possesses derivatives of all orders and that the remainder Rn 
approaches zero as a limit when n — » qo regardless of the choice 
of X in (a, b). If such is the case, the infinite series 

(11-6) /(«) ■l-/'(a)(x - o) +rM . 

+ /“(«) • • ■ 

is convergent and, in general,* it converges iof{x). 

The series given in (11-6) is called the Taylor's series expansion, 
or representation, of the function f{x) about the point x = a. 
The special form of (11-6) that is obtained when a = 0, namely, 

(11-7) m +rio)x +no) + ■ • ■ +P’'Ko) J + • • • 

is called Maclaurin's series. 

Example, Find the Taylor's series expansion of cos x in powers of 

T 

^ 2 
Since 

fix) = 003 X, ^ ( 1 ) 
f'(x) = - sin X, ( 1 ) = 

* For a further discussion of this point, see I. S. Sokolnikoff, Advanced 
Calculus, pp. 29d-298. 
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r(x) = -cos®, /"(1)=0; 
/'"(x) = sin X, J'" (I) = 1; 

= cos X, = 0; 


it follows that the result is 


cosx= -(x-|) + |j(x-|)‘-i(x-|)‘+ • • • . 


Since it is often possible to obtain a power series expansion 
of a function f{x) by some other method, the question arises as 
to the relation of such an expansion to the Taylor^s series expan- 
sion for f(x). For example, a power series expansion for the 


function 



is obtained easily by division, giving 


= l + a; + x2+ ••• + + •••, 


The reader can check the fact that the Maclaurin expansion for 
this function is identical with the power series obtained by 
division. That this is not an exceptional case is established in 
the following theorem: 

Uniqueness Theorem. There is only one possible expansion 
of a function in a power series in x — a; and, therefore, if such an 
expansion be found in any manner whatsoever, it must coincide with 
Taylor's expansion about the point a. 

Suppose that fix) could be represented by two power series in 
X — a, so that 

fix) = ao + aiix — a) + a2(x — a)^ + * • • 

+ anix - a)~ + • • • 

and 

fix) = 6o + biix — a) + biix ~ a)* -}- • • • 

+ bnix -«)»*+•••. 

Since both these expsuisions represent fix) in the vicinity of a, 
there must be some interval about the point a; = a in which 
ooth the^ expansions are valid. Then, in this interval. 
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2) an(.x - a)» = 2) ~ ®)”» 


n“0 


n — 0 


or 


^ (a„ — 6„)(a; — a)“ = 0. 


n*0 


It follows from Theorem 5, Sec. 10, that 

On — bn = 0, (n = 0, 1, 2, 


), 


or 


)• 


dn — hny — 0 , 1 , 2 , 

Hence, the two power series expansions are identical. 

Taylor's formula is frequently more useful in a slightly modified 
form. Let 

X — a ^ hy 

so that 

X = a + h. 

Then 

fix) = /(a) +f'ia)(x - a) ix-a)^+‘-- 


/"(a) 


becomes 

(11-8) /(a + h)= fia) + r(a)h + + • - 

/<»-»(a) I /<">(« + eh) 

in which 0 < ^ < 1, so that a<a + dh<a + hn 


PROBLEMS 

1. Find the expansion of each of the following functions in power 
series in x: 

(a) e*, (6) sin x, (c) cos x, (d) tan”^ x, 

(e) sin“^ ic, (/) sec a;, (gr) tan a?, W 

2. Expand 

(а) log X in powers of a; — 1; 

(б) tan X in powers of a; — j; 

(c) e* in powers of a; — 2; 

w 

(d) sin X in powers of a? — 

(c) 2 + — 3a;® + 7a;* in powers of a; — 1. 
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3, Show that sin x can be developed about any point a in a series 
(11-8) which converges for all values of h» 

4. Differentiate term by term the power series in x for sin x and thus 
obtain the power series in x for cos a?. What is the interval of con- 
vergence of the resulting series? 


6. Divide the series sin x 


x^ , «« 

® 31 + 6 ! " 


by the series 


cos X 


35* X^ 

l'2l+4i- 


, and thus obtain the series for tan x. 


6. Differentiate the series for sin"^ x to obtain the expansion in powers 

of 35 for (1 — Find the interval of convergence. Is convergence 

absolute? Investigate the behavior of the series at the end points of the 
interval of convergence. 

7. Establish with the aid of Maclaurin’s series that 


(a + 6)*** = k(l + 35 )** = A; 1^1 + tnx -1 — ^-- 2 ] — ^ x* + * • • 
where m is not a positive integer. 

This series is convergent for |x| < 1 and divergent when |x| > 1. A 
complete discussion of this series will be found in Sokolnikoff ’s Advanced 
Calculus. Some facts are: 

If X = 1, convergence is absolute if m > 0; 

If X = 1, convergence is conditional if 0 > m > —1; 

If X = —1, convergence is absolute if m> 0; 

If X = —1, series diverges when m < 0; 

If X = 1, series diverges when m ^ — 1. 


S.Letfiy)= S bny^&ndy= S a^x^^ be convergent power series. If 

n—O n—0 

<{y) is a polynomial, then the powers of y in terms of x can be determined 
by repeated multiplications and thus the expansion for f{y) in powers of x 
can be obtained. But if f(y) is an infinite series, this procedure may not 
be valid. Inasmuch as the power series in x is always convergent for x = 0 
and since the value of 2 / for x = 0 is ao, it is clear that the interval of con- 

00 

vergence of S Ky^ must include ao if the series for fiy)f in powers of x, 

n*i0 

is to converge. But if ao = 0, then f{y) surely can be expanded in 
power series in x by this method, for the point 0 is contained in the 

00 

interval of convergence of S 

n-O 

Apply this method to deriving the series in powers of x for » by 
setting 

. X* X* 

y«smx«x — + • •• 
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= l + y + |-,+ 



Explain why this method fails to produce the series in powers of x for 
log (1 + e*), where e* = y. 

12. Application of Taylor’s Formula. In this section two 
illustrations of the application of Taylor^s formula will be given, 
and in each case the remainder will be investigated to determine 
the error made in using the sum of the first n terms of the expan- 
sion instead of the function itself. 

1. Calculate the Value of sin 10®. Since 10° is closer to 0° 
than to any other value of x for which the values of sin x and its 
derivatives are known, the Maclaurin expansion for sin x will 
be determined and evaluated for x = 10° = tt/IS radian. Then 


/(*)=/(0)+/'(0)x+^** + 


• + 


/(»-i)(0) 

(n- 1)! 


|.n— 1 


^ n! 


where 0 < ^ < 

Since 

/(x) = sin X, 
/'(x) = cos X, 
/"(x) = — sin X, 
S"' {x) = — cos X, 


m = 0 ; 

/'(O) = 1; 
/"(O) = 0; 

r(0) = -1; 


therefore, 


/(»)(x) = sin ^x + /‘"HO) = sin 


x“ 


Here, 


sinx = x-J-, + ^j-f,+ ---+^/<»m 


R„{x) = ^ /<->(€) = g /‘")(0x), 0 < 0< 1, 


X" . , nir\ 


If only the terms through (or x®) are used in computing 
sin t/ 18, the error will be 
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SO that 


(ft) - (b)' ri “ (* b + 1 ') - (b)’ k b 

, <14 w 

' JL-JL- {iX 1 m {iXi - (ly 1 

“l8 18 \18/ 3!'’'\18/ 5! \18/ 7!’ 

\18j 9! 


with an error less than 


2. Compute the Value of It can be established readily, by 
expanding e* in Maclaurin’s series and evaluating for x = 1, that 
e = 2.71828 • • • . In order to compute the value of the 
expansion of e* about x = 1 will be used. The expansion is 


Six) = /(l) +/'(!)(* - 1) (*-!)*+••• 


. /'“>(«) 

n! 

Since 

fix) = e*, /(I) = e; 

fix) = /'(I) = e; 


ix - 1)». 


fMix) = e», f^"^il) = c; 

and 

/(»)({) = gt, 1 < { < 

therefore 

e* = e + e(x _ 1) + (x - 1)* + • . . + (* - l)“-‘ 

+ ^jix-l)’‘. 

Here 

^ (x - 1)» 


so that the error made in using only four terms is 
^4 ~ i! ~ 

i = e [l + 0.1 + M + |! (0.1)4 

= 1.105166e + 5:^ e*. 
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Thus, = 1.106166e with an error of (0.0001/24)e^ Since f 
lies between 1 and 1.1 and since is an increasing function, is 
certainly less than An approximate value of is 7, and the 
error is certainly less than 0.0007/24 = 0.00003. Therefore, 


= 1.1052e, 

correct to four decimal places. 

13. Evaluation of Definite Integrals by Means of Power 
Series. One of the most important applications of infinite series 
is their use in computing numerical values of definite integrals, 
such as c""®* dxj in which the indefinite integral cannot bo 
found in closed form. Moreover, the values of many tran- 
scendental functions are computed most easily by this method. 
Several examples of this use of infinite series are given in this 
section. 


Example 1. Consider 

log (1 + = /; = /; (1 + *)-* dz. 

Since 

(1 + 8 )-‘ = 1 — 2 + *»-«»+ • • • 
for \z\ < 1, it follows that 

+ *) = X* rri " X* + 



Example 2. Since 


sm" 


■s: 

-s:{ 

if \z\ < 1, therefore 


dz 


VT^ 


1 • 3 -4 4 - 1-3-5 


1 + - 2 * + Z* + 

^2 ^ 2-4 ^ 2 - 4-6 


«• + 


gin-1 a; ss aj 4- i 4- ^ ^ 4- 

^ 23 ^ 2 - 45 ^ 



It is evident that this method of obtaining the expansion of sin~^ x is 
much less complicated than the direct application of Taylor’s formula. 
Example 3. In order to evaluate the integral 


r 


2a dx 


y/2g{h — x){2ax — »*) 


I 
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expresi? it as 


fa fh dx 

■\gJo -v/te - x^\ 2o/ 

( X X 

1 — by its expansion in powers of giving 


+ Lili5fiV+ . . . 1 

^ 2 • 4 • 6 \2a/ ^ J 


dx 

\/kx x^ 


If this integral is expressed as 

/a r p dx n 

'ygL Jo Jo 2 \2a/ 


dx 


y/hx — x^ 

'o 2 *4X20/ 


\/Aa; — X* 


+ 


and each integral evaluated, there results 


]■ 


This expression gives the period of the simple pendulum. By making 
the change of variable x = h sin* <p, the integral reduces to 


I = 2 - ft* sin* 

where k* = h/2a. 

This is the form used in the discussion of the simple pendulum given 
in Sec. 71. In this illustration, h denotes the height of the pendulum 
bob and a the length of the pendulum. 

j '*! 0-9 

Q — dx cannot be evaluated by the 

usual method for evaluating a definite integral, for the indefinite integral 

0X ^ 0—X 

cannot be obtained. Moreover, the expansion for — - — > if obtained 

directly with the aid of Maclaurin's formula, would lead to an extremely 
complicated expression for each derivative. The expansion is most 
easily obtained by using the separate expansions for e» and e~*. Thus, 
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/ aj» X® N 

e. _e- = 2(x+^, + 5,+ ■ ). 


Hence, 


- 2 ( 1 + r?si. + 5^ +■•■)■ 2”«- 

Example 5. In order to evaluate the integral J'j e““ * dx, recall that 


11 I 1 I 


so that 


e«in * = 1 sin X + 


sin^ X . sin^ X 


2! 


+ 


31 


+ 


Then 


. . . sin^ic . sin®x . 

0 e^n.dx = Jg + sm * + + 


:£( 


1 + sin X 4- 


sin* X ^ sin® x 


+ 


2! ^ 3! 

which can be evaluated with the aid of the Wallis formula 


) dx 

• ) dx, 


TT IT 

r r^ ^ (n- 1)(« - 3) • • • 2 or 1 

sin-xrfx = J„ cos-xdx - 2) • • • 2 or 1 ' «> 


where a = 1, if n is odd, and a = 2 * if ^ is even. 

In order to justify the term-by-term integration, it is sufficient to 
show that the series in the integrand is uniformly convergent. That 
such is the case is obvious if one considers 

1 + 1 + + ^, + • • • 

as the Weierstrass M series. ^ 

PROBLEMS 

1. Calculate cos 10°, and estimate the maximum error committed 
by neglecting terms after x®. 

2. Find the interval of convergence of the expansion of e* in power 
series in x. Determine the number of terms necessary to compute 
accurate to four decimal places from this expansion, and compare the 
result with illustration 2, Sec. 12. 

3. Compute sin 33°, correct to four decimal places. 
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J '** dx . . . 1 

Q in power senes m x, and 

integrate term by term. Compare the result with that of Prob. 1(d), 
Sec. 11. 

6. Compute = 2(1 + H2)^9 correct to five decimal places. 

6 . Develop the power series in x for sin“i x and hence establish that 

. LUiV . . . . . 

6 2 2 3 \2/ 2 • 4 5 V2/ 

7. Show, by squaring and adding the power series for sin x and 
cos X, that 

sin^ X -h cos* x = 1. 

8. Evaluate by using series expansions of the integrands 
, s sin X dx 


/ s , /,x sin X dx 

(a) J„ sin (**) dx; (6) 

fo fo 

(e) f: cos (x*) dx; (/) X 


[2 - (1 - 2sin“|)] 

“ (^ + 2 1 ) 

(<') XI = r w r 

(i) e^^^dz* 

9, Show, by multiplication of series, that 

(1 + X + X* + • • • )' = 1 + 2x + 3x* + • • • 

= (1 - x)~*. 

10. Expand to terms in x^ 

(a) Vcosx; 

. sin X 

W ^mi; 

1 + e»’ 

11. Determine the magnitude of a, if the error in the approximation 

sin a = a is not to exceed 1 per cent. 

TT» 1 CL sin a . a® a® 

Evni: s= 0.01 and sina^a — — •••. 

« o! Ol 

♦ See form 787, B. O. Peirce, A Short Table of Integrals. 
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14. Rectification of Ellipse. Elliptic Integrals. In spite of its 
importance and apparent simplicity, the problem of finding the 
length of an elliptical arc is not usually considered in elementary 
calculus. This is because the integral that arises is incapable 
of evaluation in terms of elementary functions. However, the 
evaluation can be effected by means of series expansion of the 
integrand function, as will be shown in this section. 

Let the equation of the ellipse be 


?! 4 - = 1 

n2 


a >b. 


The length of arc from (0, b) to (xi, yi) is given by the integral 


(14-1) 


J ^X\ 
0 



Computing dy/dx and substituting its value in (14-1) gives 




Recalling the fac t that the numerical eccentricity of the ellipse 
is fc = integral given above can be written as 


(14-2) 






where fc* < 1. 

Let a: = o sin then dx = a cos 0 dO, and (14-2) becomes 


(14-3) 


= a y/l — sin^ B dB. 


The series expansion of the integrand function is most easily 
obtained by writing it as (1 — and expanding by 

use of the binomial theorem. Then (14-3) is replaced by 




L *2 sin2 B sin^ B 


and term-by-term integration* gives 

* Term-by-term integration is valid here, for the series 

1 , 1 jL. 1 . 1 • 3 • • • (2n — 3) , I 

H-2A:* + -- + ^ 

serves as a Weierstrass M series. 
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(14-4) « ■= ® “ 5 ** ^ ~ I X ~ 


1-3-6 


(2n - 3) 


2-4-6 


2n 


fc*" si 


sin^*" $ dS • 


]■ 


If (14-4) is used, it is possible to evaluate s for particular 
values of k and tp. However, the integral in (14-3) is so impor- 
tant that there are extensive tables* giving its value for many 
choices of k and <p. This integral for the value of a = 1 is 
called the elliptic integral of the second kind and is denoted by 
the symbol E{kj <p). It <p = t/2, the integral is called the com- 
plete elliptic integral of the second kind, which is denoted by the 
symbol E. 

The elliptic integral of the second kind having been defined, it 
seems desirable to mention the elliptic integral of the first kind, 
although the latter arises in considering the motion of a simple 
pendulum and will be discussed in more detail in Sec. 71. The 
elliptic integral of the first kind, F(fc, ^), has the form 


(14-5) F(k, v) - r , — 

Jo \/l — sin 2 ^ 

The complete elliptic integral of the first kind, which arises 
when ip = ir/2, is denoted by the symbol K, Values of F{ky <p) 
and of K are also tabulated, but the evaluation can be obtained 
from (14-5) by means of series expansion of the integrand. 
Thus, one has the expansion 


(14-6) F(fc, <p) = ip + 
1-3 




sin® 0 dff + g A:^ 


f 


+ 


(2n - 1) 


2-4 


2n 


si 


sin‘ 0 dd 
sin®“ 0d0 + • • 


16. Discussion of Elliptic Integrals. The elliptic integral of 
the first kind is a function defined by the integral 


(15-1) 


<p) m 


d0 I 

\/l — fc* sin* d 


< L 


* See the brief table in B. O. Peirce, A Short Table of Integrals, pp. 
121-128. 
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If sin 9 is replaced by z, (15-1) becomes 


(16-2) 


P(k,x) = 


dz 


\/(l - 2*)(1 - **«*)' 


fc* < 1. 


This is an alternative form of the elliptic integral of the first kind. 

Similarly, the same change of variable transforms the integral 
of the second kind 


(16-3) <f>) = JJ VI - A* sin* 9 d0, fc* < 1, 

into 

(16-4) M(k, x) = dz, k»<L 

It will be recalled that any integral of the type 
f R(x, V o** + bx + c) dx, 
where A is a rational function of the variables x and 


Vo** + 6* + c, 

is integrable in terms of the elementary functions, i.e., power, 
trigonometric, and logarithmic functions. It can be shown that 
the integration of integrals of the type 

(15-5) J R(x, V tt®* + bx^ + cx + d) dx 

and 

(15-6) J* R(x, V + 5x® -b ca:® -b dx -b e) dx 

requires, in general, the introduction of new functions obtained 
from the elliptic integrals. 

The evaluation of (15-5) and (16-6) can be reduced to the 
evaluation of integrals of the elementary types and the following 
new types: 

o. Elliptic integral of the first kind: 


Pik, x) s 


dz 


'o V(1 - 2^)(1 - 


or 


b. Elliptic integral of the second kind: 


de 


VI — fc* sin* 9 



or 

E(k, v) 



ifc* 8in*9d9. 
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c. Elliptic integral of the third kind: 


or 


’ “ Jo 

X 


dz 


n(n, ky <p) s 


(1 + m^) VCl - «")(! - A;V) 



(1 + n sin^ 6) \/l — siix^ 6 


The problem of reducing the integrals of expressions involving 
square roots of cubics and quartics to normal forms is not difficult, 
but it is tedious and will be omitted here.* Integrals involving 



square roots of polynomials of degree higher than the fourth lead, 
in general, to more complicated functions, the so-called hyper- 
elliptic functions. 

The graphs of the integrands of the integrals of the first and second 
kinds are of some interest (see Fig. 4). For A: = 0, 


M s Vl - k^ sin®^ 


and 


1 _ 1 

Vl' - k^ sin2 e 


both become equal to 1, and the corresponding integrals are both equal 
to For 0 < A? < 1, the curve y = 1/A6 lies entirely above the line 

y — 1 and the curve y ^ AS lies entirely below it. As if> increases, 


♦For a detailed account see Goursat-Hedrick, Mathematical Analysis, 
voi. 1, p. 226. A monograph. Elliptic Functions by H. Hancock, may also 
be consulted. 
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F{k, ip) and E{ky tp) increase continuously, F being ^tlways the larger. 
As k increases, ip being fixed, the value of F(A;, ip) increases and that of 
F(fc, ip) decreases. Also F(A;, t) = 2K and F(ik, ir) = 2F, for the curves 
are symmetrical about 6 — ir/2. If t/2 < ^ < it, it is obvious from 
the figure that 


(15-7) F(Aj, ip) =2K- F{k, TT - ip), 

E{kf = 2F — F(fc, TT — ^). 

Moreover, 

(15-8) F{kt mir ip) = 2mK + Fik, ^), 

slky mir + ^) == 2mE + E(ky ^), 


where m is an integer. 

Since the values of K and E, and of F(ky ip) and E{ky ip) for ip < ir/2, 
are tabulated, the relations (15-7) and (15-8) permit the evaluation of 
PiK ip) and E(ky ip) for all values of 

The discussion* above was restricted to values of A;* < 1. If /k® = 1, 
y = AS becomes y = jcos ^1 and y = 1/A^ becomes y = |sec ^|. 

Consider 


(15-9) 


u = 





\/(l - 2’*)(l - 



dd 

Vl - sin’ 0 


where x = sin ip. 

For a fixed value of A;, (15-9) defines u = P{x) or u = F{ip), The 
function resulting from the solution of (15-9) for ip in terms of u is 
called the amplitude of u and is denoted by am (w, mod A;), or more 
simply by ^ = am u. It will be assumed that the equation u = F{ip) 
can be solved for ip. Since p = am Uy 

a; = sin ^ s sin am w « sn w. 

Moreover, 

cos p s \/l — X* = V" 1 — sn^ u = cn u. 

Finally, 

Ap = y/l — k^x^ = dn u. 

The functions sn u, cn Uy and dn u are called the elliptic functions. 
From the definitions, it is obvious that 

am (0) = 0, sn (0) = 0, cn (0) == 1, dn (0) =* 1; 
am (— ti) = —am w, sn (— w) * — sn Uy cn (— t^) = cn Uy dn (— t^) 

= dn u. 

The elliptic functions are periodic functions and in some respects 
resemble the trigonometric functions. There exists a complete set of 

♦ See Prob. 1, at the end of this section. 
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formulas connecting the elliptic functioi^ analogous to the set for the 
trigonometric functions.* 


An interesting application of elliptic integrals to electrical 
problems is found in the calculation of the magnetic flux density 
in the plane of a circular loop of radius a carrying a steady 
current J. 

Upon applying the law of Biot and Savartf to a circular loop of 
radius a, the flux density B at any point P in the plane of the 


wire is given by 
(16-10) 


C A 

J c 


sin (r, s) da 


where C is the circumference of the loop, r is the radius vector 



Fig. 5. 


from P to an element of arc 
ds, and (r, s) is the angle 
between r and this element 
(Fig. 5). 

If the point P is at the 
center of the loop, then (r, s) 
= 90°, r = a, and the integral 
is easily evaluated to give 

#? 5 = jL ?![? — 1 
Air 2a 

a familiar result. 

If, however, the point P is 
not at the center, the evalua- 
tion of the integral is not so 


easy. Copsider the triangle RQSj where the side RQ = r dd 
makes an angle a with ds. It is clear that ds cos a = rdd; and, 
since a — 90° — (r, «), it follows that 


cos a = sin (r, «). 

Hence, 

.rdd 
ds = — 7 — ?• 
sin (r, 8) 

* See Appbl, P., and E. Lacoub, Fonctions elliptiques; Peircb, B. 0., A 
Short Table of Integrals; Gbbenhill, A. G., The Application of Elliptic 
Functions. 

tThis formula is known to engineers as Ampere’s formula. See, for 
example, E. Bennett, Introductory Electrodynamics for Engineers. The 
system of units used here is the '’rational” system of units used m M. Mason 
and W. Weaver, Electromagnetic Field. 
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53 


(15-11) 


for the magnetic flux density at P. 

Now, from triangle OQA, it is evident that 

‘s/a^ — (r sin 6)^ = r cos 6 -hA, 
which, after squaring both sides and simplifying, becomes 
+ 2rA cos 6 + {h^ — a^) = 0. 

Solving for r gives 

r = --h cos 6 ± \/A^ cos* $ + 

and, since r is always positive, the radical must be taken with the 
positive sign. Substituting this value of r in (15-11) gives 


£ p- 

4*- Jo -} 


4*- Jo -h cos e + cos* e + d* - A* 
or, upon rationalization of the denominator, 

/ cos 6 — \/h^ cos^ ^ 

Jo 

° 4^(a«^- fe*) ( X + X 

The first of these integrals is zero, and the second is an elliptic 
integral of the second kind, so that 




u r 

- A2) Jo 


^1 — sin^ 0 dBy 


where k = A/a. This integral can be evaluated for any k with 
the aid of the tables of elliptic integrals. 


PROBLEMS 


1. Prove that 


J r^ d(p ^ 1 da ^ ^ ^ 

0 \/l — Z* sin* ^ ^ Jo \/l — Z-** sin* a 

Hint: Change the variable by setting Z* sin* <p = sin* a. 
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2. Plot, with the aid of Peirce’s tables, F{kf ^), where k = sin a, 

using a as abscissa and P(fc, (p) as ordinate. Draw 10 curves on the 
same sheet of rectangular coordinate paper for ^ = 0, ^ = 10, <p = 20, 
^ = 30, ^ = 40, = 50, • • • , ^ = 90. 

3 . Plot four curves representing F{k, (p) on the same sheet of rectan- 
gular coordinate paper. Use <p as abscissa and the values of k as 0, 
Ht V3/2, and 1. 

4 . Plot the integrand of for the values of A: = 0, 

Jo Vl - k^ sin^ ip ’ 

H, and 1. Use (p as abscissa. The areas under the curves give the 
values of the elliptic integrals. 

6. Compute the value of P(0, t/2). 

6. The major and minor axes of an elliptical arch are 200 ft. and 
50 ft., respectively. Find the length of the arch. Compute the length 
of the arch between the points where x = 0 and x = 25. Use Peirce’s 
tables. 

7. Plot with the aid of Peirce’s tables E{ky <p), where k = sin a. 
Use a’s as abscissas and E{ky tp) as ordinates. Draw 10 curves on the 
same sheet of rectangular coordinate paper for ^ = 0, 10, 20, • • • , 
90. 


8. Plot on a sheet of rectangular coordinate paper the four curves 
representing E{ky <p). Use <p as abscissa. The four curves are for A; = 0, 
V3/2, and 1. 


9. Plot the integrand of \/l — fc* sin^ <p dtp for the values of 


A; = 0, and 1. Use ip as abscissa. Compare the result with that of 
Prob. 4. What can be said abou4»the relative magnitudes of F(A;, ip') 
and E{ky ip)l 


10. Show that 


kind. 




Vl + sin* ip 


is an elliptic integral of the first 


Hint: Change the variable by setting sin ^ ^ tan x. 

11. Show that 


Jo Vshi * Jo V cos X 


= ^/2 


cl -y 

Jo Vl H Sin* ip 


Hint: Set \/ cos x = cos (p. 

Note that the integral is improper but that it is easy to show its 
convergence. 

12. Show that 




e de (jr 
Vr~Fsfii> 0 fc’ ^ 



(16 


INFINITE SERIES 


65 


Hint: sin* ^ — p (1 — fc* sin* ff). 
13. Show that 


K s 


de 


fo^ Vl - li!’‘ sin* d 


k* 


/i^sy 

^ V2-4-6y 


fc« + 


if A* < 1. 


14. Find the length of one arch of the sine curve. 

16. Find the length of the portion of y = sin x lying between a; = 1 
and a? = 2. 

16. Given: 


^ (2’ 


dJd 


Find K and sn %K, 


17. Show that 


f 


de 


\/a^ cos e 


34 sin^ B 

> where a > 1, is an elliptic integral. 


18. Show that the length of arc of an ellipse of semiaxes a and h is 
given by 


T 


s = 4a \/l — sin^ 6 dO 

= 27ra where e is the eccentricity. 


16. Approximate Formulas in Applied Mathematics. It is 

frequently necessary to introduce approximations in order to 
make readily usable the results of mathematical investigations. 
For example, an engineer seldom finds it necessary to use the 
exact formula for the curvature of a curve whose equation is 
y = f{x), namely, 


(16-1) 



since in most applications the slope dy/dx is small enough to 
permit the use of the approximate formula 


(16-2) 



Many such approximations are obtained by using the first few 
terms of the Taylor’s series expansion in place of the function 
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itself. Thus, the formula (16-2) is obtained from (16-1) by 
neglecting all except the first term in the expansion of [1 + 
(dy/dxy]'“^^ in powers of dy/dx, 

1. Small Errors, The values of physical quantities determined 
by experiment are subject to errors due to inaccuracies arising in 
the measurements of the quantities involved. It is often neces- 
sary to know the size of such errors. 

Let a capillary tube contain a column of mercury. The radius 
R of the tube can be determined by measuring the length L and 
the weight W of the column of mercury. Let L be measured in 
centimeters and W in grams. Since the density of mercury is 
p = 13.6, 


R = 


Fw 

0.153 

\7rpL 



The principal error arises in the measurement of L, Let L be the 
true value, and let L' = L -h € be the observed value. Then, 
if R is the true value of the radius, let i?' = i? -f be the com- 
puted value. The error in measuring W is negligible because of 
the high accuracy of the balance. It follows that 



Since c is small compared with L, it follows that ri is approxi- 
mately given by ^ Clearly, € can be either positive or 
negative. 

2. Crank and Connecting Rod, If one end of a straight line 
PQ (see Fig. 6) is required to move on a circle, while the other 
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end moves on a straight line which passes through the center of 
the circle, the resulting motion is called connecting-rod motion. 
This kind of motion arises in a steam engine in which one end of 
the connecting rod is attached to the crank PB and therefore 
moves in a circle whose radius is the length of the crank, while the 
other end is attached to the crosshead and moves along a straight 
line. 

Let r be the length of the crank, I the length of the connecting 
rod, and s the displacement of the crosshead from the position A, 



Fig. 6. 


in which the connecting rod and crank lie in a straight line. 
Then, 

AB I + r, 

and 

AB = AQ + QD + DB ^ s + I cos (p + r cos $. 


Moreover, 

PD = Z sin ^ = r sin 6, 

so that 

. T . 

Sin ^ j sin d 

and 

Therefore, 

1 r* 

cos ^ = ^1 — sin* e. 

/ r* 

s + Z 1 

f 1 — ^2 sin* $) + r C09 6 = 1 + r 

or 

\ / 


5 = z[l — ^1—^ sin^ 1 ^ 


(l- 


If 
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be replaced by its expansion, it follows that 
* ~ ^ ^ sin® sin^ 0 + • • • j + r(l - cos «) 

= sin® ® + Ip sin* » + • * • ^ + r(l - cos B). 


if r is small compared with Z, the displacement of the crosshead 
is given approximately by r(l — cos d) . 

3. Surveying. In railroad surveying, it is fre- 
quently useful to know the amount of difference 
between the length of a circular arc and the 
length of its corresponding chord. 

Let r be the radius of curvature of the arc 
AB (Fig. 7), and let a be the angle intercepted 
by the arc. Then, if s is the length of the arc 



AB and c is the length of chord AB, s = ra and c = 2r sin 
Since 


sin a; = a; - cos f , 


where 0 < ^ < x, the error in using only the first two terms of the 


expansion is certainly less than j^jj* Then, 

(32 • 120) “ 


with an error less than 
2r 

Therefore, 

8 — c — ra 


ra^ 

1 ^’ 




with an error that is less than ra^/1920. 

4. Vertical Motion under EarWs Attraction. Let it be required 
to determine the velocity of a body of mass m that is falling from 
a height so above the center of the earth and is subject to the 
earth^s attraction alone. 

Let F be the attraction on the earth^s surface and F' be the 
attraction at a distance h from the surface (Fig. 8). Then 
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and, therefore, 




SO that 


„ kmmf , kmrn! 

p = _ and J” = 5rqrjgi- 

where m' is the mass of the earth, k is the gravitational constant, 
and r is the radius of the earth. Hence, 

F _{t + hY i 1 

F ' “ ' I I 

I ^ 

Also, let g be the acceleration at the surface of the j i 

earth and g' be the acceleration at a distance h above j 

the surface, so that F = mg and F' = mg\ It , 
follows that 

F _ g _(r + hy _s‘^ ( \ 

F' ~ g' ~ r* r®’ V J 


V.-_ 


This equation can be solved for v = ds/dt by the following 
device: Multiplying both members by 2 ds/dt and integrating 


= 

\dt ) 




where C is the constant of integration. If the initial velocity 
{ds/dt ) is zero, then C = — 2gfr2/so and hence 


^ 2grd^ - 1\ 

\dt) ^ \s sj 

ds/dt = Vy so that the eq 

- - )• 
\r + A 50/ 


But s = r + h and ds/dt = v, so that the equation becomes 


This formula can be used to calculate the terminal velocity 
{i.e.y the velocity at the earth^s surface) when the body is released 
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from any height. Thus, setting A = 0 gives 


(16-3) 




Upon denoting by ha the initial height above the earth^s surface, 
so that ^ + hoi (16-3) can be written as 


or 

(16-4) 

Now- 




r -t- Ao 
is permissible, so that 




and if y < 1, then series expansion 


ftt-'-f + CfT- •• 

T -T flo T \^/ 

Hence, if Ao/r < 1, (16-4) can be replaced by 

“*=^[ 7 -( 7 )’ + ( 7 )‘-- 



Moreover, if Ao is very small compared with r, then the powers 
of Ao/r higher than the first can be neglected* and 


= 2firAo, 

whichji^he fapiiliar formula for the terminal velocity of a body 
falling ^^b^rom a height Ao that is not too great. 

It follows from (16-3) that the square of the terminal velocity 
will be less than 2gr\\/r) = 2gr, Moreover, for large values of 
So the terminal velocity will be very close to y/ 2gr. Accordingly, 
if a body falls from a very great distance it would attain a ter- 
minal velocity (air resistance being neglected) of approximately 

'\/2gr = 6.95 miles per second. 


The results stated in the last paragraph may receive a different 
interpretation. Suppose a body were projected outward from 
the earth's surface with a velocity of more than \/2gr ^ 6.95 
miles per second. The previous discussion shows that, if air 

* Since the series is alternating, the error will be less than 2gr(ho/r)\ 
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resistance is neglected, the body would travel an infinite dis- 
tance. This velocity is called the critical velocity or the velocity 
of escape. 

It may be recalled that the earth’s rotation exerts a centrifugal 
force on a particle which is falling toward the earth and that this 
force diminishes the effect of the force due to the earth’s attrac- 
tion. For a particle of mass m on the surface of the earth at the 
equator, this centrifugal force is 


mv^ 

r 


r 


= = 


mg 

289 


dynes, 


where co = 0.00007292 radian per second is the angular velocity 
of the earth, r = 6,370,284 m., and g = 980 cm. per second per 
second. At a distance s from the center of the earth, this force is 


moih 


mgs 

28 ^' 


But the earth’s attraction at this distance is F = mg\ Since 
g' = 

p _ mgr^ 

.92 


If the particle is to be in equilibrium, 

mgs __ mgr^ 

SO that 

s* = 289r* or s = 6.6r = 26,000 miles approx. 

Thus, if all other forces are neglected, a particle would be in 
equilibrium at approximately 22,000 miles above the earth’s 
surface. This gives a very rough approximation to the extent 
of the earth’s atmosphere. The actual thickness of the atmos- 
pheric layer is supposed to be considerably smaller. 


PROBLEMS 

1. The mass of the moon is nearly one-eighty-first that of the earth, 
and its radius is approximately three-elevenths that of the earth. 
Determine the velocity of escape for a body projected from the moon. 
Acceleration of gravity on the surface of the moon is one-sixth that on 
the surface of the earth. 
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2 . Show that the time required for a body to reach the surface of the 
earth in Illustration 4, Sec. 16, is 


Hint: 


t = 




3. If the earth is considered as a homogeneous sphere at rest, then 
the force of attraction on a particle within the sphere can be shown to 
be proportional to the distance of the particle from the center. Let a 
hole be bored through the center of the earth, the air exhausted, and 
a stone released from rest at the surface of the earth. Show that the 
velocity of the stone at the center of the earth is about 5 miles per 
second. 

Hint: 


dH 


r 


where 5 is the distance of the stone from the center of the earth and r 
is the radius of the earth. 



CHAPTER II 
FOURIER SERIES 


17. Preliminary Remarks. It is frequently necessary to find 
the equation of a curve that passes through a certain number of 
given points lying in the xy-plojie. This can be accomplished in 
an infinite number of ways. Thus, if there are three given points, 
the coefficients of 

(17-1) 2 / “ Uo “h diX -f- (i2X^ 

can be chosen so that the resulting parabola will pass through the 
three given points. This is accomplished by solving the three 
linear equations, in ao, ^i, and ^ 2 , that arise when the coordinates 
of the given points are substituted in (17-1). If there are four 
given points, it is impossible, in general, to determine ao, ai, 
and a 2 so that the parabola (17-1) will pass through all four 
points, since the four linear equations in ao, ai, and a 2 will be, in 
general, incompatible. However, it will be possible to determine 
the coefficients of 

(17-2) y = CLo + aix + a 2 X^ + a^x^ 

so that the curve defined by (17-2) passes through all four points. 

The determination of the equation of a curve passing through 
a set of given points is not unique. Thus, if four points are 
given, it is possible to determine a curve whose equation is (17-2) 
which passes through them. But it is also possible to determine 
the coefficients of 

2 / = 6 o "f“ hiX ” 1 ” 62 sin x -h bzx^ 

so that the curve defined by this equation will pass through the 
given points. Obviously this curve will not coincide with that 
defined by (17-2). The type of curve can be varied at will, but the 
number of coefficients to be determined must equal the number of 
given points. 

If a curve is defined by the equation y == f{x), it is possible, 
as indicated above, to make the curve whose equation is 
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(17-3) y = ao + ai COS a: + ^2 COS 2a; + • • • + an cos nx 

+ fei sin a; + 62 sin 2a; + • • • + bn sin nx 

pass through any 2n + 1 points of y = f{x) in the interval from 
a; = 0 to X = 2ir. The question arises as to whether it is possible 
to make the curve y = /(x) coincide with that defined by (17-3) 
at all points of the interval (0, 2 w) by increasing indefinitely the 
number of terms in (17-3). It is already known that it is possible, 
under rather severe restrictions, to represent /(x) by an infinite 
series in powers of x. This was accomplished with the aid of 
Taylor’s series. The analogous problem of representing /(x) by 
an infinite trigonometric series was developed by Fourier and 
will be discussed in the succeeding sections. 

Whereas the representation of a function /(x) in Taylor’s series 
demands that f(x) possess derivatives of all orders, the develop- 
ment in a trigonometric series 


(17-4) 


2) («n cos nx + bn sin nx) 


y^f(x) 


is possible for a much larger group of functions. In fact, many 

periodic* functions having a 
finite number of ordinary dis- 
continuities can be represented 
by infinite series of trigonomet- 
ric functions. The term or- 
dinary (or finite) discontinuity 
is used to describe the situation 
that arises when the function 
/(x) suffers a finite jump at some 
point X = xo (see Fig. 9). 
Analytically, this means that the two limiting values of /(x), as 
X approaches Xo from the right-hand and the left-hand sides, exist 
but are unequal; i.e., 

lim/(xo + e) 9 ^ lim/(xo — c). 

€-►0 •—►o 



Xo-€ Xf, Xo+€ 

Fig. 9. 




In order to economize on space, these right-hand and left-hand 
limits are written as/(xo+) and/(xo— ), respectively, so that the 
foregoing inequality can be written as 

f{Xo+) 9^ f{XQ-). 

* A function is said to be periodic of period a if f{x) f(x + a). 
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Inasmuch as each term of (17-4) is a periodic function of 
period 2?r, it is necessary to restrict the discussion of the repre- 
sentation of functions by series of the type (17-4) to those func- 
tions which have period 2t, Or, what amounts to the same 
thing, the problem of representing a non-periodic function can be 
restricted to some inteiwal of width 27r, and the function defined 
outside this interval so that it is periodic. For the present, it 
will be assumed that the interval in which f(x) is considered is 
the interval (— tt, tt) and that outside this interval the function 
is defined by the equation f(x + 2m) = /(x). Of course, any 
interval (a, a + 2m) would do equally well. 

The theory of Fourier series is one of the most beautiful 
developments of analysis, and it serves as an indispensable instru- 
ment in the treatment of nearly every physical problem. Solu- 
tions of such important problems as sound vibration, propagation 
of electric currents and wireless waves, heat conduction, and 
mechanical vibrations give but a mere indication of its value. 

The following section contains the celebrated theorem giving 
conditions on f{x) that are sufficient to permit its representation 
by a Fourier series, and also a derivation* of the formulas for 
the Fourier coefficients (that is, the coefficients in the trig- 
onometric series). 

18. Dirichlet Conditions. Derivation of Fourier Coefficients. 

Theorem. Let f{x) be a function defined arbitrarily in the 
interval —m^x^Tj and outside this interval defined by the 
equation /(x + 27r) = /(x). If /(x) has a finite number of 'points 
of ordinary discontinuity and a finite number of maxima and 
minima in the interval —t^x^t, then it can be represented by 
the series 

CO 

^ (a* cos kx + bk sin kx), 

with 

1 f' 

o* = - I fit) cos kt dt, 

i)fc = i f fit) sin kt dt, ik — 0, 1,2, ‘ ' ), 

TT J-ir 

* For a more extended treatment, see I. S. Sokolnikoff, Advanced Calcu-. 
lus, Chap. XI. 
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which converges at every point x = xo of the interval to the value* 

/(aJo+) +/(a:o“~) 

2 

The restrictions imposed upon the function f(x) in this theorem 
are known as the Dirichlet conditions. 

The following demonstration that the Fourier coefficients au 
and bh have the form given in the theorem assumes that the 
Fourier series development of the function /(x), namely, 


(18-1) f{x) = ^ + 


can be integrated term by term in the interval (—tt, tt). The 
proof that if a function fix) satisfies the Dirichlet conditions 
then its Fourier series expansion actually converges to fix) is 
too involved to be given here.f 

In order to determine uo, multiply (18-1) by dx and integrate 
term by term from — tt to tt. Since 

cos nx dx = sin na; dx = 0 for n = 1, 2, • • • 


and, hence, 

J*’'^/(x) dx = aoTT, 

there results 

(18-2) 

1 r' 

ao = - 1 f{x) dx, 

^ «/--X 


The coefficient an of the general cosine term can be obtained 
by multiplying both members of (18-1) by cos nx dx and per- 
forming term-by-term integration from — tt to tt. Since, for all 
integral values of m and n, 


and 



sin mx cos nx dx 
cos mx cos nx dx 


= 0 


= 0 , 


for m 9 ^ n, 


• If f(x) is continuous at the point x — Xqj then/(xo+) = /(®o— ) = /(a^o), 
so that at all points of continuity the series converges to f{x). At the 
points of ordinary discontinuity, it converges to the arithmetic mean of the 
values of the right- and left-hand limits. 

t Knopp, K., Theory and Application of Infinite Series, p. 366; Cabslaw, 
H. S., Fourier’s Series and. Integrals, p. 207. 
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f(x) cos nx dx = an cos^ nx dx = anw. 

Therefore, 

1 

(18-3) an = - I fix) cos nx dx. 

^ J-ir 


It should be observed that (18-3) becomes (18-2) for w = 0. 

Similarly, by multiplying (18-1) by sin nx dx and performing 
term-by-term integration from — tt to tt, one obtains 


(18-4) 


bn 



sin nx dx. 


It can be shown that if the values of the function are not equal 
at the end points of the interval (— tt, tt), that is, if /(— tt) 5 *^ /(tt), 
then at these end points the Fourier series expansion for fix) 
converges to M[/(— ^) + fi^r)]. 

The student will convince himself that, if the function fix) 
is defined in the interval from 0 to 27r, then the coefficients an 
and bn in (18-1) are given by the formulas 


1 r 

an = - 
^ Jo 

19. Expansion of Functions in Fourier Series. This section 
contains some illustrative examples of expansion of functions, 
satisfying the Dirichlet conditions in the interval (— tt, tt), in the 
series 


1 

fix) cos nx dx and 6n = - I fix) sin nx dx. 


(19-1) 


90 

^ (an cos nx + bn sin nx). 


n-l 


where the coefficients an and 6« are given by the formulas 

1 T’" 

(19-2) an = - I fix) cos nx dx 

J-v 

and 

(19-3) 6n = - I fix) sin nx dx. 

Illustrative Example 1. Expand fix) = x in Fourier series in 
the interval — tt x ^ ir. Calculating the coefficients an and 
bn gives 
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1 f' 

■ = - I X coanxdx = 0 , 

^ —T 

K ^ r • ^ 2 

E)n = “ I X sinnx ax = cos twr. 

X J-x n 


Hence, 

X = 2[( — K cos x) sin a: + cos 2x) sin 2a: 

+ (—3^ cos 3x) sin 3a: + • • • ] 


= 2 ^sir 


sin 2a: . sin 3a: 
2 3 


In this particular case, only the sine terms remain. It may be 
noted that whenever the function /(a:) is an odd function, that is, 
when /(—a:) = —f(x), then an = 0, for w = 0, 1, 2 • • • , since, 
for such a function, 

J(x) cos nxdx= fix) COS nx dx. 

Similarly, if /(a:) is an even function, that is, when/(— a:) = /(a:) 
then hn = 0, forn = 1, 2, 3, • • • , since 

f{x) sin nx dx = --Jj fix) sin nx dx, 


so that the function would be represented by a series of cosine 
terms. 

If in the foregoing illustration the first four terms be plotted by 
composition of 


y = 2 sin a:, y = — sin 2x, y — % sin 3a:, y = sin 4a:, 
the curve 

y =5= 2 sin a: — sin 2a: + % sin 3a: — sin 4a: 

is obtained. It is represented on Fig. 10. As the number of 
terms is increased, the approximating curves approach y = a: as a 
limit for all values of x, —x < x <ir, but not for a: = ±x. 
Since the series has period 2x, it represents the discontinuous 
function shown in Fig. 11 by a series of parallel lines. It should 
be noted that each term of the series is continuous and the func- 
tion from which the series was derived is continuous, but the 
function represented by the series has finite discontinuities at 
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a? = ± (2fc + l)7r. At such points the series converges to zero, 
which is one-half the value of the sum of the right- and left-hand 
limits. 



Illustrative Example 2. Develop f{x) in Fourier series in the 
interval (— tt, tt), if 

fix) = 0, for — IT < X < 0, 

= V, for 0 < a; < T. 

Now 


ao = ^ ^ 0 • daj + J* ir = t, 

1 r*" 

On == - I IT cos nx dx = 0, 

’T Jo 
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1 r* 1 

6„ = -J irsinna:da: = -(l — cos nv). 

The series is then 

jT , /sin a; , sin 3x , sin 5x , \ 

The graph of f(x) from —r to v consists of the a:-axis from 
— T to 0, and the line AB from 0 to ir (see Fig. 12). There is a 
finite discontinuity for x = 0. For x = 0 the series reduces 



to v/2, which is equal to half the sum of lim /(O — «) and 

€-♦0 

lim /(O + e). It may be obsersred from the series that every 

e — >0 

approximation curve will pass through the point (0, ir/2). The 
figure shows the first, second, third, and fourth approximation 
curves, whose equations are • 

IT ^ I o • '^ \ c%( • I sin 3a;\ 

2 / = 2 ' = 2 / = 2 + ^ ^ 3 ~ / 

y = -+2\mix + -^ + —^)> 
as well as the graph of /(x). 

At a; = ±ir the series reduces to t/2, and again every approxi- 
mation curve gives this same value for the ordinate at ±ir. 
This value is one-half the sum of /(~"7r+) and/(7r— ). 

Illustrative Example 3, Let /(x) be defined by the relations . 

/(X) = —IT, if — TT < X < 0, 

= X, if 0 < X < tt; 
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then the Fourier coefficients for f{x) are given by 


o„ = ^ J* fix) cos na: da; = ^ ^ J* —v cos nx dx 


-TT dx + 


X'**) 


cos nx dx 


+ I X cos nxdx , 


-ij: 


_1/ 

T \ ny 

__ 1 / cos mr — l \ 

“ v\ J’ 

1/ fo 

fix) sin nxdx ~ ~ I J "*■ 


+ 1 a; sin nx dx 


7r\n n 


cos nr cos nr 

n ) 


= - (1 — 2 cos utt). 
n 


Therefore, 


V TT 2 2 cos 3a: 2 cos bx 

M- 

, « . sin 2a: , 3 sin 3a: sin 4x , 3 sin 5x 

+ 3 sin a, _ + 

When X = 0, the series reduces to 


_^_?fi + i + i + 

.4 IT \1* ^ 3* ^ 6* ^ 


which must coincide with (see Fig. 13) 


Thus, 


/(0+) +/(0-) ^ _T 

2 2 


. 1 . 1. . . . ^ ^ 

4 tVi* / 2 
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j* -r 32 -r gj T 


Also, for X — ±T, the series gives 

— - + + — + — + • • • ^ 


since 


= 0 , 


/("»•+) +/(ir — ) 
2 


= 0 . 


§!• 


This example suggests the use of Fourier series in evaluating 
sums of series of constants. 



1. Show that 


PROBLEMS 


2. If 


then 

/(x) 

3. If 

then 


x*=|’ + 4 2 (-1)”. 


n — 1 


COS nx 


(— TT ^ Tt). 



for — TT < a? < 0, 
for 0 < X <irt 


TT 2 ^ cos ( 2 n — l)a^ V ^ ^ 

4 ^ 2/ (2n - 1)2 2/ n 

n-l n«l 


fix) = 0 
= sin X 


for 

for 0 ^ w, 


1 _ 2 A cos 2nx . 
“ 5f ir 2/ 4n* — 1 

n — l 


1 . 

^ sm X. 


fix) 
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_«**■ — 1 /I -A cos nx .A n sin nx\ 

V V2 Z 1 + ” 2/ 1 + n> A 

n -1 n*«l ' 


6. Deduce from Prob. 1 that 


6. Show that 


cos ax = 


n*l 


Sin ira 


Xj ^ ‘ n’* “ 12‘ 

, A . . . 2q: sin Tra 

n»l 


if —TT ^ TT. 

7. Deduce from Prob. 6 that 



8. Deduce from the expansion of f{x) = a; + a;* in Fourier series in 
the interval ( —tt, tt) that 

2^ 6* 

n=»l 

9. Expand x sin a; and x cos a; in Fourier series ip the interval (0, 2ir). 

10. Find the Fourier series expansion for/(x), if ^ 

fix) =1 for -TT < X < 

= 0 for ^ < X < IT. 

20. Sine and Cosine Series. The Fourier expansion for f(x) 
in (— TT, tt) has the form (19-1), in which the coefficients an and 6« 
are given by (19-2) and (19-3). As previously observed (Sec. 
19), if fix) is an even function, (19-1) reduces to a series contain- 
ing only cosine terms; and if fix) is an odd function, (19-1) 
reduces to a series containing only sine terms. Now suppose 
that it is desired that fix) be expanded in a Fourier series which 
will be used for the interval 0 to ir only. In that case, it is 
frequently convenient to obtain the expansion in terms of sines 
alone or in terms of cosines alone. For this purpose, define 

Fix) s fix) fovO < X <T 
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F(x) s /(—a;) for — tt < o; < 0, 

so that F{x) is an even function identical with/(x) inO < x < t. 
For an even function: 


1 T' 

bn = ~ I F{x) sin nx dx 

if r f' 1 

= - J F(a:) sin nx dx + F{x) sin nx dx\ 

= i ”” Jo { — nx){ — dx) + F(x) sinna:da:j 


F(a:) sin nxdx | F(x) sin nx dx 


-if - r 

IT L JO 


F{x) sin nx dx I F{x) sin nx dx 


]=«, 


-i[r 


F{x) cos nx dx 


Fix) cos nx dx + I F(x) cos nx dx | 
»■ Jo J 

Fix) cos nx dx + f Fix) cos nx dxl 


= - I F(x) cos nx dx. 

^ Jo 

Hence, in the expansion of Fix) in the interval — tt < x < tt, only 
the cosine terms appear. Moreover, F(x) is identical with /(x) 
for 0 < X < TT. Therefore,* 

(20-1) / (x) = ^ + ai cos X + az cos 2x + • • • + an cos nx + • • • 
in the interval (0, t), where 


(20-2) 


2 r 

an = - I 
TT Jo 


/(x) cos nx dx. 


Similarly, if F(x) be defined so that 
Fix) s fix) 


for 0 < X < TT 


F(x) s ^fi^x) for — TT < X < 0, 

* If fix) has a finite discontinuity at the point x »» xo; then the left-hand 
member of (20-1) is defined to be HlfixQ+) +/(xo-")I. 
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then the an all vanish and 

(20-3) / {x) = 6i sin a: + 62 sin 2a; -H • • • + 6n sin wa; + • • • , 
where 

(20-4) 6n = “ I f{x) sin nx dx. 

^ Jo 

Thus, /(a;) can be represented in the interval 0 < a; < t by either 
(20-1) or (20-3). Frequently, one series is more desirable than 
the other. 

Example. As has been determined already (see Illustrative Example 
1, Sec. 19), the expansion for f{x) = a; in a sine series is 

. / . sin 2a; , sin 3a; \ 

® = 2 Vsm a; - - ■ ). 

This series represents /(a;) = a; in the interval (— x, ir). If one is inter- 



ested in the values of the function in the interval (0, tt), the same func- 
tion can be expanded in a series of cosines. 

In fact, in the interval (0, tt), 

IT 4 ^ cos X , cos 3a; , cos 5x , \ 

® = 2 “ iF vnr + +•••;. 

since 

^=§Jo ard® =ir 
and 

2 f » 2 

Z In xcosnx dx = [(-1)’^ - 1]. 

TT Tl*7r ^ 

The sine series represents the odd function shown in Fig. 14 and the 
cosine series the even function in Fig. 15. The two graphs are identical 
in the interval (0, t). 

PROBLEMS 

1. Show that if c is a constant, then, in 0 < a; < ir, 

4y^.„ ^ , sin Sx , sin 5a; 
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2. Give sine and cosine developments of y = a; sin x in the interval 
(0,ir). 

3. Show that, in (0,ir), 




I [(t “ p) - ^sin 2a: + (^ - sin 3x 


— sin 4x + 



4. Prove that, if /(x) is any function of x, it can be expressed as the 
sum of an even function of x and an odd function of x. 

6. Show that, if /(x) = x for 0 < x < t/2 and /(x) = — x for 

v/2 < X < TT, then 


. IT 2 /cos 2x , cos 6x . cos lOx . \ 

/(») = 4-jl,-ir- + -35- + — p-+ • • • ;• 

6. Show that 

1., (2 sin I). -5)52^ ilO<».<x, 

7. Find the expansion in the series of sines, if 


fix) = ja:, 


for 0 ^ X 


= i (it - X), 


for 2 ^ ^ 


8, Expand /(x) = e» in the series of cosines in the interval (0,7r). 

21. Extension of Interval of Expansion. The methods devel- 
oped up to this point restrict the interval in which /(x) can be 
expanded in a Fourier series to (— ir, ir). In many problems, 
it is desirable to develop f{x) in a Fourier series that will be 
valid over a wider interval. In order to obtain an expansion 
that will hold for the interval (— Z, 0, change the variable by 

replacing x by - 2 . Then /(x) = / 
w 

Fourier series .in 2 , 




can be developed in a 


(2M) /(i,) 

in which 



>, On COS m + ^ bn 
nta n-1 


COS run dz 


sin ne, 
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6n 



sin nz dz. 


The expression (21-1) will be valid for — ir < z < ir; but 
z = irx/l so that (21-1) becomes 


(21-2) fix) = ^ + 2 ^ 

ns“l n«l 

Also, 

a„ = i cos nzdz = ^ J ^ f(x) cos da; 

and 

, 1 P . , 1 .. V . mrx , 

On = ~ I / 1—1 sin nzdz = I f{x) sin — y- aa;. 

Example. Develop f{x) in Fourier series in the interval ( — 2, 2), 
if f(x) = 0 for —2 < (c < 0 and f{x) = 1 for 0 < a; < 2. Here 


do 


dn 

6n 



) = ^' 




2 . nirx 1 \ 

1 • cos dX J 

1 . nTTX 

= — sin — 

’ rnr 2 

2 

o = «’ 


1 • sin ox j 

= — (1 — cos rwr). 
rnr 


Therefore, 

V 1,2/. Tra; , 1 . STra? 

/(*) = o + z(sin-jr + oSm- 7 r 


7r\ 


, 1 . bwx , 
+ 5““T- + 



PROBLEMS 

1. The expansion of S{x) is desired for 0 < a; < Z. If F(x) = fix) 
for 0 < a; < Z and Fix) = —fi—x) for — Z < a? < 0 [that is. Fix) is 
defined as an odd function], show that the expansion of Fix) and fix) 
for 0 < a; < Z is 


where 

r 2 N . nirx j 

“ iX ““ T 

If (pix) s fix) for 0 < a; < Z and ^(x) = /(— x) for — Z < x < 0 [that 
is, vj(x) is defined as an even function], show that the expansion of ^(x) 
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&ndfix) for 0 < « < Z is 

do I nirx 

2 + “■ 

where 

2 A/ V nirx j 
an = j I /(x) cos -j- dx. 

2. Using the results of the preceding problem, obtain the sine and 
cosine expansions of the following functions: 


(a) f{x) = 1 in the interval (0, 2) ; 

(b) f(x) = a; in the interval (0, 1) ; 

(c) f ix) = x^ in the interval (0, 3). 


3 . Expand f(x) = costx in the interval (—1, 1). 

4. Expand 

fix) ^ H - X, if 0 < a; < H, 

= a: - if K < X < 1, 

in the series of sines. 

6. Find the expansion in the series of cosines, if 

fix) =0, if 0 < X < 1, 

= 1, if 1< X < 2. 


6. Expand fix) = jx] in the series of cosines in the interval (—1, 1). 

7. Show that the series 



2n7rx 

I 


represents }4l x when 0 < x < Z. 

8. Find the expansion in the series of cosines, if 

fix) = 1 when 0 < x < tt, 

= 0 when x < x < 2x. 


22. Complex Form of Fourier Series. The Fourier series 


(22-1) /(a:) = ^ ^ 

n — 1 . 

with 

On = i J* /(O cos ntdJl, ~ ^ J ^ 

can be written, with the aid of the Euler formula* 
(22-2) = cos « -f- x sin «, 

* See See. 73. 
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(22-3) 


/(«) = X 


where the coefficients c» are defined by the equation 


(22-4) 




dt. 


The index of summation n in (22-3) runs through the set of all 
positive and negative integral values including zero. 

The equivalence of (22-3) and (22-1) can be established in 
the following manner: Substituting from (22-2) in (22-4) gives, 
for n > 0, 


f(t) cos nt dt 


1 

Cn = ^ I f(t) (cos nt — i sin nt) dt 

1 f’" 

= ^ J /(O cos ntdt — 

dn • 

~ 2 

A similar calculation forn < 0 gives 


_ Qn I hn 
C^n - 2 ^ 2 ' 




sin nt dt 


cio 

Co - Y 


Now (22-3) can be written in the form 

00 00 

/(x) = Co + 2) + X 

n = l 

Making use of the expressions for the Cn just found gives 


f(^) = t + 2 + 2 


2 ^ 


®o , c*"* + e“’"* . , 

= 2 + 2 — 2 * 2 


gin* 

2 
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Recalling that 

e*“ + = 2 cos M and e‘“ — e“*“ = Bi sin u 

gives 

00 

/W = ^ + 2 (®» CO® nx + bn sin nx), 

n *1 

which establishes the identity of (22-3) with (22-1). 

PROBLEM 

Show that the Fourier series in the interval {-‘1,1) can be written in 
the form 

n=* + * inirx 

fix) = ^ c„e * 

n a — w 

where 

1 M inirt 

J -I 

23. Differentiation and Integration of Fourier Series. Some 
general rules concerning differentiation and integration of infinite 
series were given in Sec. 8. It may be added here that, if the 
Fourier series represents a function f{x), then the term-by-term 
integral of the series will converge to the integral of f{x). Thus, 
in Sec. 20 it was found that the function x can be expanded in a 
cosine series in (0, ir) to give 

. V IT 4 / cos X , cos 3x , cos bx , \ 

(23-1) • j. 


The term-by-term integral of this series gives 



On account of the presence of the term Trx/2 in the right mem- 
ber of thm equation, the resulting series is not a Fourier series. 
However, if the sine series development for x (Illustrative Example 
1, Sec. 19) be substituted in this term and the like terms collected, 
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the resulting series will be the development of x^/2 in (0, «-) in a 
sine series. 

The series resulting from term-by-term differentiation of a 
Fourier series converges more slowly than the original series, 
and, in fact, it may diverge. Thus, term-by-term differentiation 
of 

rt / . sin 2a: , sin 3a: \ 

x = 2[smx- - ... J 

gives the series 

2(cos X — cos 2x + cos 3x — • • * )> 

which is divergent, as was observed in Sec. 8. Accordingly, great 
caution must be exercised in differentiating Fourier series 
termwise.* 

24. Orthogonal Functions. A set of continuous functions 
(24-1) Ui{x)j U2{x), • • • , Un{x), • • • , 

which do not vanish identically in the interval a x ^ b, is 
said to be orthogonal with respect to the interval (a, h) if the 
functions Ui{x) satisfy the relations 

(24-2) Ui{x)Uj{x) dx = 0, if i j. 

For i = j, (24-2) becomes 

[wi(x)]2 dx s c?, 


where c\ certainly is not zero. 

If each of the orthogonal functions Ui{x) be divided by Ct, 
there will be obtained a system of normal orthogonal functions, 


v,{x) = V2(x) = • • • , 

Cl C2 



that are characterized by the property that 

(24-3) Vi(x)vj(x) dx = 0, if i j, 

= 1, if i = j. 

Consider a set of normal orthogonal functions v»(x), and 
assume that an arbitrary function f(x) can be expanded in a 
series 

* Some important theorems in regard to this will be found in I. S. Sokol- 
nikoff, Advanced Calculus, Sec. 108 . 
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(24-4) /(*) = aiViix) + OiViix) +•••-}- o«t;«(a:) -)-••• 

oa 

which can be integrated term by term. Multiplying both sides 
of (24-4) by Vj(x) and integrating term by term between the 
limits a and b yield 

£/(x)vy(x)dx = a, £ Vi(x)vi(x) dx, 

which, by virtue of (24-3), gives the formula 

(24-6) ai = f(x)vi(x) dxy (i = 1, 2, 3, • • • )• 

The numbers a** are known as the Fourier coefficients of the 
function/(a:) associated with the system of normal and orthogonal 
functions 

Vl(x), V 2 (x)y • • • , Vnix), • • • . 

The set of functions 

1 cos x sin X cos 2x sin 2x cos nx sin nx 

— — ; — ; — * * * , -pr^} * * * 

v27r y/TT y/T \/w y/w y/'K y/w 

is obviously a normal orthogonal set in the interval (— tt, tt), 
and the development of a function f{x) with the aid of this 
particular set of orthogonal functions is precisely the Fourier 
development of f(x). Among many other useful sets of orthog- 
onal functions are the functions of Bessel and Legendre, which 
are of frequent occurrence in applied mathematics and will be 
used in the discussion of some important problems in Secs. 101, 
113, and 114. 



CHAPTER III 


SOLUTION OF EQUATIONS 

Students of engineering, physics, chemistry, and other sciences 
meet the problem of the solution of equations at every stage of 
their work. This chapter gives a brief outline of some of the 
algebraic, graphical, and numerical methods of obtaining the 
real roots of equations with real coefficients, of types that occur 
frequently in the applied sciences. It also contains a short 
summary of those parts of the theory of determinants and the 
theory of matrices that are immediately applicable to the solution 
of systems of linear equations. 

26. Graphical Solutions. The subject of the solution of equa- 
tions will be introduced by considering a simple problem that 
any engineer may be called upon to solve. 

It is required to design a hollow cast-iron sphere, 1 in. in 
thickness, that will just float in water. It is assumed that the 
air in the cavity is completely exhausted. The specific gravity 
of cast iron will be denoted by p, for convenience. 

By the law of Archimedes, the weight of the sphere must 
equal the weight of the displaced water. This gives the con- 
dition on the radius of the sphere, namely, 

= %irp[x^ — (x - 1)®]. 

Simplifying gives 

(26-1) x^ - 3px2 + 3px - p = 0. 

It will be convenient to remove the second-degree term in (26-1). 
To accomplish this, let x = j/ + fc, giving 

y* + dV*k + Zyk* + - 3p(y* + 2yk + A;*) + 3p(y + A;) - p = 0, 

or 

y® + (3A? — 3p)y* + (3 Aj* — 6pk + 3p)y + A;® — 3pAj* + Zpk — p » 0. 

Choosing A; = p makes the equation reduce to 

(26-2) y* -h (3p - 3p*)y - 2p» -h 3p* - p = 0. 

83 
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For cast iron, p — 7.6, and (25-2) becomes 

(25-3) - 146.26^ - 682.5 = 0. 

If (25-3) is solved, the solution of (25-1) is also determined, since 
= y + 7.5. 

A graphical method of solution will be used. The solution 
of (26-3) is equivalent to the simultaneous solution of the system 


(26-4) 


2 = 

z = 146.26y + 682.6. 


The accompanying figure (Fig. 16) represents the graphs of the 
two functions of (25-4) ; since they inter- 
sect at y = 14.0, this value gives an 
approximate solution of (25-3). The cor- 
responding solution of (25-1) is a; = 21.5. 
From the graph, it is clear that there is 
only one real solution of (25-4) and hence 
of (25-3). 

This graphical method can be applied to 
any cubic equation. The general fourth- 
degree equation (quartic) can also be reduced to a form that is 
convenient for graphical methods of solution. 

Consider the quartic 



+ ax^ + bx^ + cx + d = 0. 

Let X — y + k, as in the cubic equation. This substitution 
gives 

y4 ^ y3(4^ + a) -h y2(6fc2 + 3aA; + b) 

+ y{Ak^ -f- 3ajfc2 -j- 2bk + c) + + ak^ + bk^ + ck + d 0. 

In order to remove the term in y®, choose fc = — This reduces 
the equation to the form 

y^ -f Ay* + By + C = 0. 

If A > 0 , the further transformation y = y/A z is made, and the 
equation is reduced to 

AV + AV + B VA« + C = 0, 


5^ + «* + + ff = 0. 


or 
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The solutions of this equation are the same as the solutions of 
the simultaneous system 

u = 

u = —'pz — q. 

The graphs of these two functions are easily plotted, and the 
solutions can be read from the graph. In case A < 0, the 
transformation would hey = \/ — A 2 , which leads to the equation 

_ ^2 _|_ ^ ^ — 0 

and the graphical solution of the system 

u = z^ — 2 ^, 
u = —pz — q. 

This method of solution for the real roots of an equation is 
also applicable to many transcendental equations. In order to 
solve 

Ax — B sin X = 0, 

write it as 

ax — sin x = 0, 

and plot the curves of the simultaneous system 

y = sin X, 
y = ax. 

Similarly, the equation 

a* — = 0 

can be solved graphically by plotting the curves of the equivalent 
simultaneous system 

y = a*, 

y = x^. 

PROBLEMS 

1. Solve graphically 

(a) 2* ~ = 0, 

\h\ - X - I ^ 0, 

(c) X® — re — 0.5 = 0, 
id) e* + X = 0. 

2. Find, graphically, the root of 

tan X — X = 0 


nearest 
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26. Algebraic Solution of Cubic. The graphical method of solution 
is perfectly general, but its accuracy depends upon the accurate con- 
struction of the graphs of the equations in the simultaneous systems. 
This is often extremely laborious and, at most, yields only an approxi- 
mate value of the roots. 

In the case of the linear equation aa; + 6 = 0, where a 5^ 0, the solu- 
tion is X = —h/a. For the quadratic equation aa;^ + 6a; + c = 0, 


where a 0, there are two solutions given bya; = 


— 6 ± \/6^ — 4ac 
2a 


The question naturally arises as to the possibility of obtaining expres- 
sions for the roots of algebraic equations of degree higher than 2. 
This section will be devoted to a derivation of the solutions of the 
general cubic equation 


a^x^ -1- dxX^ “h a^x -j- eta == 0, 


flo 5^ 0. 


Dividing through by ao gives 
(26-1) x^ “h hx^ “i- ca; "h d = 0, 

and the a;* term can be removed by making the change of variable 


The resulting equation is 
(26-2) 


6 

a; = 2/-3 

+ vy + 


where 


and 


p = c 


3 


q = d — 


3 27’ 


In order to solve (26-2), assume that 
(26-3) 2/ = A + 5, 

so that 

2 /» = A« -f- 53 + ZAB{A + B). 

Substitute in this last equation for A -f jB, from (26-3), and there is 
obtained the equation 

(26A) 2/» - SABy - (A» -h B^) = 0. 

A comparison of (26-4) with (26-2) shows that 

ZAB = —p and A* -f- JB* = — g, 
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(26-6) A»B» “ - § and A* + B* = -q. 

If is eliminated by substituting from the second of Eqs. (26-5) into 
the first, there appears the quadratic equation in A\ 

(A*)* + - g = 0, 

♦vhose roots are 



The solution for yields precisely the same values. However, in 
order to satisfy Eq. (26-5), choose* 




-q + yji 


, 2 .. 4 ?“ 


?“ + 


27 


(26-6) 


= 




If the values of y are to be determined from (26-3), it is necessary to 
find the cube roots of and Recall that if x® = a®, then the solu- 
tions for X are given by a, wa, and co^a, where w = — ^ -f i and 


0)2 = — 2 ^ ^ complex roots of unity. Hence, if one cube 

root of A® be denoted by a and one cube root of R® by the cube roots 
of A® are 


whereas those of R® are 

<aa, 

and 

Ci)*a, 


up, 

and 



It would appear that there are nine choices for y, but it should be 
remembered that the values must be paired so that SAB = — p. The 
only pairs that satisfy this condition are a and jS, (oa and co^jS, and o)^ 
and cujS. Hence, the values of y are 

(26-7) 2/1 = Of -h jS, 2/2 = coa + 0)®j8, yz = + o)j9, 

where 



* The opposite choice for the values of A® and R* simply interchanges 
their role in what follows. 
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The solutions of (26-1) can be obtained from the values given in 
(26-7) by recalling that x = y — 6/3. 

The expressions for a and are quite complicated, and when the 
quantity under the square-root sign has a negative value the values of 
a and cannot, in general, be determined. This is the so-called 
irreducible case of the cubic, which can, however, be solved by using a 
trigonometric method. This method will be 
^ described later in the section, but first it is im- 

I portant to find a criterion that will determine 

/ in advance which method should be used. 

/ In order to determine the character of the 

y ^ roots of (26-2), whose coefficients are assumed 

^ to be real, consider the function 

^ fiy) + py + q 

Fig. 17. and its maximum and minimum values. Since 

riy) = 32/2 + p, 


it appears that, if p > 0, then f'{y) is always positive and f(y) is an 
increasing function. In this case the graph of f {y) has the form shown 
in Fig. 17, and there is evidently only one real value for which /( 2 /) = 0. 

If p < 0, however, f{y) is zero when 2/ = ± V-p/S. Since 
f\y) = 62/, it follows that y ^ P /3 gives a minimum value to 

/(y), whereas y = — V—p/3 furnishes a maximum value. The cor- 


responding values of f{y) are 

fl+gp"' 

'/-i 

and 



J ^ 

ff-gP^ 

^ 3’ 


The graph of fiy) will have the 
appearance of one of the curves in 
Fig. 18. 

It is evident that f{y) = 0 will 
have only one real root if the graph 
of f(y) has the appearance shown by (1) or (5), that is, if the maximum 
and minimum values of f(y) are of the same sign. Hence, 

(« + ipV^)(9-fpV- |)>0, 

or 

8 ’ + > 0 , 
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is the condition that (26-2) have only one real root. It may be observed 
that this condition is automatically satisfied if p ^ 0. It should be 
noted that, if p = 0, Eq. (26-2) reduces to + 3 = 0, which obviously 
has only one real root. 

If (26-2) has three real and distinct roots, then the graph of f(y) must 
have the appearance shown in (3), and it follows that the maximum 
and minimum values must be of opposite sign. Hence, 

+ H7 P^<0 

is the condition for three real and unequal roots. 

If O'* + either the maximum or the minimum value of 

f{y) must be zero [see (2) and (4)], and (26-2) will have three real roots, 
of which two will be equal (a so-called double root). 

The expression 

(26-8) A = -27g2 - 4p3 

is called the discriminant of the cubic equation (26-2), for its value 
determines the character of the roots of the equation. The discriminant 
for (26-1), obtained by replacing p and q in (26-8) by their values in 
terms of 6, c, and d, is 

(26-9) A = ISbcd - - 4c« - 27^^. 

It may be worth noting that the discriminant of any algebraic equa- 
tion, with leading coefficient unity, is the product of the squares of 
the differences of the roots taken two at a time. Inasmuch as 

ixi - X2y{^2 - XiYixz - XiY = (yi - 2/2)H2/2 - ys^iys - 2/l)^ 

the discriminant has the same value for (26-1) and (26-2). 

In view of the definition of A, it follows that 

if A < 0, one root is real and two are complex; 
if A = 0, all the roots are real and two are equal; 
if A > 0, the three roots are real and unequal. 

Example, Consider the cubic equation 

+ 3x2 + 9a; - 1 q. 

From (26-9), it follows that A = —2592, and hence there will be one real 
root and two complex roots. Setting x = 2 / — 1 yields the reduced cubic 

2/^ 4- 6i/ — 8 =* 0, 

and substituting p — 6 and q = —8 in (26-6) gives A* = 4 -h 2 \/6 and 
s 4 — 2 \/6. Therefore, the solutions for y are 

•^4 +2V6 + - 2 V6, a aJ^4 +2\/6 + «» "^4 - 2 V6, 

and w* ^^4 + 2 Vc + « "^4-2 VS. 
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The solutions of the original equation can now be obtained by recalling that 

« « y - 1. 

The discussion of the solution of the cubic equation will be con- 
cluded by giving the derivation of the expressions for the roots in the 
case where the roots are real and unequal (that is, when A s — 27^^ 
- 4p* > 0). 

Let 


— 2 = r cos 0 


and 


Then* 


and 


a = (r cos 0 + ir sin 0)^ = ^cos ^ 4- i sin ^ 


P = (r cos 0 — ir sin 0)^ = ^cos | ^ sin 

If it is noted that 


and 


CO = cos ~ + t sin ^ 
o o 


< 0 * = cos ^ — i sin —t 
V 3 


it is easily checked that the expressions for 

yi = a + ^, 2/2 = coa + co^/S, 2/3 s co^a + cojS 

become 

(26-10) = 2rM cos t y, = 2rH cos i±^, 

o 3 


Since 


2/8 = 2rH cos 

3 


and 



the values of yi, j/j, and y% can be obtained directly from the coeffi- 
cients of (26-2) or (26-1). 

* By De Moivre’s theorem (cos S + i sin #)* = cos n9 + * sin »9. 
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Example, Determine the real roots of 

~ 3x* + 3 =» 0. 

Here 

A 4(-27)(3) - 27(9) > 0, 

and the roots are all real and unequal. Since p = — 3 and 5 = 1, it follows 
that r = 1 and cos ^ — K- Hence, 


and 

rt 2w ^ Stt 14ir 

2 / 1=2 cos -y » 2 / 2=2 cos -g-» 2/8 = 2 cos —• 

The solutions of the general quartic equation 

+ hx^ + cx^ + dx + e = 0 

can be found, but the methods of obtaining the expressions for the 
roots depend upon the solution of an auxiliary cubic equation. More- 
over, these expressions are, in general, so involved that they are prac- 
tically useless for computation.* It has been shown that the ordinary 
operations of algebra are, in general, insufficient for the purpose of 
obtaining exact solutions of algebraic equations of degree higher than 4. 
However, it is possible to obtain the expressions for the solutions of the 
general equation of the fifth degree with the aid of elliptic integrals. 

The reader should not confuse the problem of obtaining expressions 
for the exact solutions of the general algebraic equation with that of 
calculating numerical approximations to the roots of specific equations 
which have numerical coefficients. The latter problem will be dis- 
cussed in Secs. 28 and 29, and it will be shown that the real roots of 
such equations can be computed to any desired degree of accuracy. 
Moreover, if the roots are rational they can always be determined 
exactly. 

PROBLEMS 

Determine the roots of the following equations: 

(а) 2 /® - 22 / - 1 = 0 ; 

(б) 2 /" - 146.252/ - 682.5 = 0; 

(c) a;® — — 5a; — 3 = 0; 

(d) X® — 2x^ — X + 2 = 0; 

( 6 ) X® - 6x2 + 6 x - 2 = 0; 

(/) X® + 6x2 + 3x + 18 = 0; 

(g) 2x® -h 3x2 + 3x + 1 == q. 

* See Dickson, L. E., First Course in Theory of Equations, pp. 5Q-54; 
Burnside, W. S., and A. W. Panton, Theory of Equations, vol. 1, pp. 121-142. 
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27. Some Algebraic Theorems. The student of any applied 
science is usually interested in obtaining numerical values, 
correct to a certain number of decimal places, for the roots of 
equations. Unless the roots are rational, the expressions for the 
exact roots, provided that they can be found at all, are usually 
complicated and the process of determining numerical values 
from them is tedious. Accordingly, it is distinctly useful to 
consider other methods of finding these numerical values. 
Horner^s method, Newton^s method, and the method of inter- 
polation are the ones most frequently used; they will be dis- 
cussed in Secs. 28 and 29. However, all these methods are 
based on the assumption that a root has first been isolated, that 
is, that there have been determined two values of the variable 
such that between them lies one and only one root. In many 
practical problems the physical setup is a guide in this isolation 
process. This section contains a review of some theorems* from 
the theory of equations that provide preliminary information as 
to the character and location of the roots. 

Theorem 1. (Fundamental Theorem of Algebra.) Every 
algebraic equation 

J{x) = + • * * + dn^ix + an = 0 

has a root 

It should be noted that this theorem does not hold for non- 
algebraic equations. For example, the equation e® = 0 has 
no root. 

Theorem 2. (Remainder Theorem.) If the polynomial 

f{x) s aox^ + aix^-^ + • • • + an-ix + On 

is divided by x — b until the remainder is independent of Xj then 
this remainder hxis the value f{b)T 

Theorem 3. (Factor Theorem.) If f{b) = 0, then x — b is 
a factor of the polynomial f{x) and b is a root of f{x) = 0. 

This theorem follows directly from Theorem 2. In many 
cases the easiest way to compute the value of /(b) is to perform 
the division of f{x) by x — 6. This is a particularly useful 

* Those students who are not already familiar with these theorems and 
their proofs will benefit by referring to H. B. Fine, College Algebra, pp. 
425-453, and L. E. Dickson, First Course in the Theory of Equations, Chap. 

II. 
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method when the factor theorem is being used for the purpose of 
determining the roots of f{x) = 0. For if a; — 6 is a factor of 
f(x), it follows that f{x) = (x — b) g(x), where g{x) is a poly- 
nomial of degree one less than that of f{x). Obviously the roots 
of g{x) = 0 will be the remaining roots of f{x) = 0, so that only 
g{x) =0 need be considered in attempting to find these roots. 
Moreover, when f{x) is divided hy x — b the quotient is g(x). 
If synthetic division is used, the computation is usually quite 
simple. 

Example. If f{x) = x^ + 2x'^ + 2x + 1 is divided by a; + 1, the 
quotient is a;^ + a; + 1 and the remainder is zero. Hence, a; = —1 is 
a root of f{x) = 0 and the remaining roots are determined by solving 

+ a; + 1 = 0. 

Theorem 4. Every algebraic equation of degree n has exactly 
n roots if a root of multiplicity m is counted as m roots. 

A root b of f(x) = 0 is said to be a root of multiplicity m ii 
(x — 6)"* is a factor of f(x) but (a: — 6)”*+^ is not a factor of f{x). 

It follows from Theorems 3 and 4 that the polynomial of 
degree n can be factored into n linear factors, so that 

f(x) = aox^ + aix^-'^ + • • * + an-ix + an 
= ao(a; — a:i)(a; — x^ • • • {x — Xn)- 

Theorem 6. 'If 

f(x) = aox^ + aia;'*-^ + • • • + an-ix + an 

has integral coefficients and if f(x) = 0 has the rational root b/c^ 
where b and c are integers without a common divisor, then b is an 
exact divisor of an and c is an exact divisor of ao. 

Example. Consider the equation 

f{x) s 2a;® + X* + aj — 1 = 0. 

The only possible rational roots are ± 1 and + Since /(I) = 3, 
/(—I) = —3, /(— K) = and/(K) = 0, it follows that H is the 
only rational root. As a matter of fact, if f{x) is divided by x — M the 
quotient is 2x® + 2x + 2 whose factors are 2, x — w, and x — 
where co and co® are the complex roots of unity.* 

Theorem 6. Given f{x) s x’* + aix^^ + • • • + an^ix + an 
= 0. If f{p) andf{b) are of opposite sign, then there exists at least 

* See Sec. 26 and the example following Theorem 9 of this section. 
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one root of f{x) == 0 between a and b. Moreover , the number of such 
roots is odd. 

Graphically this means that y = f{x) must cross the x-axis 
an odd number of times between a and 6. 

Example. If f(x) = 8a;® — 12a;® — 2a; + 3 =* 0, 

/(-I) - -15, m = 3; /(I) = -3, /(2) = 15. 

Since /(—I) is negative and /(O) is positive, there is at least one root 
between —1 and 0. Similarly, there is a root between 0 and 1, and 
another between 1 and 2. 

Theorem 7. (Descartes^ Rule of Signs.) The number of 
positive real roots of an algebraic equation f{x) =0 with real coeffi- 
cients is either equal to the number of variations in sign of f{x) or less 
than that number by a positive even integer. The number of negative 
real roots of /(a;) = 0 is either equal to the number of variations in 
sign of f{’-x) or less than that number by a positive even integer. 

Example. f{x) = 8a;® — 12a:® — 2a; + 3 has two changes in sign, and 
therefore there are either two or no positive roots of f(x) = 0. Also, 
/(—x) s —8a;® — 12a;® + 2a; + 3 has only one change in sign, and f{x) 
must have one negative root. 

Theorem 8. Every algebraic equation of odd degree, with real 
coefficients, and leading coefficient positive has at least one real root 
whose sign is opposite to that of the constant term. 

Example. Since f{x) = 8a;® — 12a;® — 2a? + 3 = 0 is of odd degree 
and the constant term is positive, it follows that there must be at least 
one negative root. 

Theorem 9. If an algebraic equation f{x) = 0 with real coeffi- 
cients has a root a + bi, where b 9 ^ 0 , and a and b are real, it also has 
the root a — bi. 

Example. Thus, a;® — 1 = 0 has the root — M + H V3 i, and there- 
fore it has the root — K — K t. This theorem states that imaginary 
roots always occur in pairs. 

PROBLEMS 

1. Mnd all the roots of the following equations: 

(a) a;® + 2a;® — 4a; — 8 *= 0; 

(5) 2x« - a;® - 5a; - 2 » 0; 

(c) 4a;^ + 4a;® + 3a;® — a; — 1 = 0; 

(d) 2a;® - 3a;® - 3a; - 2 « 0. 
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2 . Isolate the roots of the following equations between consecutive 
integers: 

(а) a:® — 2a;^ — a; + 1 = 0; 

(б) 2x5 + 4x2 - 2x - 3 = 0; 

(c) x* + 5x2 + 6x + 1 = 0; 

{d) x^ - 5x2 + 3 = 0. 


28. Homer’s Method. Many readers are already familiar with 
Horner’s method of determining the value, to any desired number of 
decimal places, of the real roots of algebraic equations. However, the 
development given here is somewhat different from that used in the 
texts on algebra, in that it depends on Taylor’s series expansion. 
Suppose that the equation is 

(28-1) fix) = aox^ + aix"“^ + • • • + an^ix + «„ = 0 

and that it is known that the equation has a root between c and c + 1 , 
where c is an integer. If/(x) is expanded in Taylor’s series in powers of 
X — c, there will result* a polynomial in x — c, namely, 


fix) =/(c)+/'(c)(x-c) + 


ric) 

2! 


(x - c)2 + 


in - 1)! 


(x — c)^“^ + 


n\ 


(x — c)»‘. 


Now, let X — c = Xi and 


mo) 

r\ 


— i4n— r. 


Then (28-1) is replaced by 


(28-2) /i(xi) = An + An-\X\ + • • • + AiXi” ' + AoXi" = 0. 

Since (28-1) had a root between c and c + 1 and since Xi = x — c, 
it is evident that (28-2) has a root between 0 and 1. By the use of 
Theorem 6, Sec. 27, this root can be isolated between d and d -f- 0.1, 
where d has the form a/10 and 0 < a < 9. Moreover, /i(xi) = /(xi + c ) ; 
and it follows that, if f\ has a root between d and d -|- 0.1, then / has a 
root between c + d and c -|- d + 0.1. It should be noted that c may 
be negative but that d will always be positive or zero. 

The function /i(xi) can be expanded in Taylor’s series in powers of 
xi — d; and, if X 2 = Xi — d, there will be obtained an equation 


/2(X2) = Bn + Bn-.\X2 + • * ’ + 5iX2’*'"^ + B^X^^ = 0. 

But /i(xi) = 0 had a root between d and d + 0.1 ; and since X 2 = Xi — d, 
it follows that ftix^ - 0 will have a root between 0 and 0.1. 

This process can be continued as long as desired, each step deter- 
mining another decimal place of the root of the original equation (28-1). 

* Since fix) is a polynomial of the nth degree, the derivatives of order 
higher than n are all zero. 
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The solution of a specific equation may help to clarify the procedure. 
Let it be required to find the values of the real roots of the equation 

Fix) s - 3a;* - 6® - 3 « 0. 

Since there is only one variation in sign, Fix) has at most one positive 
root. Fi—x) has three variations, and so there will be at most three 
negative roots. The only possibilities for rational roots are ±1 and 
±3. Since — 1) = 0, it follows that a; = —1 is a root. Moreover, 
if Fix) is divided by a; + 1, the quotient is/(a;) = a;® — 3a; — 3. Hence, 
the remaining roots of Fix) =0 are the three roots of 


fix) = a;® — 3a; — 3 = 0. 


It is easily checked that fix) — 0 has no rational roots. Also, 
A = 108 — 243, so that]there is only one real root which, since /(2) = — 1 
and /(3) = 15, must lie between 2 and 3. Therefore, fix) will be 
expanded in powers of a; — 2. Since 


fix) = X® - 3x - 3, 
fix) = 3x* - 3, 
f\x) = 6x, 
rix) = 6 , 


/(2) = -1, 
/'(2) = 9, 
/"(2) = 12, 
/"'( 2 ) = 6 , 


the expansion becomes 

/(x) = -1 + 9(x - 2) + ^ (x - 2)* + |j (X - 2)». 


Replacing x — 2 by Xi gives 

fiixi) = -1 + 9xi + 6xi* + xi® = 0. 

Since the real root of this equation lies between 0 and 1, the Xi* and Xi® 
terms do not contribute very much to the value of /i(xi). Hence, a 
first approximation to the root can be obtained by setting 9xi — 1 = 0. 
This gives xi = K = 0.111 • • • , and suggests that the root probably 
lies between 0.1 and 0.2. It is easy to show that/i(0.1) = —0.039 and 
/i(0.2) = 1.048; there is thus a root between 0.1 and 0.2, and it is 
evidently closer to 0.1. Therefore, fix) = 0 has a root between 2.1 
and 2.2. 

Expanding /i(xi) in powers of xi — 0.1 gives 

/i(x,) = -0.039 + 10.23(xi - 0.1) + ^ (xx - 0.1)» + | (xi - 0.1)^ 

and replacing x\ — 0.1 by X 2 yields 

hix%) = -0.039 + 10.23X2 + 6.3X2* + Xa® = 0. 

Now 10.23x2 — 0.039 = 0 gives the approximation X 2 = 0.0038, and 
testing 0.003 and 0.004 reveals that /2(0.003) * —0.008253273 and 
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/2(0.004) = +0.002020864. Thus, the root lies between 0.003 and 0.004 
and is closer to 0.004. If it is desired to determine the root of /(x) = 0 
to three decimal places only, this value will be 2.104. If more decimal 
places are desired, the process can be continued. It should be noted 
that in each succeeding step the terms of the second and third degree 
contribute less, so that the linear approximation becomes better. 

PROBLEMS 


1. Apply Horner’s method to find the cube root of 25, correct to three 
decimal places. 

2. Determine the real roots of a;® — 2a: — 1 = 0 by Horner’s method. 

3. Determine the root of a:^ + a:* — — a: + 5 = 0, which lies 

between 2 and 3. 

4. Determine the real root of 2x^ — 30:^ + a: — 1 = 0. 

5. Determine the roots of x^ — Sx^ + 3 = 0. 

6. Find, correct to three decimal places, the value of the root of 
a;® + 3a;® — 2x^ + a; + 1 =0, which lies between — 1 and 0. 

7. A sphere 2 ft. in diameter is formed of wood whose specific gravity 
is %. Find to three significant figures the depth h to which the sphere 

will sink in water. j^The volume of a spherical segment ('-!>] 

The volume of the submerged segment is equal to the volume of the 
displaced water, which must weigh as much as the sphere. Since water 
weighs 62.5 lb. per cubic foot. 


and, since r = 1, 


i'-i) 


62.5 = I JIT’ I 62.5 


h* - 3h^ + % = 0 . 


29. Newton’s Method. Horner’s method of obtaining a 
numerical solution of an equation is probably the most useful 
scheme for solving algebraic equations, but 
it is not applicable to trigonometric, ex- 
ponential, or logarithmic equations. A 
method applicable to these types as well as 
to algebraic equations was developed by Sir 
Isaac Newton sometime before 1676. 

Newton applied his method to an alge- 
braic equation, but it will be introduced 
here in the solution of a problem involving 
a trigonometric function. 

Let it be required to find the angle subtended at the center of 
a circle by an arc whose length is double the length of its chord 
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(Fig. 19). Let the arc BCA be an arc of length 2BA, Let 2x 
be the angle (measured in radians) subtended at the center of the 
circle. Then, arc BCA = 2xr and BA = 2 DA = 2r sin x. If 
arc BCA = 2BA, then 2xr = 4r sin x, or a: — 2 sin a: = 0. 

The graphical solution of equations 
of this type was discussed in Sec. 25. 
A first approximation can be obtained 
by graphical means. If y = x and 
y = 2 sin X are plotted, it appears 
from the graph (Fig. 20)- that they 
intersect for a; lying between 108° and 
109°, or, expressing this in radians, 

1.8850 < X < 1.9024. 



If xi = 1.8850 be chosen as the first approximation, the question 
of improving this value will be discussed first from the following 
graphical considerations. 

If the part of the curve y = a; — 2 sin a: in the vicinity of the 
root be drawn on a large scale, it will have the appearance shown 
in Fig. 21. It is clear from the graph that adding to Xi the 



distance AE, cut off by the tangent line to the curve at xi = 
1.8850, will give a value X 2 which is a better approximation to the 
actual root xo. But AE is the subtangent at Xi and is equal to 

— where f{x) == a; — 2 sin x. Thus,* 


X2 = Xi 


f(^i) 

r(x^y 


* See, in this connection, Prob. 8, at the end of this section. 



$29 SOLUTION OF EQUATIONS 99 


Similarly, upon using 0:2 as the second approximation and 
observing that — is the subtangent EF^ the third approxi- 

J \X2) 

mation is found to be 

r - r - 

* f{Xi)’ 

and in general the nth approximation Xn is given by 

f(Xn-l) 


(29-1) 




Since Xi = 1.885, the formula gives 


X2 = Xi’ — 


f(xi) 

T(xi) 

In a similar way, 


xi - 


Xi — 2 sin xi 
1 — 2 cos Xi 


)• 


, 1.8850 - 1.9022 , 

= 1.8850 - 1 ^ 0 6180 = 


Xz = X2 


fM _ 1 1*8956 - 2 sin 1.8956 

f(x 2 j ~ ^ 1 - 2 cos 1.8956 


1.8955. 


It follows that the angle subtended by the arc is 3.7910 radians. 

The use of Newton^s method requires some preliminary 
examination of the equation. It may happen that the equation 
is of such a character that the second approximation to Xo will be 
worse than the first. A careful examination of the following 
sketches of four types of functions, ■ 
sketched in the vicinity of their 
roots, reveals the fact that some 
care must be exercised in applying 
Newton^s method. For all four 
figures, it is assumed that xo has 
been isolated between Xi and x[. 

The graphical interpretation of the 

correction — K as the subtan- 
f i^i) ^ ^ 

gent must be kept in mind throughout this discussion. If Xi is 
used as the first approximation, then xz will be obtained as the 
second approximation by using Newton^s method; if is used, 
then X 2 will be obtained. 

In Fig. 22, both X 2 and xi are closer to xo than Xi or x[. In this 
case the method would work regardless of which value is chosen 
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as the first approximation. In Fig. 23, Z 2 is better than the 
first approximation xi, but xj is worse than xj. It appears from 



Fia. 23. 


the figure that this occurs because the curve is concave down 
between xi and x(, and hence /"(x) < 0, whereas f(xi) < 0 and 



f(x[) >0. A similar situa- 
tion would obtain if the curve 
is concave up, so that /"(x) 
> 0 (Fig. 24). The reader 
will readily convince himself 
from an inspection of Fig. 23 
that caution must be ex- 
ercised in the choice of the 
first approximation if the 
curve has a maximum (or a 


minimum) in the vicinity of Xo. 

If the curve has the appearance indicated in Fig. 25, then it is 
evident that the choice of either Xi or x'l as the first approximation 
will yield a second approxima- 
tion which is worse than the 
first one. This is due to the 
fact that the curve has a point 
of inflection between xi and 

x'l. 

From the foregoing discus- 
sion, it is apparent that New- 
ton's method should not be 
applied before making an investigation of the behavior of the first 
and second derivatives of /(x) in the vicinity of the root. The 
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conclusions drawn from this discussion can be summarized in the 
following practical rule for determining the choice of the first 
approximation: If f(x) andf'{x) do not vanish in the given interval 
(x\y x[) and if the signs of f(xi) and f(x[) are opposite j then the first 
approximation should he chosen as that one of the two end points for 
which f(x) andf\x) have the same sign. 

It can be proved* that if the single-valued continuous function 
f(x) is of such a nature that/(a:) 

= 0 has only one real root in 
(xi, x[) and both /'(a;) andf'(x) 
are continuous and do not 
vanish in (xi, x[)y then repeated 
applications of Newton’s 
method will determine the value 
of the root of f(x) = 0 to any 
desired number of decimal 
places. 

The cases to which Newton’s 
method does not apply can be 
treated by a method of interpolation {regula falsi) that is appli- 
cable to any equation. 

Let X be the value of x for which the chord AB intersects the 
x-axis. From similar triangles (Fig. 26), 

X ~ Xi _ x[ — X 
-f{xi) ■“ /(x'l) ’ 

Solving for x gives 

_ xifjxj) - x;/(xi) 

The value x is clearly a better approximation than either Xi 
or x'l. 

PROBLEMS 

1. Solve Prob. 7, See. 28, by Newton’s method. Also, apply the 
method of interpolation. 

2. Determine the angle subtended at the center of a circle by a chord 
which cuts off a segment v/hose area is one-quarter of that of the circle. 

3. Find the roots of e* — 4x = 0, correct to four decimal places. 

4. Solve X — cos X = 0. 

* See Weber, H., Algebra, 2d ed. vol. 1, pp. 380-382; Coate, G. T., On the 
Convergence of Newton's Method of Approximation, Amer. Math. Monthlyj 
vol. 44, pp. 464-466, 1937. 
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6 . Solve X = tan x in the vicinity ot x - 

6 . Solve X + €* = 0 . 

7. Solve — X — 1 =* 0 . 

8 . Show with the aid of Taylor^s series that, if x = xi is an approxi- 
mate root of /(x) = 0 , then the nth approximation is, in general, 
determined from the formula (29-1). 

Hint: fix) - fix,) +f'ix,)ix - x,) + ■ ■ ■ - Xi)" + • ' • : 

and if /(X 2 ) ~ 0 , then 

0 = /(xi) +f\xA)(x 2 - Xi). 

30. Determinants of the Second and Third Order. The solu- 
tion of systems of linear equations involves the determination of 
the particular values of two or more variables that will satisfy 
simultaneously a set of equations in those variables. Since the 
discussion is simplified by using certain properties of deter- 
minants and matrices, the remainder of this chapter is devoted 
to some elementary theory of determinants and matrices and its 
application to the solution of systems of linear equations. 

Consider first a system composed of two linear equations in 
two unknowns, namely, 

(30-1) I °‘*t 

' ) a 2 X + 622/ = /C2. 

If y is eliminated between these two equations, there is obtained 
the equation 

(30-2) (aib2 — a2bi)x = kib2 — ^ 2 ^ 1 ; 

and if X is eliminated, there results 

(30-3) (0162 — (i2bi)y = aik2 — a2ki. 

If the expression ai 62 — 0261 is not zero, the two equations 
(30-2) and (30-3) can be solved to give values for x and y. That 
the values so obtained are actually the solutions of the system 
(30-1) can be verified by substitution in Eqs. (30-1). 

The expression aib 2 — a 2 bi occurs as the coefficient for both 
X and y. Denote it by the symbol 

(30-4) <1162 — ct2bi = ^ ^ • 

(I 2 02 

This symbol is called a determinant of the second order. It is 
also called the determinant of the coefficients of the system 



§30 


SOLUTION OF EQUATIONS 


103 


(30-1), for the elements 6f its first column are the coeflScients of 
X and the elements of its second column are the coefficients of y. 
Using this notation, (30-2) and (30-3) become 


(30-5) 


0>l 61 


ki bi 


ai hi 


Oi hi 

0,2 62 

X = 

ki 62 

f 

02 62 

y = 

O2 ^2 


The definition (30-4) provides the method of evaluating the 
symbol. If 


D = 


ai hi 
(l2 h^ 




the unique solution of (30-1) can be written as 


h\ 61 

I 1 

loi hi 

A/2 62 

rk 

II 

I O2 ^2 

T\ 


If D = 0, aihi = ajbi or aija^ = hijh^^ But if the correspond- 
ing coefficients of the two equations are proportional, the two lines, 
whose equations are given by (30-1), are parallel (if aija^ 5*^ 
^1/^2) or coincident (if aija^ = 61/62 = ^1/^2). In the first 
case, the determinants appearing as the right-hand members of 
the equations in (30-5) will be different from zero and there will 
be no solution for a; and y. In the second case, these deter- 
minants, as well as D, are zero and any pair of values a;, y that 
satisfies one equation of the system will satisfy the other equa- 
tion, also. 


Example 1. For the system 


2x - 32/ = -4 2 -3 

3a; — 2 / = 1> 3—1 


-2 -h 9 = 



7, 


Example 2. For the system 


but 


2a; — 32/ = 4 
6a; — 92/ = 5, 


D = 


2 -3 
6 -9 


= 0 , 


2 

6 




The two lines whose equations are given are parallel. 
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Example 3. For the system 

2x - .3y = 4 
6 a; - 92/ = 12, ^ ~ 

2-3 4 

6 “ -9 “ 12* 

The two lines are coincident. 

Consider next the system of three linear equations in three 
unknowns, 

{ aix + biy + ciz = fci, 
a2X + h2y + C 2 Z = k2, ^ 

azx + hzy + C32 = hz. 

If these equations are multiplied, respectively, by 

62C3 — hzC 2 i bzCi — biCzj biC 2 — 62^1, 

and the resulting equations are added, the sum is 

(30-7) ((Z162C3 — aibzC 2 + ci 2 bzCi — (i 2 biCz + (izbiC 2 — azb 2 Ci)x 

== kibzCz — kibzCz "1" kJbzCx — kzhiCz ~h kjbiCz — k^bzCi, 

The coefficient of x in (30-7) can be denoted by the symbol 

CLl bi Cl 

(30-8) D ^ fl 2 C 2 — (Z 162 C 3 — ciJbzC2 “|- (iJ)zCi — ajbiCz 

az bz Cz “h azbiCz — CLzbzCi. 

This symbol is called a determinant of the third order. It is 
also the determinant of the coefficients of the system (30-6). 
Using the notation of (30-8), Eq. (30-7) can be written as 

di hi Cl ki bi Cl 

^)x — ~ d2 ^2 C2 X kz bz Cz 
Qz bz Cz kz bz Cz 

Similarly it can be shown that 

di hi Cl di fci Cl 

d2 h2 C2 y — d2 kz Cz , 

ds bz Cz ds kz C3 

di hi Cl di hi Aji 

d 2 h 2 C2 Z — d 2 h 2 kz 
d3 hs Cz ds hs kz 




and 
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If D 7 *^ 0, the unique solutions for y, and z can be obtained as 


(30-9) 


X = 


ki bi Cl 

Jc2 ^2 

kz bz Cz 


D 


y = 


cbi ki Cl 
0/2 ^2 ^2 
€Lz kz Cz 


D 


z = 


Q/i hi ki 
0/2 ^2 
OLz bz kz 


D 


In order to show that the values of x, y, and z, given in (30-9), 
actually satisfy Eqs. (30-6), these values can be substituted in 
the given equations. 

If D = 0, the three equations (30-6) are either inconsistent or 
dependent. A detailed analytic discussion of these cases will be 
given in Sec. 35. Since the three equations of (30-6) are the 
equations of three planes, a geometrical interpretation will now 
be given. 

If the three equations are inconsistent, the tjiree planes are 
all parallel, or two are parallel and are cut by the third plane in 
two parallel lines. In either case, there is obviously no solution 
for Xj y, and z. If the equations are dependent, all three planes 
intersect in the same line or all three planes coincide. In either 
case there will be an infinite number of solutions for x, y, and z. 


Example, For the system 


Therefore, 


X = 


3a; - 

y - 

- z = 2, 


3 

-1 -1 

X — 

2y - 

- 3z = 0, 

D = 

1 

-2 -3 

Ax + 

y + 2z = 4. 


4 

1 2 

2 -1 

-1 



3 

2 -1 

0 -2 

-3 



1 

0 -3 

4 1 

2 

9 

y = 

4 

4 2 

9 



9 


= 2 . 




2 , 


1. Evaluate 


PROBLEMS 


1 2 3 


2 

0 -3 


4 -2 1 

2-1 4 

3 -1 -2| 

f 

1 

-1 

4 2 

1 -2 

, and 

6 0-1 
2 3-3 


2. Find the solutions of the following systems of equations by using 
determinants: 
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(а) 6aj — 42/ ~ 3, 

2x + Sy^ 7; 

(б) 2aJ + 32/ — 2« =s 4, 

jc + 2/ - 0 = 2, 

3a; — 62/ + 3« = 0; 

(c) 3a; - 22/ « 7, 

3y + 22 = 6, 

2a; + 32 = 1; 

(d) 3a; + 22/ + 22 « 3, 

a; — 42/ + 22 = 4, 
2a; + 2/ + « = 2. 


31. Determinants of the nth Order. Determinants of the 
second and third orders were defined in the preceding section. 
These are merely special cases of the definition of the determinants 
of any order n. Instead of a symbol with 2^ or 3^ elements, the 
determinant of the nth order is defined as the symbol, with n 
rows and n columns, 


D ^ 


dll di2 
d2l U22 


din 

d2n 


\dnldn2 • • • dnn\ 


which stands for the sum* of the n! terms ( — iydk^iak ^2 * • • 
where ki, ^ 2 , • * • , fcn are the numbers 1, 2, • • • , n in some 
order. The integer k is defined as the number of inversions of 
order of the subscripts ki, k 2 , • • • , kn from the normal order 
1, 2, • • • , n, where a particular arrangement is said to have k 
inversions of order if it is necessary to make k successive inter- 
changes of adjacent elementsf in order to make the arrangement 
assume the normal order. There are n! terms since there are n! 
permutations of the n first subscripts. Moreover, it is evident 
that each term contains as a factor one and only one element from 
each row and one and only one element from each column. 

* This sum is sometimes called the expansion of the determinant. 

t It should be noted that it is not necessary to specify that the inter- 
changes should be of adjacent elements, for it can be proved that, if any 
particular arrangement can be obtained by k interchanges of adjacent 
elements and also by*fc' interchanges of some other type, then k and A;' are 
always either both even or both odd. Hence, the sign of the term is inde- 
pendent of the particular succession of interchanges. 
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Example. Consider the third-order determinant 

ail ai2 ai8 

D = ^21 a22 a23 • 

aai ^32 a33 

The six terms of the expansion are, apart from sign, 

aiia22a33, aiia32a23, a2iai2a33, 
a2ia32ai3, a3iai2a25 , a3ia22ai3. 

The first term, in which the first subscripts have the normal order, is 
called the diagonal term, and its sign is positive. In the second term 
the arrangement 132 requires the interchange of 2 and 3 to make it 
assume the normal order; therefore. A; == 1, and the term has a negative 
sign. Similarly, the third term has a negative sign. The fourth term 
will have a positive sign, for the arrangement 231 requires the inter- 
change of 3 and 1 followed by the interchange of 2 and 1 in order to 
assume the normal order. Similarly, it appears that the fifth term will 
have a positive sign. In the sixth term, it is necessary to make three 
interchanges (3 and 2, 3 and 1, and 2 and 1) in order to arrive at the 
normal order; hence, this term will have a negative sign. As a result of 
this investigation, it follows that 

D = aiia22a33 aiia32a23 a2iai2ais d" aziaz^iz d* azia-i^zz azia^zaiz 

It is evident that if k is equal to zero or an even number the 
term will have a positive sign, whereas if fc is odd the term will 
be negative. 

PROBLEM 

Find the signs of the six terms involving an in the expansion of the 
determinant 

an ai2 an au 
a^i a22 a23 a24 
aai a82 ass as4 
a4i a42 a43 a44 

32. Properties of Determinants. 1. The value of a determi- 
nant is not changed if in the symbol the elements of corresponding 
rows and columns are interchanged. 


If 

D s 

Oil ai2 • 

©21 a22 * 

• • din 

* * d2n 



\anl dni • 

dnn 
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then the determinant formed by interchanging the corresponding 
rows and columns is 



Oil O21 * 

* ’ Onl 

D ' s 

012 022 * 

* • dnZ 


ain 02n * 

• * Onn 


Any term • • * CLk^n of Z>, where ki^ k 2 f • • • ,kn are 

the numbers 1, 2, * • • , n in some order, will correspond to a 
term ( — * • * dk^n of Z)', for each determinant must 
contain every possible term that is a product of one and only 
one element from each row and each column. But the number 
of inversions is the same for the term of D as it is for the term of 
D', owing to the fact that the corresponding first subscripts are 
the same. It follows that each term of D occurs also in D\ and 
conversely each term of D' occurs also in D. 


Example. If 


then 



2. An interchange of any two rows or of any two columns of a 
determinant will merely change the sign of the determinant. 

If D is the original determinant and D" is the determinant 
having the fth andjth rows of D interchanged, then the expansion 
of D" will have the first subscripts of each term the same as those 
of the corresponding term of D, except that i and j will be inter- 
changed. Since it requires one interchange to restore i and j to 
their original order in each term, the sign of every term will be 
changed. Thus, D" = — D. 


Example. If 




2 5 
1 -1 
-3 -2 


3 

4 
1 


- 66 , 




2 6 3 

-3 -2 1 

1 -1 4 


66 . 


then 
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3. If any two rows or any two columns of a determinant ate 
identical^ the value of the determinant is zero. 

For, by property 2, if these two rows (or columns) were inter- 
changed, the sign of D should be changed. But since these two 
rows (or columns) are identical, D remains unchanged. There- 
fore, D =?= — D, and hence D = 0. 

Example. If 


then 



-1 2 

4 3 

5 -2 


D - 0. 


4. If each element of any row or any column he multiplied by 
m, the value of the determinant is multiplied by m. 

This follows from the definition of the determinant. Since one 
and only one element of any row or column occurs in each term, 
each term will be multiplied by m and therefore the value of the 
determinant is multiplied by m. 


Example 1. If 


and 



3 

4 
1 


= -66 




2 5 
1 -1 
-6 -4 


3 

4 
2 




which has each element of the last row twice the corresponding element 
of the last row of D, then 


S = -132 and D = 2D. 


Example 2. If 


then 




3 

2 

4 


1 

2 

. 4 

2 

9 

2 

1 

= 2-3 

3 

2 

1 


-6 

3 

-1 


-2 

3 

-1 


5 . From properties 3 and 4 , it follows that the value of a deter-- 
minant is zero if any two rows or any two columns have corre- 
sponding elements proportional. 
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6. The product of two determinants D and Z)', both of order n, 
is the nth-order determinant D" which has as the element in its 
ith row and jth column the sum 


n 

5) ^hbkj = + ai2?)2j + • * • + Clinbnij 


which is formed by multiplying each element a^ of the ith row of D 
by the corresponding element bkj of the jth column of D\ 

Thus, if 


then 


D = 


an ai2 
U21 U22 


and 


D' = 


bn bi 2 

b21 &22 


D D' ^ D" 


anbn + «i2?>2i 
^21611 + U22i>21 


cinbi2 + ^12622 
a>2lbl2 “h ^22^^22 


Example. The product of the following determinants is easily found 
by expanding the product determinant: 


sin X cos X 1 


sin X sec x esc a? 


sec X tan x 1 

• 

cos X — tan x — cot x 


CSC a; cot a; 1 


- 1 - 1 - 1 


0 

tan a? — sin a: — 1 

— cos x cot x 

tan a; + sin a; 

- 

1 0 

sec X CSC X ~ 2 

cos X cot X 

sec a; CSC X — 2 

0 


= 2 cos* x(2 — sec x esc x). 


33. Minors. The method of evaluating a determinant by the 
use of the definition of Sec. 31 is exceedingly tedious, especially 
if n ^ 4. There are other schemes for this evaluation, and these 
require the definition of the minors of a determinant. The 
simplest of these schemes will be described and used here. 

If, in the determinant D, the ith row and the jih column be 
suppressed, the resulting determinant Aij (which is of order one 
less than the order of D) is called the minor of the element at,-, 
which is in the ith row and jth column. 


ail O12 ttis 014 




021 U22 O 28 024 
Osi 032 OsS U84 
041 C 42 ^48 au 


Example, If 
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A2S == 


flu ai2 ai4 
Ctai ^32 (Z 34 . 


|a41 042 O 44 I 


From the definition of a determinant, it is evident that an A a is 
composed of all the terms of D which contain the element ai, as a 
factor, except for the possibility that all the signs may be reversed. 
Then the expression ( — ly^anAn is composed of all the terms of 
D containing an as a factor; ( — l)**a 2 iA 2 i is composed of all the 
terms containing 021 as a factor; (— I)**a3ii4.3i is composed of all 
the terms containing a 3 i as a factor; etc. But D is composed of all 
the terms containing an, 021 , nsi, • • • , ani as a factor, and so 


D = { — ly^anAn + ( — l)*'a2iA2i + * • • + ( — l)*"aniAni. 


It can be proved* that fci = 1 4- 1. ^2 = 2 + 1, fcs = 3 + 1, 
• • * , Ain = n + 1, so that 


D = anAn anA^i + * * * + (“ 


In the above development for D the elements an, « 2 i, • * • , flni 
are the elements of the first column of Z). Similarly, the value 
of D can be formed by taking the elements of any other column 
or of any row. 

Using the tth column gives 

D = {^ly^aiiAii + ( — I)*2a2ti4.2» + • • • + ( — l)*”ani 2 lni, 

where fci = i + 1, *2 = i + 2, • • • , = z + n. Similarly, 

using the ith row gives 

D = { — ly^anAii + ( — I)**a* 2 -Ai 2 -}-••• 4 - { — ly^ainAin, 

where ki = i + 1, k2 = i + 2, • • • , kn = i + n. It may be 
observed that each K is equal to the sum of the subscripts of its 
aij and is thus equal to the sum of the number of the row and the 
number of the column in which this clement occurs. This 
development is known as the expansion hy minors, or the simple 
Laplace development. 

Since the term cofactor is frequently used in applications of 
this type of development, it will be defined here. The cofactor 
Cif of an element a^ is defined as the signed minor, that is, 

Cn = 

* Dickson, L. E., First Course in Theory of Equations, pp. 101-127; 
Fine, H. B.,College Algebra, pp. 492-619. 
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Thus, the expression for D can be written as 


D=X 




or as 


Example. 




13 4 0 6 
0 5 2 1 
0 3 4 0 
12 7 1 



5 2 1 


4 0 6 


4 0 6 


4 0 6 

= 3 

3 4 0 
2 7 1 

- 0 

3 4 0 
2 7 1 

0 

+ 

5 2 1 
2 7 1 

- 1 

5 2 1 
3 4 0 


.a[-. 


+ 4 


0 


5 2 
2 7 


2 1 
4 0| 


- 0 


5 1 
3 0| 


+ 6 


5 2 
3 4 


= -9(2 - 7) + 12(5 - 2) - 4(0 - 4) - 6(20 - 6) 
= 13. 


0 


Here, the first expansion is made by using the elements of the first 
column, for it contains two zeros (the third row is an equally good 
choice). The expansion of the first third-order determinant is made by 
using the elements of the second row, but the third column could be 
used to equal advantage. In the expansion of the last third-order 
determinant the first row was chosen, but the third row and the second 
and third columns provide equally good choices. 


The following theorem is given here because of its frequent use 
in many fields of pure and applied mathematics: 

n 

Theorem. The sum S aaCkf is zero, if k 9 ^ i. 

Each term of this sum is formed by taking the product of the 
cofactor of an element of the kth row by the corresponding 
element of the ith row. This is the expansion of a determinant 
whose ith and fcth rows are identical and whose value is accord- 


ingly zero. 
Example. 


n 

Similarly, it follows that S aifiik = 0, if ft 5*^ j. 


Let 




3 -1 2 

1 2 -1 

4 -3 -2 


Then, 


Cu « -7, Cl, - -2, Ci8 = -11 
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Similarly, 



-28 + 6 + 22 


0 . 


3 

02jCiy =* —7 — 4 -f- 11 = 0. 

J-1 


By using the theory of determinants, the solution of a system 
of n non-homogeneous linear equations in n unknowns can be 
obtained. The rule for effecting the solution will be stated but 
not proved. * The proof for the cases when n = 2 and n = 3 has 
already been given in Sec. 30. 

Cramer^ s Rule, Let 


■ CLnX\ + anX2 + • * * + CLlnXn = kly 

(33-1) a2lXi + a22X2 + * * * + Cl2nXn = ^2, 


(^n\X\ “b (ln2.X2 “f" • • • -f" (XnnXn — 

be a system of n equations in the unknowns Xi, X 2 , * • • , Xn, 
such that the determinant 



dll di2 * ■ 

din 

D = 

U 2 I d22 ’ ‘ 

d2n 


dnl dn2 • ' 

dnn 


of the coefficiehts is not zero. The system (33-1) has a unique 
solution given by 

Xi ^ ^ 

where Di is the determinant formed by replacing the elements 
CLiiy a 2 iy aziy • • • , ttni of the zth column of D by ki, ^ 2 , fcs, * • • , 
kn, respectively. 

Example, Solve, by Cramer^s rule, the system 

3x + y + 2z ^ 3j 
2x — — z = —3, 

x + 2y + « = 4. 

* Dickson, L. E., First Course in Theory of Equations, pp. 114-115. 
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3 1 2 

D = 2 -3 -1 = 8 
1 2 1 

3 1 2 3 3 2 

3-3-1 2 -3 -1 

4 2 18 1 4 1 16 ^ 

8 ~ 8 ”^’ 8 

3 1 3 

2 -3 -3 
12 4-8 

^ 8 8 "■ 


PROBLEMS 

1. Evaluate 


1 

-2 

0 

3 


2 

-1 

3 

0 

1 

3 

2 

-1 

-1 

4 

1 

2 


1 

4 

0 

^ , and 
0 

0 

4 

-3 

2 

2 

0 

-1 

-3 

1 

1 

1 

3 

-3 

1 

0 

1 

3 

1 

0 

1 


-1 

0 

2 

0 

1 

2 

0 

-4 


2. Solve, by Cramer’s rule, the systems 


(a) X + 2y + Sz = (b) 2x y Sz “ 2, 

2x — y + z = Qj 3aj — 2!/ — 22 = 1, 

3x + 2/ — 2 = 4. X — y + z ^ —1, 

(c) a; + 2i/ = 1, (d) 2a; + t/ + 32 + it; = -2, 

2a; — 2 / — 22 = 3, 5a; + 3i/ — 2 — w; = 1, 

— a; + 2 / + 32 = 2. a; — 21/ + 42 + Sii; = 4, 

3a; — 2 / + 2 = 2. 

34. Matrices and Linear Dependence. In order to discuss 
the systems arising in the succeeding sections, it is convenient to 
give a short introduction to the theory of matrices. * 

An m X n matrix is defined as a system of mn quantities an 
arranged in a rectangular array of m rows and n columns. If 
m ^ n, the array is called a square matrix of order n. The 
quantities an are called the elements of the matrix. Thus, 

I dll di2 • • • ain 

d21 022 ' * • 02n 

or 


fOml Om2 * * ’ Ow»nll 
* For detailed treatment see M. Bocher, Introduction to Higher Algebra, 
pp. 20-53; L. E. Dickson, Modern Algebraic Theories, pp. 30-^3. 


Oil 

Oi2 

' * Oln' 

021 

022 

* • 02n 



\Oml 

Om2 

Ofnn 1 


Here 

and 

X = 
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where double bars or parentheses are used to enclose the array of 
elements. If the order of the elements in (34-1) is changed or if 
any element is changed, a different matrix results. Any two 
matrices A and B are said to be equal if and only if every element 
of A is equal to the corresponding element of JB, that is, if 
an = bij for every i and j. 

If the matrix is square, it is possible to form from the elements 
of the matrix a determinant whose elements have the same 
arrangement as those of the matrix. The determinant is called 
the determinant of the matrix. From any matrix, other matrices 
can be obtained by striking out any number of rows and columns. 
Certain of these matrices will be square matrices, and the 
determinants of these matrices are called the determinants of the 
matrix. For an m X n matrix, there are square matrices of 
orders 1, 2, • • * , p, where p is equal to the smaller of the 
numbers m and n. 

Example. The 2X3 matrix 

^ _ /an an an\ 

\a 21 ^22 ^ 23 / 

contains the first-order square matrices (an), (U12), (U23), etc., obtained 
by striking out any two columns and any one row. It also contains 
the second-order square matrices 

( an ai2\ /an ai3\ /an ai3\ 

021 0,22/ \021 O23/ \022 O23/ 

obtained by striking out any column of A. 

In many applications, it is useful to employ the notion of the 
rank of a matrix A. This is defined in terms of the determinants 
of A. A matrix A is said to be of rank r if there exists at least 
one r-rowed determinant of A that is not zero, whereas all deter-- 
minants of A of order higher than r are zero.* 

Example. If 

( 1 0 1 3\ 

2 1 0 - 2 , 

-1 -1 1 5 / 

In case an m X n matrix contains no determinants of order higher than r, 
obviously r is the smaller of the numbers m and n, and the matrix is said 
to be of rank r. 
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the third-order determinants are 

1 0 1 

2 1 0*0, 

-1 -1 1 


Since 


there is at least one second-order determinant different from zero, 
whereas all third-order determinants of A are zero. Therefore, the 
rank of A is 2. 

It should be observed that a matrix is said to have rank zero 
if all of its elements are zero. 

The notion of linear dependence is of importance in connection 
with the study of systems of linear equations, and it will be con- 
sidered next. 

A set of m, m 2, quantities fi^ / 2 , fz, * • * ,fm (which may be 
constants or functions of any number of variables) is said to be 
linearly dependent if there exist m constants Ci, 02 , • • • , which 
are not all zero, such that 

(34-2) Clfi + C2/2 + • • • + Crnfm * 0. 

If no such constants exist, the quantities /» are said to be linearly 
independent 

Example, If the /* are the polynomials 

Vj ^ 2x2 ~ 3x2/ + 4z, 

y, ~ + 2x2/ - 3z, 

fz(^i y, «) “ 4x2 ^ xy - 2z, 

and if the constants are chosen as ci * 1, C 2 = 2, cs = —1, then 
Ci/i + cst/i + Czfz = 0. 

Therefore, these three polynomials are linearly dependent. 

It is evident that, whenever the set of quantities is linearly 
dependent, at least one of the /» can be expressed as a linear 
combination of the others. Thus, from (34-2), if ci 9 ^ 0, then 

fl =* ^ 2/2 + + • “ * + Umfmt 
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where 

C2 C3 , 

(I2 — — — f (I3 ) etc. 

Cl Cl 

The definition of linear dependence requires the existence of at 
least one constant c* 9 ^ 0, and therefore the solution for fi is 
assured. 

Obviously, in most cases it would be extremely difficult to 
apply the definition in order to establish the linear dependence 
(or independence) of a given set of quantities. In case the 
quantities /»• are linear functions of n variables, there is a simple 
test which will be stated without proof. * 

Theorem. The m linear functions 

fi = anxi + ai2X2 + • * * + ainXny (f = 1, 2, • • • , m), 

are linearly dependent if and only if the matrix of the coefficients 
is of rank r < m. M or cover ^ there are exactly r of the fi that form 
a linearly independent set. 

If m > n, obviously r < m, and it follows that any set of m 
linear functions in less than m unknowns must be linearly 
dependent. 

The fact that the polynomials 

fi^ 2x - Sy + 42, 

/a = x + 2y - 32, 

/a = 4x + y - 22, 

are linearly dependent can be determined by observing that the matrix 
of the coefficients. 



is of rank 2. 

36. Consistent and Inconsistent Systems of Equations. A set 

of equations that have at least one common solution is said to 
be a consistent set of equations. A set for which there exists no 
common solution is called an inconsistent set. 

The question of consistency is frequently of practical impor- 
tance. For example, in setting up problems in electrical net- 
works, there are often more conditions than there are variables, 

* Dickbon, L. E., Modern Algebraic Theories, pp. 65-60; Bocher, M., 
Introduction to Higher Algebra, pp. 34-38. 
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This leads to a system in which there are more equations than 
there are unknowns. It is important to have a method for 
testing whether all the conditions can be satisfied simultaneously. 

Theorem 1. Consider a system of m linear equations in n 
unknowns, 


(35-1) 


aiiXi + ai2X2 + * • • + ainXn == ki, 
a2lXi + ^22^2 + * * * + a2nXn = AJ2> 


\ afn\Xx ”f” ^m2^2 "f" * * * "I" ^mnXn — km, 


where at least one ki jA 0. If the matrix of the coefficients is of 
rank r, Eqs, (35-1) are consistent provided that the rank of 


K ^ 


'ail 

ai2 * * 

• Oi„ ki \ 

^21 

^22 * * 

a2n ^2 I 


am2 • 

amn kml 


is also r. 

The matrix K is called the augmented matrix. The proof of 
this theorem will be found in any standard work on higher 
algebra. 


Example 1. Consider the system 


Since 


2x + Zy = 1, 
a; - = 4, 

4a; — y — 9. 



is of rank 2, the equations are consistent if 



is also of rank 2. This condition is satisfied; for the determinant of K 
is zero, and there exists a second-order determinant of K that is dif- 
ferent from zero. 

Example 2. The system 

2a; + 3y = 1, 
a; — 2y = 4, 

4a; — 2 / = 6 
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is inconsistent, because 



is of rank 3, whereas the matrix A is of rank 2. 

In the case in which there are n equations in n unknowns, the 
theorem on consistent equations shows that if the determinant of 
A is zero, so tha t the rank of is r < n, then the rank of K must 
be r also, if the set of equations is to be consistent. If the rank of 
K is greater than r, the set of equations is inconsistent. This 
provides the analytic discussion that should accompany the 
geometric discussion given forn = 3 in Sec. 30. 

If the set of equations is consistent and the rank of A is r, 
then it can be shown that n — r of the unknowns can be given 
arbitrary values, and the values of the remaining r unknowns are 
determined uniquely in terms of those n — r arbitrary values. 
These n — r unknowns cannot be chosen at random, for the 
m X r matrix of the coefficients of the remaining r unknowns must 
have rank r if these unknowns are to be uniquely determined. 

Example 3. Solve the system 

X - y A- 2z ^ 3, 

X y - 2z = 1, 

X Sy — Qz = —1. 

Since A and K are both of rank 2, the equations are consistent. If 
either y or z is chosen arbitrarily the matrix of the coefficients of the 
remaining variables will have rank 2. If z = kj the equations to 
be solved are 

X — y = 3 — 2A;, 

X y = I 2k, 
x -f- 32/ = — 1 + 6A;. 

Solving the first two for x and y gives x = 2 and y — 2k — 1, These 
values are seen to satisfy the third equation. Therefore, the solutions 
X 2, y = 2k ^ I, z = k satisfy the original system for all values of k. 

The preceding discussion has dealt with non-homogeneous 
linear equations. In case the ki are all zero, the system becomes 
the set of homogeneous equations 

+ ^12X2 + • * • + CllnXn = 0 , 

(I21X1 + a22X2 + • • ’ + 0,2nXn = 0 , 


OmlXl + am2X2 + ’ * ’ + dmnXn = 0 . 
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Obviously, Xi ^ — • • • = Xn =« 0 is a solution of (35-2). 

It may happen that there are other Solutions. If ai, 02 , • * • , Un 
is a solution of (35-2), it is evident that fcai, A;a 2 , • • • , kant where 
fc is an arbitrary constant, will be a solution, also. The condition 
for solutions different from the Xi — X 2 = • • • = Xn = 0 solu- 
tion will -be stated without proof. 

Theorem 2. The system (35-2) mil have a solution different 
from the solution Xi — • • • = Xn — Oy if the rank of the 

matrix of the coefficients is less than n. 

It follows that if the number of equations is less than the 
number of unknowns, that is, if m < n, there are always solutions 
other than the obvious zero solution. If m = n, there exist 
other solutions if the determinant of the square matrix of the 
coeflBcients is zero. As in the case of the non-homogeneous 
system, if the m X n matrix of the coefficients is of rank r, then 
n — r of the unknowns can be specified arbitrarily and the 
remaining r unknowns will be uniquely determined, provided that 
the rank of the matrix of the remaining unknowns is r. 

Example 4. Consider the system 

2x — y + Sz = Of 
X + 2y — z = Of 
Sx + iy + z — 0, 

Here 

2-1 3 

\A\= I 2 -1 = 10. 

3 4 1 

Therefore, a? == 0, y = 0, 2 = 0 is the only solution. 

Example 5. Consider 

3a; — 2y = 0, 
x + Ay ^ Of 
2a; — y - 0, 

for which the matrix of the coefficients is of rank 2. Since the number 
of unknowns is 2, a; = 0, y =* 0 is the only solution. 

Example 6. Consider 

2a; — y + 32 * 0, 

3a; + 2y -h 2 * 0, 
a; + 3y — 22 = 0, 

5a; + y + 42 « 0. 
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Here, 



which is of rank 2. Since the number of unknowns is 3, the system 
has solutions other than a: = 0, y = 0, 2 = 0. Let z = kj and solve 
any two of the equations for x and y. If the first two are chosen, 
X = — A; and y — k. By substitution, it is easily verified that x ~ — A?, 
y — kf z — k satisfies all four equations for any choice of k. 

Example 7. Consider 

2x — 4i/ + 2 = 0, 

3x + y — 22 = 0. 

For this system, 



which is of rank 2. Since the number of unknowns is greater than the 
number of equations, there exist other solutions. Let 2 = A;, and solve 
the two equations for x and y. There results x = ^k and y = ^ik. 
Thus, X = Hk, y = Hk^ 2 = A; is a solution for any choice of k. 
Example 8. Consider 

X — y + 2 = 0, 

2x + 32/ + 2 = 0, 

3x + 2y + 22 = 0. 

Here, 

1 -1 1 

|A| = 2 3 1=0. 

3 2 2 

Since the determinant of A is zero, there are solutions different from 
X = 0, 2/ = 0, 2 = 0. Let z — k^ and solve any two of the equations. 
If the first two are chosen, x = —%ky y = Hky z — k. It is verifiable 
by substitution that these values satisfy all three equations, whatever 
be the choice of Aj. 

PROBLEMS 

1 . Investigate the following systems and find solutions whenever the 
systems are consistent: 

(a) X - 22/ = 3, (6) 2x + 2/ - « = li 

2x + 2/ = L X - 22/ + 2 = 3, 

3x — 2/ = 4. 4x — 32/ + 2 = 5. 

(c) 3x + 22/ = 4, (d) 2x - 2^ + 32 = 4, 

X — 32/ * 1, X + 2/ - 32 = -1, 

2x + 62/ “ —1. 5x — 2/ + 3z = 7. 
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2. Investigate for consistency, 
they exist. 

(a) X + 3^/ — 2z = Of 
2a; — y + « == 0. 

(c) 3x -2y + 2 = 0, 
jc + 2y — 22 = 0, 

2a; — 2 / + 22 = 0. 

(e) 4x — 22/ + 2 = 0, 

2a; — 2 / + 32 = 0, 

2a; — 2 / — 22 = 0, 

Oa; — 32/ + 42 = 0. 


and obtain non-zero solutions when 

(6) a; - 22 / = 0, 

3a; + 2 / = 0, 

2a; — y — 0, 

(d) 2a; — 42/ + 32 = 0, 

X -h 2y — 2z = 0, 

3x — 2y + 2 = 0. 

(/) a; + 22 / + 2z = 0, 

3a; - y + 2 = 0, 

2a; + 32/ + 22 = 0, 

X + 4y — 2z = 0, 



CHAPTER IV 

PARTIAL DIFFERENTIATION 


36. Ftmctions of Several Variables. Most of the functions 
considered in the preceding chapters depended on a single 
independent variable. This chapter is devoted to a study of 
functions depending on more than one independent variable. 

A simple example of a function of two independent variables 
X and y is 

2 = xy, 

which can be thought to represent the area of a rectangle whose 
sides are x and y. Again, the volume t; of a rectangular parallel- 
epiped whose edges are x^ y, and namely, 

V = xyz, 

is an example of a function of three independent variables x^ 
y, and 2 . A function w of n independent variables Xi, X 2 , • • • , 
Xn can be denoted by 

U = /(xi, X2, • • • , Xn). 

A real function of a single 
independent variable x, say 
y = /(^)j can be represented 
graphically by a curve in the 
, xy-plane. Analogously, a real 
function z = /(x, y), of two 
independent variables x and 
y, can be thought to represent a surface in the three-dimen- 
sional space referred to a set of coordinate axes x, y, z (Fig. 
27). However, one must not become too much dependent 
on geometric interpretations, for such interpretations may 
prove to be of more hindrance than help. For instance, 
the function v = xyZy representing the volume of a rectangular 
parallelepiped, depends on three independent variables x, y^ 

123 
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and 2 and hence cannot be conveniently represented geometrically 
in a space of three dimensions. 

Corresponding to the definition of continuity of a function 
of a single independent variable x (see Sec. 7), it will be said that 
a function z — f{x, y) is continuous at the point (xo, j/o) provided 
that a small change in the values' of x and y produces a small 
change in the value of z. More precisely, if the value of the 
function z — f(x, y) at the point (a:o, yo) is Zo, then the continuity 
of the function at the point (xo, yo) means that* 

(36-1) lim/(x, y) = /(xo, yo) = Zo. 

X— >850 

In writing the left-hand member of (36-1), it is assumed that the 
limit is independent of the mode of approach of (Xj y) to (xq, yo)' 
The statement embodied in (36-1) is another way of saying 
that 

/(a:, y) = fixo, yo) + *, 

where lim € = 0; that is, if the function /(x, y) is continuous at 

X— >xo 

(xo, yo), then its value in the neighborhood of the point {xoy yo) can 
be made to differ from the value at the point (a^o, yo) by as little as 
desired. 

If a function is continuous at all points of some region R 
in the xy-plane, then it is said to be continuous in the region R, 
The definition of continuity of a function of more than two 
independent variables is similar. Thus, the continuity of the 
function u = /(x, y, z) at the point (xo, yo, zo) means that 

lim fix, y, z) = /(xo, yo, Zo), 

X— ♦Xo 
V-*V9 

independently of the way in which (x, y, z) approaches (x©, yo, z^). 

PROBLEM 

Describe the surfaces represented by the following equations: 

(a) X + 2y ~ 3, (5) X — y + « =* 1, (c) x « 2, (d) a; * y, 

(e) 2x - 3y + * 1, (/) x* - y» = 0, (y) y* + « 25, ‘ 

(A) y* =» 2x, (f) X* + y® — lOx « 0, {j) X* + y® + 2® “ 1, 

(Ji) X* + 18* » y, (Z) X* + 2 y* + =* 0 , (?n) x* + y* = 2*, 

♦ For details, see I. 8. Sokolnikoff , Advanced Calculus, Chap. III. 
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37. Partial Derivatives. The analytical definition of the 
derivative of a function y = f(x), of a single variable x, is 


dx 


lira ^ 

Ax -->0 


lim 

Ax-*0 


This derivative can be interpreted geometrically as the slope of 
the curve represented by the equation y == f{x) (Fig. 28). 





jylay. * 






Fia. 28 . 


It is natural to extend the definition of the derivative to 
functions of several variables in the following way: Consider 
the function z = /(x, y) of two independent variables x and y. 
If y is held fast, z becomes a function of the single variable x 
and its derivative with respect to x can be computed in the 
usual way. Let AZx denote the increment in the function 
^ = f(^) y) when y is kept fixed and x is changed by an amount 
Lx] that is, 

= Six + Ax, y) - fix, y). 

Then, 


I* Ar. 
hm — 
A »-.0 Ax 


lim /(^ + Ax, y) - fix, y) 
A *-,0 Ax 


is called the •partial derivative of z •with respect to x and is denoted 
by the symbol dzfdx, or Zx, or/,. 
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Similarly, the partial derivative of z with respect to y is defined 
by 

^ = lim y + Ay) - fix, y) 

In general, if = /(xi, X 2 , • • • , Xn) is a function of n inde- 
pendent variables Xi, X 2 , * • • , Xn, then du/dXi denotes the 
derivative of u with respect to Xi when the remaining variables 
are treated as constants. Thus, if 

z = x^ + x^y + y^y 

then 

= 3x^ + 2xy and = x^ + Sy^. 

dx dy ^ 

Also, if u = sin (ax + by + cz), then 

a cos (ax + by + cz) (both y and z held constant) ; 
b cos (ax + by + cz) (both x and z held constant) ; 
c cos (ax + by + cz) (both x and y held constant). 

In the case of z = /(x, y), it is easy to provide a simple geo- 
metric interpretation of partial derivatives (Fig. 27). The equa- 
tion z = /(x, y) is the equation of a surface; and if x is given the 
fixed value xo, z == /(xo, y) is the equation of the curve AB on 
the surface, formed by the intersection of the surface and the 
plane x = xo. Then dz/dy gives the value of the slope at any 
point of AB. Similarly, if y is given the constant value yo, 
then z = /(x, yo) is the equation of the curve CD on the surface, 
and dz/dx gives the slope at any point of CD. 

PROBLEMS 

1. Find dz/dx and dz/dy for each of the following functions: 

(a) z = y/x; (6) z = x®y + tan-^ {y/x)\ (c) z = sin xy + x; 

(d) z * 6* log y; (e) z = x^y + sin“^ x. 

2. Find du/dXy dw/dy, and du/dz for each of the following functions: 

(a) w * x*y + yz — xz*; (6) u = xyz + log xy; 

(c) w » z sin""^ {x/y)] (d) = (x* + y* + z*)^; 

(e) w « (x* + y* + z*)“^. 


du 

dx 

— 

^y 

du 

dz 
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38. Total Differential. In the case of a function of one 
variable, y = f(x), the derivative of y with respect to a: is defined 
as 


lim ~ 
A*-,o Aa: 


/'(*), 


so that Ay/ Ax = f'(x) + e, where lim * = 0. Therefore, 

Ax— >0 

f(x + — f(x) ^ Ay = f\x) Ax + € Ax, 

where e is an infinitesimal which vanishes with Ax. Then, 
fix) Ax = fix) dx 

is defined as the differential dy. 

For the independent variable x, the terms increment^' and 
differential” are synonymous (that is. Ax = dx). However, 
it should be noted that the differential dy (of the dependent 
variable y) and the increment Ay differ by an amount e Ax (see 
Fig. 28). 

The differential of a function of several independent variables 
is defined similarly. Let z = fix, y), and let x and y acquire the 
respective increments Ax and Ay, Then, 

Az = fix + Ax, y + Ay) - /(x, y). 

If 2 = fix, y) is a continuous function, then, as Ax and Ay 
approach zero in any manner, Az also approaches zero as a limit. 
It will be assumed here that fix, y) is continuous and that 
df/dx and df/dy are also continuous. 

The expression for Az can be put in a more useful form by 
adding and subtracting the term/(x, y + Ay), Then, 

= fix + Ax, 2/ + ^y) - fix, t/ + A2/) + fix, 2/ + Ay) ~ /(x, y). 
But 

lim + Ax, y + Ay) - fjx, y + Ay) ^ dfjx, y + Ay) 

^,*—.0 Ax dx 


so that 

f(x + Ax,y + Ay) - f{x, y + Ay) = + *i] 

where lim «i = 0. Moreover, 

Ax-»0 


lim y + Ay) ^ dfjx, y) 
A]/— ♦o ^x ^ dx 
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since the derivative is continuous. Therefore, 


y(a:» y + Ay) ^ df{x, y) 


dx 


dx 


«*> 


where lim cj = 0. 
In like manner, 


f(x, y + Ay) - fix, y) = + «' Ay, 


where lim c' = 0. It follows that 

AV-*0 


Az = Ax + Ay + « Ax + «' Ay, 

ax ay 

in which e = €i + € 2 . 

The expression 

Ax + ~ Ai/ = ^ dx + ^ di/ 

5x dy ^ dx dy 

is defined as the differential of z and denoted by dz. In 
general, if w = /(xi, X2, • * • ,Xn), the total differential is given by 


(38-1) 


du = 


dxi 


df ^ 
5^2 


dxi + ^ dx2 


+ 




The expression for the total differential is called the principal part 
of the increment Aw, and is a close approximation to Aw for 
sufficiently small values of dxi, dx 2 , * • • , and dxn. As in the 
case of a function of a single independent variable, the differential 
of each independent variable is identical with the increment of 
that variable, but the differential of the dependent variable 
differs from the increment. 

If all of the variables except one, say Xi, are considered as 
constants, the resulting differential is called the partial differ- 
ential and is denoted by 


dxiU 


= ^ dxi. 
dXi 


The partial differential expresses, approximately, the change 
in u due to a change Ax< = dXi in the independent variable x,-. 
On the other hand, the total differential du expresses, approxi- 
mately, the change in w due to changes dxi, dx 2 , • * • , dXn in all 
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the independent variables xi, X 2 , * * • , Xn. It may be noted 
that the total differential is equal to the sum of the partial 
differentials. Physically, this corresponds to the principle of 
superposition of effects. When a number of changes are taking 
place simultaneously in any system, each one proceeds as if it 
were independent of the others and the total change is the sum 
of the effects due to the independent changes. 

Example 1. A metal box without a top has inside dimensions 
6 X 4 X 2 ft. If the metal is 0.1 ft. thick, find the actual volume of 
the metal used and compare it with the approximate volume found 
by using the differential. 

The actual volume is AF, where 

A7 = 6.2 X 4.2 X 2.1 - 6 X 4 X 2 = 54.684 - 48 = 6.684 cu. ft. 
Since V = xyz, where x = 6, y = 4, 2 = 2, a. 

dV = yzdx + xz dy + xy dz \ 

= 8(0.2) + 12(0.2) + 24(0.1) = 6.4 cu. ft. ^ \ 

/ot»y 

Example 2. Two sides of a triangular piece of ^ ^ 
land (Fig. 29) are measured as 100 ft. and 125 ft., Fia. 29. 
and the included angle is measured as 60°. If the 
possible errors are 0.2 ft. in measuring the sides and 1° in measuring 
the angle, what is the approximate error in the area? 

Since A = ^ixy sin a, 

dA = K (2/ sin a dx + x sin a dy + xy cos a da), 
and the approximate error is therefore 

dA = I [l25 (^) (0.2) + 100 (^) (0.2) 

+ 100(125) (I) jIq] = 74.0 sq. ft. 

PROBLEMS 

1. A closed cylindrical tank is 4 ft. high and 2 ft. in diameter (inside 
dimensions). What is the approximate amount of metal in the wall 
and the ends of the tank if they are 0.2 in. thick? 

2. The angle of elevation of the top of a tower is found to be 30°, with 
a possible error of 0.5°. The distance to the base of the tower is found 
to be 1000 ft., with a possible error of 0.1 ft. What is the possible error 
in the height of the tower as computed from these measurements? 
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3. What is the possible error in the length of the hypotenuse of a 
right triangle if the legs are found to be 11.5 ft. and 7.8 ft., with a 
possible error of 0.1 ft. in each measurement? 

4. The constant C in Boyle’s law pv = C is calculated from the 
measurements of p and v. If p is found to be 5000 lb. per square foot 
with a possible error of 1 per cent and v is found to be 15 cu. ft. with a 
possible error of 2 per cent, find the approximate possible error in C 
computed from these measurements. 

6. The volume y, pressure p, and absolute temperature T of a perfect 
gas are connected by the formula pv = RTy where is a constant. If 
T = 500®, p == 4000 lb. per square foot, and v = 15.2 cu. ft., find the 
approximate change in p when T changes to 503® and v to 15.25 cu. ft. 

6. In estimating the cost of a pile of bricks measured as 6 X 50 X 4 
ft., the tape is stretched 1 per cent beyond the estimated length. If 
the count is 12 bricks to 1 cu. ft. and bricks cost $8 per thousand, find 
the error in cost. 

7. In determining specific gravity by the formula s = ^ 

where A is the weight in air and W is the weight in water, A can be read 
within 0.01 lb. and W within 0.02 lb. Find approximately the maxi- 
mum error in s if the readings are A = 1.1 lb. and W = 0.6 lb. Find 
the maximum relative error As/s. 

8. The equation of a perfect gas is pv — RT, At a certain instant 
a given amount of gas has a volume of 16 cu. ft. and is under a pressure 
of 36 lb. per square inch. Assuming R == 10.71, find the temperature 
T. If the volume is increasing at the rate of H cu. ft. per second and 
the pressure is decreasing at the rate H lb. per square inch per second, 
find the rate at which the temperature is changing. 

9. The period of a simple pendulum with small oscillations is 

If T is computed using Z — 8 ft. and g — 32 ft. per second per second, 
find the approximate error in T if the true values are Z = 8.05 ft. and 
g « 32.01 ft. per second per second. Find also the percentage error. 

10. The diameter and altitude of a can in the shape of a right circular 
cylinder are measured as 4 in. and 6 in., respectively. The possible 
error in each measurement is 0.1 in. Find approximately the maxi- 
mum possible error in the values computed for the volume and the 
lateral surface. 

39. Total Derivatives. Thus far, it has been assumed that x 
and y were independent variables. It may be that x and y 
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are both functions of one independent variable t, so that z 
becomes a function of this single independent variable. In 
such a case, z may have a derivative with respect to t 
Let z = /(x, y)y where x = <p{t) and y = ^(0; these functions 
are assumed to be differentiable. If t is given an increment Ai, 
then Xj y, and z will have corresponding increments AXf Ay, 
and Az, which approach zero with A^. As in the case when x and 
y were independent variables, 

, Az = -^Ax + -^ Ay + (1 Ax + et Ay. 

Then, 

^ -1- 4- ^ 4_ ^ 

At dx At dy At A< A< 

and 

m n ^ 4- 

dt dx dt dy di 
Moreover, from (39-1) it appears that 

dz^gdx+^ydy 

gives the expression for the differential in this case as well as 
when X and y are independent variables. 

The general case, in which 

^ = f(^h X2y • • • y Xn) 

with 

Xi = ^l(0> ^2 ~ ^2(t')y * * * ; 

can be treated similarly to show that 

cfe _ ^ ^ 4- ... 4- ^ 

dt dxi dt dX2 dt dxn dt 

and 

+ 5 /"^* + ” ■ + 

In case ^ = x (39-1) becomes 

^ 4- ^ ^ 4- 

dx dx dy dx dx ^ dy dx 

This formula can be used to calculate the derivative of a func- 
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tion of X defined implicitly by the equation fix, y) = 0. Let 
z = fix, y), so that 

dz _ df , df dy 
dx ^ dx^ dy dx 

Since z = fix, y) = 0, it follows that dz/dx = 0 and 

dx ^f/^y 

provided that df/dy 0 . 

As an example, let 

+ y'^x — 1 = 0 

define y as an implicit function of x. Then, 

4- y^ 

dx x^ + 2xy 

for all values of x and y for which the denominator does not 
vanish. 

It was noted that the total differential of a function 
z = fixi, X2, • • • , Xn), where Xi = (piit), 

is given by 

It will be proved next that the same formula can be applied to 
calculate the differential even when the variables Xi are functions 
of several independent variables ti, fe, • • • , <m. Thus, consider 

« = X 2 , • • • , Xn), where Xi = fpiih, fe, • • • , Un). 

In order to find the partial derivative of fixi, X 2 , • • • , Xn) with 
respect to one of the variables, say h, the remaining variables are 
held fixed so that fixi, X2, * * * j Xn) becomes a function of the 
single variable fc; Then, 

df _ ^ ^ 4 . ... 4 . 

Sti dxi dti dxz dti dXn 

^=-^^ 4 .-^^ 4 -... 4 - JiL ^ 

dtz dxi dtz dXz dtz dXn dtz % 


df _ d^ , ^ d^ 4 . ... 4 . ££ 2 . 

dtfn dXi dtfn dXz dtfn dXn dtm 
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If the first equation of (39-2) is multiplied by dti, the second 
by dtt, etc., and the resulting equations are added, there results 




+ 

+ 

+ 


or 

(39-3) 


df 1 

'dxi 

dxi^ 

\dti' 

df i 

(dXt 

dX2 ^ 

[dh 


( dXn 

dXn 

\dk 




^ dh 

+S'^+ 


dt^ " 1 “ 




-i- Ai 


+ 


dx, 

dL 

df 


) 
) 


OXn 


This establishes the validity of the formula (38-1) in all cases 
where the first partial derivatives are continuous functions, 
irrespective of whether the independent variables are X\, X 2 , 
* * • , Xn or t\y ^2> * * * > tm* 

An important special case of the formula (39-3) arises in 
aerodynamics and other branches of applied mathematics. 
Consider a function u = /(x, y, z, t) of four variables x, y, and L 
The total differential of u is 


(39-4) 


du 


du j 4 du j X du j . du 


Let it be supposed that x, y, and z are not independent variables, 
but functions of the variable t. In such a case, tt will depend on 
t explicitly, and also implicitly through x, y, and z. Dividing 
both members of (39-4) by dt gives 


dt ~ dx dt"^ dy dt dzdt'^ dt 

On the other hand, if the variables x, y, and z are functions of t 
and of some other set of independent variables r, s, • • • , one 
must replace dxfdt, dy/dt, and dzfdt in the right-hand member 
of (39-5) by dx/dt, dy/dt, and dz/dt, respectively, and du/dt in 
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the left-hand member by du/dL The partial derivative with 
respect to t which appears in the left-hand member differs from 
that appearing in the right-hand member, since the latter is com- 
puted from u= f(x, 2/, z, t) by fixing the variables Xy y, and z 
and differentiating the resulting function with respect to t. In 
order to indicate the distinction between the meanings of the 
two partial derivatives with respect to ty one can write 

Du _ dudx . du dy . du dz . du 
dt ^ dx dt dy dt dz dt dt 

The fact that the total differential of a composite function 
has the same form irrespective of whether the variables involved 
are independent or not permits one to use the same formulas for 
calculating differentials as those established for the functions 
of a single variable. Thus, 

d{u + v) ^ du + dVy 

— V du udo, 
etc. 


Example 1. 11 u = xy + yz + zx, and x = t,y = e-*, and z = cos t, 

du , , .dx , . , .dy , , , .dz 
^ = (y + «)^ + (a: + *)-^ + (x + j/)-^ 

— (6“‘ + COS 0(1) + (< + COS 0("“^~0 + (^ + sin 0 

= -f cos < — cos t — t sin t ~~ e~* sin t: 

This example illustrates the fact that this method of computing du/dt 
is often shorter than the old method in which the values of a:, y, and z 
in terms of t are substituted in the expression for u before the derivative 
is computed. 

Example 2. If /(x, y) — where x = r cos ^ and y — r sin 

then 


df df dx df dy . ■ s , „ , ^ 

= — 2r* cos ^ sin ^ + 2r^ cos ^ sin ^ = 0. 

Also, 


df — 2r dr or df — 2xdx + 2y dy. 
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Example 3. Let z = e*», where x = log (m + v) and y = tan"* (u/v). 
Then, 


dz 

Hence, 

Similarly, 


dz 


dx 


1 ^ dy V 

;r- = aje**', 5 - = — ; — f and 5 - = -t— ; — ;• 

ay ’ du u + V du 


dz dz dx dz dy ye^v 


du dxdu dydu u + v 


dz 


ye^v 


xe^yu 


dv u + V v'^ + 

The same results can be obtained by noting that 
dz = ye*y dx + xe*y dy. 


But 

and 

Hence, 

But 


dx du + ^ dv du H do 

du dv u + V u + V 


, dy 

dy = dM + ^ dt. = 


du — rrn — ; dv. 


+ u^ 


+ u^ 


. du + dv , V du udv 

dz = ye*y ; h xe^y — 5-^ — 5— 

( ye*v , xe^yv \ , / ye*y xe^yu \ . 

u + + uy ^ + t; + uy 


, dz , .dz 
dz = ^dM + ^dr; 


and since du and dv are independent differentials, equating the coeffi- 
cients of du and dv in the two expressions for dz gives 


du " u + V 


yexy 


+ 


xe^yv 
v^ + u^ 


and 


dz ye*y xe^yu 

dv u + + u^' 


PROBLEMS 

1. If = xyz and a; = a cos d, y = a sin d, 2 = kB, find du/dS, 

2. If u = X* — y* and y = r sin d and x = r cos 6, find du/dr and 
du/dB, 

3. If u — xy — yz and x = r4-s, y — r s, z — t, find du/dr 
du/dSf and du/dt. 
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4. If z ** X ^logy/u"^ + v^y and y = tan“^ find dz/du and 
dz/dv, 

b. li z f{x + Uyy + v)y show that dz/dx *= dz/du and dz/dy 
« dz/dv. 

6 . If w » a?®2/ + yH + zHy verify that 


du . du : du , , , 

ai + + “ (» + i/ + 8). 


7. (a) Find du/diy if = e* sin yz and x - t^y y - t — ly z - \/t, 

(6) Find du/dr and du/dBy \iu — — \y^ and x — r sec By y = r tan 0. 

8. (a) Find du/dx and du/dxy if u - x^ + y^ and y == tan x. 
(6) Given V = /(x, t/, «), where x = r cos y = r sin 2 = L Com- 
pute dV/dfy dV/dBy dV/dt in terms of dV/dXy dV/dy, and dV /dz, 

9. If / is a function of u and v, where u = \/x^ + and v = tan“' 

find df/dXy df/dyy and \/ (df/dxy -f (df/dyy. 


40. Euler’s Formula. A function /(xi, X 2 , * • • , Xn) of w 
variables Xi, X 2 , • • * , Xn is said to be homogeneous of degree m 
if the function is multiplied by when the arguments xi, 
« 2 i * • * , Xn are replaced by Xxi, Xx 2 , * * * , Xxn, respectively. 
For example, /(x, y) = x V -y/x^ + is homogeneous of degree 
1, because the substitution of Xx for x and \y for y yields 

XxVa/^* + Again, /(x, is homo- 

y ^ 

geneous of degree —1, whereas f{x,y,z) == z^/^x^ + y^ is 
homogeneous of degree 

There is an important theorem, due to Euler, concerning 
homogeneous functions. 

Euler’s Theorem. If u — /(xi, X 2 , • * • , Xn) is homogeneous 
of degree m and has continuous first partial derivativeSy then ^ 

+ ■ ■ ■ +®»3^ = m/(*x,x*, • ••• ,x„). 

The proof of the theorem follows at once upon substituting 
*1 = Xasi, x'i = \xt, • • • , *» = Xa!„. 

Then, mnce f(xi, xi, • • • , ®») is homogeneous of degree m, 

O * X’y(a;i, ®j, • • • , «»). 
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Differentiating with respect to X gives 

^ ^ + ■ ■ ■ + “ mX«-y(xi, X*, • • • , x»). 

If X is set equal to 1, then xi = x[, 0:2 = ajJ, * * * , Xn = and 
the theorem follows. 

PROBLEM 


Verify Euler^s theorem for each of the following functions: 

(a) y, «) = + 2 x^ 2 ; 

(b) fix, y) = Vy* ~ sin“i 

(e) /(®, a) = 

(/) fix, V) = 

-y/x + y 

y ’ 


(y) /(x, y) = 


(X) /(x, y) — J.2 _ yl' 


41. Differentiation of Implicit Functions. It was noted in 
Sec. 39 that the derivative of a function of x which is defined 
implicitly by the equation /(x, y) = 0 could be calculated by 
appl 3 dng the expression for the total derivative. This section 
contains a more detailed discussion of this method. 

The equation /(x, y) = 0 may define either x or y as an implicit 
function of the other. If the equation can be solved for y to 
give y = <p{x), then the substitution of y = <p{x) in /(x, y) = 0 
gives an identity. Hence, fix, y) = 0 may be regarded as 
a composite function of x, where x enters implicitly in y. If 
u = fix, y), then 

d« = |^dx+|dy = 0, 

so that 


(41-1) 


dx ~ df/dy 





It will be observed that this discussion tacitly assumes that 
fix, y) = 0 has a real solution for y for every value of x. If 
(41-1) is applied formally to x* + y* = 0, it is readily checked 
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that ^ — 5. This result is absurd for real values of x and y, 

dx y 

inasmuch as the only real values of x and y that satisfy a;* + 2 /^ ~ 0 
are x = 0 and y = 0. 

Example 1. Find dy/dx^ if ^x^y^ + x cos y = 0. Here, 

M /if 

^ = 9x^y^ + cos y, ^ = Qx^y — x sin y, 

so that 

% _ 4- cos y 

dx ~ Qx^y — X sin y 

The relation /(a;, y^z) =0 may define any one of the variables as 
an implicit function of the other two. Let x and y be independent 
variables. Then f{x, y, z) — 0 defines z as an implicit function 
of x and 2 /, and 

But 

Therefore, by substitution, 

I + I + g (I <ir + I <!») . 0. 

This can be written as 


\dx ^ dzdx) . dz dy) ^ 

Since dx and dy are independent differentials and the above rela- 
tion holds for all values of dx and dy, it follows that 

dx ^ dzdx 


and 


If df/dz ^ 
(41-2) 
Example 2. 


dy ^ dz dy 

0, these equations can be solved to give 

dz __ dffdx dz __ Ofidy 

^ df/dz ^ ““ df/dz' 

If a?? + 2y* — 3xz « 0, then, by (41-2), 
dz 2x — Zz dz 4y 
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Frequently, it is necessary to calculate the derivatives of a 
function that is defined implicitly by a pair of simultaneous 
equations 


(41-3) 


iSix, y, z) = 0, 
I vix, y, z) = 0. 


If each of these equations is solved for one of the variables, say 
Zj to yield 

z = F{Xj y) and z = $(a:, y), 

then one is led to consider the equation resulting from the 
elimination of z, namely, 

Fix, y) - 4>(a;, y) = 0. 

This equation may be thought to define y as an implicit function 
of X, and one can apply the method discussed earlier in this 
section to calculate dyjdx. 

However, the elimination of one of the variables from the 
simultaneous equations (41-3) may prove to be difficult, and it 
iS simpler to use the following procedure: The differentiation 
of (41-3) gives 

dj ^ ^ dx dy dz 0 
^ dx dy ^ dz 

and 

These equations can be solved for the ratios to give 

y V V y V V 

dx-dvdz ^ ay, 

ax.ay.az 

By Bz Bz Bx Bx By 

from which the derivatives can be written down at once. 

Example 3. Let 

/(a?, I/, 2) ^ X* + 2/2 4- z’* - a® = 0 
and 

Vf z) = — 2/* — 2^2 — 6* = 0. 


2y 

2*1.1 

2z 

2x 

2x 

2y 

-2y 

-44-1 

— 4« 

2x ' 

‘2x • 

~2y 


= — Ayz: 12xz : — Sxy. 


Then, 
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Hence, i 

dy I2xz dz —Sxy . 


Another important case arises from a consideration of a pair 
of simultaneous equations 


(41-4) 


f(x, y, u, v) = 0, 
<P(X, Vy Uy v) = 0, 


which may be thought to define u and v as implicit functions 
of the variables x and y. 

Differentiating (41-4) gives 


(41-5) 


= ||da: + = 0. 


But, since u and v are regarded as functions of x and y, 


j du j , dw j 

du^^dx + ^^dy 

and 

, dv j , dv j 

dt; = — ax 4- -r- dv, 
dx ^ dy ^ 

Substituting for du and dv in (41-5) gives 


\dx ^ dudx^ dv dx) ^ \dy 
,dx dudx^ dv dx) \dy 


, . 


+ 


d<p du , 
du dy 




Since the variables x and y are independent, the coefficients 
of dx and dy must vanish, and this leads to a set of four equations 
for the determination of du/dx, dujdy, dv/dx, and dv/dy. Thus, 
one obtains 



df 



dx 

dv 


difi 

dip 

du _ 

dx 

dv 

dx ~ 

£ 

df 


du 

dv 


dip 

dip 


du 

dv 



§41 


PARTIAL DIFFERENTIATION 


141 


and similar expressions for du/dy^ dv/dx, and dv/dy. It is 
assumed in the foregoing discussion that all the derivatives 
involved are continuous and that 


Example 4. 
then 


If 


du dv 
d<p d<p 
du dv 


9^0. 


X -|- 2^3 ^ ^3 ^ j;3 

+ y -- u* + V* = 0, 


1 

3t;2 

3a;2 

4t;® 

3^2 

3t;2 

--4t6® 

4t;® 


I2{u'h)^ + vV) 


PROBLEMS 

1. Obtain dv/dx^ du/dy^ and dv/dy in Example 4, Sec. 41. 

2. Compute dy/dXj if x^ + y^ — Zxy = 1. 

3. Find dy/dx if 

^2 _ y2 ^ 2-2 = 0, 

x2y — + xz^ — a® = 0. 

4. Obtain du/dx and dv/dy y if 

— xy + V = 0 , 

v&» — + w = 0. 

6. If a; = /(w, v) and y = g{Uy v), then differentiation with respect 
to X gives 

dudx^ dv dx 

n - ^ ~ -L 

du dx dv dx* 


from which du/dx and dv/dx can be computed. Consider the pair of 
equations 

X = U^ -- V^y 
y = UOy 

and obtain du/dXy du/dyy dv/dx, and dv/dy. 

6. Apply the method outlined in Prob. 5 to find du/dXy dv/dx, 
du/dy, and dv/dy, if 


(n\ <»“«+», 

W iff=-Su + 2v: 

) 2a5 *= v* — 
y -= uv. 
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7. U X — r cos d and y — r ain d, find dr/dx and d$/dx, 

8 . li w = uv and 



+ V + X 

— u — y 


0 , 

0 , 


one can obtain dw/dx as follows: Differentiation of w with respect to 

X gives ^ ^ values of du/dx and dv/dx can be 

calculated from (a) by the method of Prob. 5. Find the expressions 
for dw/dx and dw/dy, 

9. If 2 = uv and 

— a; — 2/ = 0, 

^ + Sx + y = Of 

find dz/dx, 

10. If 2 = and 

x = — v^f 

y = uVf 

find dz/dx. 

11. li z - and 

w = r cos 6, 

V — r sin 6^ 

find dz/dr and dz/dd, 

12. If r = (a?* + y^)^^ and d = tan“^ find dr/dx and dB/dx, 

13. (o) Find dy/dx^ if a; sec y + x^y’^ = 0. 

(6) Find dz/dx and dz/dy^ if x^y — sin 2 + 2 ® = 0. 

14. Let u^x + y + z = 0 and v = x^ + y^ + z^ — = 0. Find 

dx:dy:dz. 

16. Find du/dx j dv/dx y du/dy^ and dv/dy, if 


u^ + v^ + 2/* — 2x = 0, 

u^ + v^ — x^ + 3y 0. 

16. Find dw/dx and dw/dyy if w? = u/v and 


X = U + Vy 

2/ = 3t4 + 2t;. 

17. Show that = 1 and ^ = -1, « /(*,!/,*) = 0. 

Note that, in general, dz/dx and dx/dz are not reciprocals. 

18. Find du/dXy if 

w* — — a;* + - 0, 

w + t; — 2a; = 0. 


19. Prove that 


dudy , dv dy 
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' if F{x, y, «, ») = 0 and G(x, y, u, v) = 0. 


20. Show that (g)’ + 1/r* (g)' = (|)* + if * = r cos 0 

and y = r sin 6. 


42. Directional Derivatives. The relation expressed in (39-1) 
has an important special case when x and y are functions of the 
distance s along some curve C,“ which goes through the point 
{Xj y). The curve C may be thought to be represented by a pair 
of parametric equations 

X = x{s), 

y = vis), 


where x and y are assumed to possess continuous derivatives with 


respect to the arc parameter s. 

Let P (Fig. 30) be any point of 
the curve C at which /(x, y) is 
defined and has partial deriva- 
tives df/dx and dfjdy. Let 

Q{x -h Aa;, y + Ay) 

be a point close to P on this 
curve. If As is the length of the 
arc PQ and A/ is the change in , 
Ay, then 




lim 

A«— >0 


As 


gives the rkte of change of / along C at the point {x, y). But 


and 

dx 

ds 

Therefore, 

(42-1) 


ds dx ds dy ds ^ 


lim — = cos Of, 
As 


ds 


lim ^ = sin a. 
As->o As 


^ ^ cos a + ^ sin a, 

ds dx dy 


and it is evident that df/ds depends on the direction of the curve. 
For this reason, df/ds is called the directional derivative. It 
represents the rate of change of / in the direction of the tangent to 
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the particular curve chosen for the point (x, y). 


da dx 


If a = 0, 


which is the rate of change of / in the direction of the x-axis. 
a = 7r/2, 


ds dy 


If 


which is the rate of change of / in the direction of the 2 /-axis. 

Let z = /(x, y), which can be inter- 
preted as the equation of a surface, be 
represented by drawing the contour* 
lines on the x^z-plane for various values 
of z. Let C (Fig. 31) be the curve in 
the xy-plane corresponding to the value 
z — and let C + AC be the neighbor- 
ing contour line for » = 7 + Ay. 
Then, Af/As = Ay /As is the average 
rate of change of / with respect to the 
distance As between C and C + AC. 







Fig. 31. 


Apart from infinitesimals of higher order, 


An 

As 


= cos 




where An denotes the distance from C to C + AC along the nor- 
mal to C at (x, y), and ^ is the angle between An and As; hence, 
dn/ds = cos Therefore, 


(42-2) 


da dn ds dn 


This relation shows that the derivative of / in any direction may 
be found by multiplying the derivative along the normal by the 
cosine of the angle ^ between the particular direction and the 
normal. This derivative in the direction of the normal is called 
the normal derivative of /. Its numerical value obviously is the 
maximum value that df/ds can take for any direction. In applied 
mathematics the vector in the direction of the normal, of magni- 
tude d//dn, is called the gradient. 

Example. Using (42-1), find the value of a that makes df/ds a 
maximum, considering x and y to be fixed. Find the expression for 
this maximum value of 
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Since df/da = /« cos a + /» sin a, 

£{t) ° -/*«««+/» cos «. 

The condition for a maximum requires that 


tan ai = 7 * 


or ai = tan 


-1 


Using this value of a\, 


dn 


VU+f,^ 
= V7JT77. 


+ fv 




The relation (42-2) can be derived 
directly by use of this expression for 
df/dn. If a (Fig. 32) gives any di- 
rection different from the direction 
given by ai, then 

df M , p • 

^ cos a +fy sin a. 



C+dC 


Fiq. 32. 


But a = ai — so that 




= /x(cos ai COS ^ 4“ sin ai sin -l-/y(sin a\ cos ^ — cos ai sin ^). 


Since 


cos ai = ^ 




and 


sm ai = 


A 




^ — cos ^ + /> — 


Vfi+n 


vn+n 

sin ^ 


+ /i 


fv .f. _ / /» 


cos ^ — /y 


= COS ^ = V5l"+72 cos 4 > 

vn +fi 


sin yf/ 


^ , 


»OBl.BUS 


1 . Find the directional derivative of /(x, y) = + sin xy at 

(l,ir/2), in the direction of the line making an angle of 45® with the 
x-axis. 
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2. Find 




if a; = r cos B, y ^ t sin B, and / is a function of the variables r and B. 

3* Find the directional derivative of fix, y) = x^y -f e*'* in the 
direction of the curve which, at the point (1, 1), makes an angle of 

30° with the a;-axis. 

4. Find the normal derivative of 
f(x, y) = a?2 -f 2/2. 

43. Tangent Plane and Normal 
Line to a Surface. It will be re- 
called that 



Ax -{- By + Cz == D 


Fio. 33. 


is the equation of a plane, where 
the coefficients A, 5, and G are called the direction components 
of the normal to the plane. If a, and y (Fig. 33) are the direc- 
tion angles made by the normal to the plane from the origin, then 


cos a = 


Therefore, 


A ^ B 

Va* + b* + c* Va* + £* + c* 

c 

COS 7 = — =====• 

Va^ + B^ + C^ 


COS a:cos /3:cos y = A :B:C, 


If the plane passes through the point (xo, j/o, Zo), its equation can 
be written as 

A{x - Xo) + B{y - yo) + C{z - Zo) = 0. 

There is also a normal form for the equation of a plane, entirely 
analogous to the normal form for the equation of the straight 
line in the plane. This form is 


or 


X QOS a + y C08 + z cos y = p, 

B 


+ 


Va^ + b^ + c^ Va^ + b^ + c^, 


y + 


Va^ + b^ + c^ 

^ D 

" VA^ +~B^ + 
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in which p = D / + B* + C* is the distance from the origin 
to the plane. 

Consider a surface defined by 2 = f(x, y), in which x and y 
are considered as the independent variables. Then, 


_ 


(43-1) 




3/ 


dy 




If xq and yo are chosen, Zo is determined hy z = y)» Let 

Ax = a; — and Ay = y — yo, and denote dz hy z — Zq, Then 
(43-1) becomes 


(43-2) 


“ dxl 


(xo, 1 / 0 ) 


(x - Xo) + 


(.xo. uo) 


(y - yo), 


which is the equation of a plane. If this plane is cut by the 
plane x = Xo, the equation of the line of intersection is 


2 


Zo 


dy 


k*o» yo) 


(y - yo), 


and this is the tangent line to the curve z = /(xo, y) at the point 
(xo, 2 / 0 , Zq). Similarly, the line of intersection of the plane defined 
by (43-2) and the plane y ^ yo is the tangent line to the curve 
z = /(x, yo) at (xo, yo, Zq), The plane defined by (43-2) is called 
the tangent plane to the surface z = /(x, y) at (xo, yo, Zo). 

The direction cosines of the normal to this plane are propor- 
tional to 


dx 


} 

(ajo, yo) 


dy 


; 

(a-o, yo) 


- 1 . 


The equation of the normal line to the plane (43-2) at (xo, yo, 2o) 
is therefore 


(43-3) 


X — Xq ^ y — yo ^ z — Zo 




df 

-1 

dx 

(« 0 , Vo) 

dy 

k*o, yo) 


This line is defined as the normal to the surface at (xo, yo, zo). 
Figure 34 shows the difference between dz = BP' and Az == BQ. 
P(^Q, yo, Zo) is the point of tangency and /2(xo + Ax, yo + Ay, Zo) 
is in the plane z = Zo- PP' is the tangent plane. 
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In case the equation of the surface is given in the form 
F(x, y, z) = 0, 

the tangent plane and the normal line at (xo, yo, zo) have the 
respective equations 


(43^) f 


and 

143-5) 


“ ^o) + T- 

|(»o, yo, zq) 


|(xo, yo, zo) 
+ 


( 1 / - 2/o) 

(aro, yo, zo) 


(Z — Zo) - 0 


X 

— xo y 

1 

0 

1 

1 

- Zo 

dx 

(sco, yo, zo) ^y 

dF 

(xo, yo, zo) 

|(xo, yo, zo) 


These equations follow directly from (41-2). 



Example 1. At (6, 2, 3) on the surface -h = 49, the 

tangent plane has the equation , 


2x1 (x ~ 6) + 22/1 (2/ - 2) + 2z\ (2 - 3) = 0 

1 ( 6 , 2 . 3 ) 1 ( 6 , 2 , 3 ) 1 ( 6 , 2 , 3 ) 

or 

6x + 22/ + 33 = 49. 

The normal line is 

g - 6 _ y -2 _ g ~ 3 
12 ~ 4 ■“ 6 ’ 


Example 2. For (2, 1, 4) on the surface g = x® + — 1, the tangent 

plane is 

g - 4 « 2x1 (x - 2) + 2y\ (y - 1) 

1(2, 1) 1(2, 1) 

or , ' 

4x + 22/ — g = 6. 
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The normal line is . 

X - 2 _ y - \ _ z - 4: 

4 “ 2 -1 ' 

PROBLEMS 

1. Find the distance from the origin to the plane x + y + z = 1 . 

2. Find the equations of the tangent plane and the normal line to 

(а) 2x^ + 32/2 + 4z2 = 6 at (1, 1, H); 

4/2 2;2 

(б) j + l -|g = 1 at (4,3,8); 

(c) lat (xqj yof 2o)j 

id) x* + 22/2 - 22 = 0 at (1, 2, 3). 

3. Referring to (43-4), show that 

^ dF dF dF 

cosa:cos^:cos7 = 

where cos a, cos cos y are direction cosines of the normal line. 

4. Show that the sum of the intercepts on the coordinate axes of any 
tangent plane to x^^ + y^^^ + z^'^ == is constant. 

44, Space Curves, It will be recalled that a plane curve C 
whose equation is 

(44-1) y = fix) 

can be represented in infinitely many ways by a pair of parametric 
equations 

(44-2) X = x(0, 

y = y(0 

so chosen that when the independent variable t runs continuously 
through some set of values h ^ t ^ h the corresponding values 
of X and y, determined by (44-2), satisfy (44-1). 

For example, the equation of the upper half of a unit circle 
with the center at the origin of the cartesian system, 

y =* VI - icS 

can be represented parametrically as 

x = cos 
y « sin tf 


(0 ^ ^ ^ tt), 
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or 

X t, 

u = 1 ), 

or 

X = 2t, 

y = Vl — (0 ^ 3^). 

Similarly, a space curve C can be represented by means of a 
set of equations 

! x = x{t), 

y = 2/(0, 

2 = z{t) 

so selected that when t runs through some set of values the 
coordinates of the point P(x, z)y defined by (44-3), trace out 
the desired curve C. 



It will be assumed that the functions in (44-2) and (44-3) 
possess continuous derivatives with respect to t, which implies 
that the curve C has a continuously turning tangent as the point 
P moves along the curve. 

Let Pixoj yoy Zo) (Fig. 36) be a point of the curve C defined by 
(44-3) that corresponds to some value to of the parameter ty 
and let Q be the point {xo + Aa;, yo + Ay, zo + As) that cor- 
responfjs to ^ + A^. The direction ratios of the line PQ 
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joining P and Q are 

Ac' Ac* Ac At* At* At 

If At is allowed to approach zero, Ax, Ay, and Az all teiid to zero, 
so that the direction ratios of the tangent line at P(xo, yo, Zo) 
are proportional to {dx/dt)t^to-idy/dt)t^to'»{dz/dt)t^to- Hence, the 
equation of the tangent line to C at P is 

X — xo _ y — yo _ z — Zo 

x'ito) ” V(M “ ^'(^o) ' 

where primes denote derivatives with respect to t. 

Example, The equation of the tangent line to the circular helix 

X = a cos t, 

?/ = a sin t, 
z = at, 

\/3 a va 

X! == ^ XI. 

a \/3 a 

’2 2 ® 

The element of arc ds is given by 

{dsy = {dxr + (dt/)2 + {dzy, 

so that the length of a space curve C can be calculated from 



The length of the part of the helix between the points (a, 0, 0) and 
(0, a, ira/2) is 



46. Directional Derivatives in Space. There is no essential 
difficulty in extending the results of Sec. 42 to any number of 
variables. Thus, if w = /(x, y, z) is a suitably restricted func- 
tion of the independent variables x, y, and z, then the directional 
derivative along a space curve whose tangent line at some point 
P(x, y, z) (Fig. 35) has the direction cosines cos (x, s), cos (y, s), 




162 MATHEMATICS FOR ENGINEERS AND PHYSICISTS §46 
/ 

and cos ( 0 , «) is 

g ^ cos ix, S) + cos (y, s)+^ cos (a, s) 

_ dudx . dudy 

dx ds dy ds Bz ds 

The magnitude of the normal derivative to the surface u = const, 
is given by 

S = >/(!)■ +(|)’+(l)’ 

The vector that is normal to the surface u = const, and whose 
magnitude is dujdn is called the gradient of u, 

PROBLEM 

1. Find the equation of the tangent line to the helix 

X == a Gostf y = asintf 2 = at, 

at the point where t = t/A, Find the length of the helix between the 
points < = 0 and t = t/A. 

2. Find the directional derivative of / = xyz at (1, 2, 3) in the 
direction of the line that makes equal angles with the coordinate axes. 

3 . Find the normal derivative of / = + 2 /* + 2 * at (1, 2, 3). 

4 . Show that the square root of the sum of the squares of the 
directional derivatives in three perpendicular directions is equal to 
the normal derivative. 

6. Express the normal derivative (45-1) in spherical and cylindrical 
coordinates, for which the equations of transformation are 

(а) a? = r sin 6 cos <Pj y = r bid. $ sin <p, z = r cos 6; 

(б) a? = r sin $,.y = r cos 2 ~ 2 . 

6. What is the direction of the curve x = t, y — z — t^ eX the 
point (1, 1, 1)? 

7 . Show that the condition that the surfaces /(a;, 2 ) = 0 and 

g(aj, y, 2 ) « 0 intersect orthogonally is that 

^ 

Bx Bx By By ^ Bz Bz 
B. Show that the surfaces 

I j . y^ z^ . 

xyz ^ I and *2 + ^ “ I “ ^ 

intersect at right angles. 
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9, Find the angle between the normals to the tangent planes to 
the surfaces + 2 /* + 2 ** ~ 6 and 2x^ + Sy^ + ^ 9 at the point 

( 1 , 1 , 2 ). 

10. Show that the direction of the tangent line to the curve of inter- 
section of the surfaces /(«, y^z) — 0 and g(x, y^z) — 9 is given by 


cos a : cos : cos y 


f. f. 


f. 

u 


h 

9v 9^ 





9* 9v 


\ 

Hint: Let (xo, yo, Zo) be a point on the curve of intersection, and find 
the line of intersection of the tangent planes to the surfaces at the 
point (xo, 2 / 0 , Zo). 


46. Higher Partial Derivatives. The partial derivatives fx^, 
‘ j of f(xiy X 2 , • • • , Xn) are functions of Xi, X 2 ^ • • • , 
Xn and may have partial derivatives with respect to some or all 
of these variables. These derivatives are called second partial 
derivatives of /(xi, 0 : 2 , * ‘ , Xn)^ If there are only two inde- 

pendent variables x and y, then f(x, y) may have the second 
partial derivatives 



\ _ a*/ _ 

/ 

dx \dx 

/ ~ dx^~ 

Jxx, 

±(^ 

\ dj 

/ 

dy \dx 

) dy dx 

— Jxui 

± 

dx \dy 

iii 

= fyx, 

dy \dy, 

/ 

fyy 


It should be noticed that fxy means that df/dx is first found and 


then 




\dx) 


is determined, so that the subscripts indicate the 


order in which the derivatives are taken. In 


ay ^ ± 

dy dx dy \dx/ 

the order is in keeping with the meaning of the symbol, so that 
the order appears as the reverse of the order in which the deriva- 
tives are taken. 

It can be proved* that, if fxy and fyx are continuous functions 
of X and y, then fxy == fvxj so that the order of differentiation is 
* See SoKOLKiKOFF, I. S., Advanced Calculus, Sec. 31. 
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immaterial. Similarly, when third partial derivatives are 
found, fsyx ~ ^xxy ~ fyxs and fvyy ~ Sv^ “ fvv^i ^ thesc denva* 
tives are continuous. 


Example. If /(x, y) = e**', then 
/» * fy = xe»y, 


1. Verify that 


ay 


f*x = 

fxy = fyx 

PROBLEMS 

ay 


e*y{xy + 1), fyy = 


for 


' dxdy dy dx 
(a) f = cos xy\ (b) f = sin^ x cos y, (c) / = 

2. Prove that if 

(o) fix, y) = log (x* + 2 /*) + tan-» then |^ + = Oj 
(x« + y* + then II + g + g = 0. 


W /(», y, *) 

3. If It = a;2 + 2 /* and 

i y 

ix — r cos a, 


4. If w = /(», 2/) and 


a: = s + 3^, ^ j a*^^ , a^w 

J d*« , dhi 

r sin 0, ar* d0^' 


6 . Use the results obtained in Prob. 6(6), Sec. 41, in order to show 
ax* ” “ m*)/(m* + vy. Find and 

6. If w = wiXf y)f where x = x{Uf v), y = 2 /(w, v), dx/du « a 2 //ai;, 
and dx/dv — ’-dy/du^ show that 

du^ ^ dv^ ■" W a^v L W/ J’ 

7. Show that the expressions 

= (ai) + (ai?) = ai-* + ^* 

upon change of variable by means of a? « r cos a and 2 / “ r sin a become 




d^z , 1 a*^ . 1 ag 
“ air^ 


|.2 ^^2 + ’^3^’ 


8. If 7 *s f(x + ct) + g(x — ct)f where / and g are any functions 
possessifig continuous second derivatives, show that 
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HT 


9. Show that if a: = cos 0 and y = s' sin 0, then 

dW , dW ^ , /dW , d»7\ 

d** + ® Adr* + ddV‘ 

10. If 7i(ic, 2 ) and F 2 (a;, y, «) satisfy the equation 


show that 


^27 527 527 

V27 = zJl 4. z-JL 4. ZjL — 0 
^ ^ dx’^^ dy^ ^ dz^ 


U s Vi(x, y, z) + (x^ + y^ + z^)V 2 (x, y, z) 
satisfies the equation 

^2 ^2 Q 2 

V^’U-O, where V= . + jp + jj.- 


47. Taylor’s Series for Functions of Two Variables. This sec- 
tion contains a formal development of a function of two variables, 
f{Xy y)y in sl series analogous to the Taylor ^s series development 
of a function of a single variable. It is assumed that the series 
obtained here converges to the value of the function /(a:, y), 
but the analysis of the conditions under which this convergence 
will occur is too involved to be discussed in this book. 

Consider /(x, i/), which is a function of the two variables x and 
y, and let it have continuous partial derivatives of all orders. 
Let 


(47-1) X = a + at and y = b + pt, 

where a, b, a, and are constants and Hs a variable. Then 

(47-2) f(x, y) = /(a -f «<, 5 + /3<) = F{t). 

If F{t) is expanded in Maclaurin’s series, there results 

(47-3) Fit) = F(0) + F'(0)t + t^ + <»-!-•••. 

From (47-2) and (47-1), it follows that 

F'(t) =Mx,y)f^+Mx,y)^^ 

== /*(», tf)« + y)fi- 

Then, 

F''(t) = [fUx, y)a +/,.(a:, y)fi] ^ + UUx, y)<x + y)fi] ^ 
“ fxx{Xf y)o£® -f- 2fxy{Xy y)afi + fyyixy 
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and 

r"(t) = [/.»(*, yW + 2/».(x, y)a^ -^fyysix, y)fi^] 

'+ y)a^ + y)afi + fyyyix, j/)/3»] ^ 

= f^sxix, y)a* + 3/,„(x, y)a^P + Sf^^ix, y)ap^ 

' +/»!«/(». 2 /)/ 3 *- 

Higher order derivatives of F(t) can be obtained by continuing 
this process, but the form is evident from those already obtained. 
Symbolically expressed, 


F"w - (4 + s i)V(., ») -.•!(+ ^ ^ ^ 0. 

f"..(0 - (4 + ^ !;)*/(., rt - g + 3..^ jgp 

+ *“^’al0 + ^'0' 

Then, 

F<-><0 - (4 + fl !;)*/(*, »)-<■• g + CJ.-J 44 

where 


■■ r!(n — r)! 

Since t = 0 gives ® = o and y = b, it follows that 

F(0) = /(o, 6), F'(0) = «/.(a, 6) + ^l/,(a, b), • • • . 

Substituting these expressions in (47-3) gives 
Fit) as fix, y) = /(a, b) + [a/*(o, b) + d/»(o, b)]t 

+ + ^oifffxuiO', b) + P%yia, b)] + • • • . 

Since at == x -- a and fit y — b, the expansion becomes 
(47-4) fix, y) = fia, b) + /*(o, b)(* - o) + /,(o, b)(j/ - b) 

+ ^, {/«(o, b)(a: - o)* + ^fsyia, b)(x - a) iy - b) 


+ /«(o, b)(y - b)*] + • • • . 
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This is Taylor's expansion for a function /(x, y) about the point 
(a, 6). 

Another form that is frequently used is obtained by replacing 
{x — a) by A and {y — b) by fc, so that a? = a + A and y == 6 + A. 
Then, 

(47-5) f{a + hy b + k) = /(a, b) 4- /*(u, 6) A + fyia, b)k 

+ 6)A^ + 2fxu(o>9 b)hk + fyy{a, b)k^] + • • • . 

This formula is frequently written symbolically as 
f(a + h,b + k)= f(a, b)+(h~ + k ^f(a, b) 

V 

Example, Obtain the expansion of tan”^ - about (1, 1) up to the 

y 

third-degree terms. Here, f(Xy y) == tan”^ so that 

X 


/(», y) = tan-i % 

/(l, 1) = taxi- 

Mx,y) - 

ed, 1) = - ; 


> 

11 

/»»(*! y) - 0^2 ^ 

/**(!, 1) = 2> 

/»»(*, y) - (j.2 ^ 

/*v(l> 1) “ 

fvv(Xy y) — ^^2 ^ y2^2^ 

/v»(l> 1) “ 


Then, 

^ » 3 - § (» - 1) + 5 (J/ - 1) 

+ • • • • 
PROBLEMS 

1. Obtain the expansion for xy^ + cos xy about (1, 7r/2) up to the 
third-degree terms. 
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2. Expand /(a;, y) = at (1, 1), obtaining three terms. 

3. Expand e* cos y at (0, 0) up to the fourth-degree tferms. 

4. Show that, for small values of x and y^ 


and 


c* sin y = 2 / + xy (approx.), 

y2 

e* log (1 + 2/) = 2/ + - -2 (approx.). 


6. Expand /(«, y) x^y + x^ + 1 about (0, 1). 

6. Expand (1 — about (0, 0) up to the third-degree 

terms. 

7. Show that the development obtained in Prob. 6 is identical with 
the binomial expansion of [1 — {x^ + y^)Y'^- 

48. Maxima and Minima of Functions of One Variable. A 

function f(x) is said to have a maximum at a: = a, if 

A+ s /(a + h) - f{a) < 0, 

and 

A*“ s /(a — A) — /(a) < 0, 


for all suflBiciently small positive values of A. If A+ and A*" are 
both positive for all small positive values of A, then f(x) is said 
to have a minimum at a; = a. 

It is shown in the elementary calculus that, if the function 
fix) has a derivative at x = a, then the necessary condition for a 

maximum or a minimum is the 
vanishing of /'(x) at the point 
X = a. Of course, the function 
fix) may attain a maximum 
or a minimum at x = a with- 
out having /'(a) = 0, but this 
can occur only it fix) ceases to 
exist at the critical point (see 
Fig. 36). 

Let it be supposed that fix) has a continuous derivative of 
order n in some interval about the point x = a. Then it follows 
from Taylor^s formula that 

A+s/(a-fA) -/(a) 



(n — 1)! n! ' 




where 0 < ^2 < 1. Let it be assumed further that the first 
n — 1 derivatives of f(x) vanish at x = a but that (a) is not 
zero. Then 


A+ = + eih) 

n\ 



Since is assumed to be continuous in some interval 

about the point a: = a, + Sih) and — ^ 2 ^) will 

have the same sign for sufficiently small values of ft. Conse- 
quently, the signs of A‘^ and A“ will be opposite unless n is an 
even number. But if f(x) is to have a maximum or a minimum 
at X = a, then A+ and A~ must be of the same sign. Accordingly, 
the necessary condition for a maximum or a minimum of f{x) 
B>t X = a is that the first non-vanishing derivative of /(x), at 
X = a, be of even order. Moreover, since both A^ and A" 
are negative if /(x) is a maximum, it follows that (a) must be 
negative. A similar argument shows that, if /(x) has a minimum 
at X = a, then the first non-vanishing derivative of /(x) at x = a 
must be of even order and positive. 

If the first non-vanishing derivative of /(x) at x = a is of odd 
order and /" (a) = 0, then the point x = a is called a point of inflection. 

Example. Investigate /(x) = x® — bx* for maxima and minima. 
Now, 

/'(x) = 5x^ — 20x®, 
which is zero when x = 0 and x = 4. Then, 

/"(x) « 20x» - 60x2, = 0 , /"(4) = 320; 

/'"(x) - 60x2 _ i20x, /'"(O) = 0; 

/^(x) * 120x - 120, /*^(0) = -120. 

Since /"(4) > 0, /(4) = —256 is a minimum; and since f^(0) < 0, 
/(O) » 0 is a maximum. 
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PROBLEMS 

1. Examine the following for maxima and minima: 

(а) y » + 1; 

( б ) y = - 6 )®; 

(c) y — X + cos X. 

2. Find the minimum of the function y == where x > 0. 

Hint: Consider the minimum of log y, 

3. Show that a? = .0 gives the minimum value of the function 

4 

2 / = e® + + 2 cos X. 

4. Find maxima, minima, and points of inflection, and sketch the 
curves, for the following: 

(а) 2 / = 3a; + 4 sin a; + sin 2x] 

(б) 2 / = 3a; — 4 sin x..+ sin 2a;; 

(c) y — 6x + S sin x + sin 2a;. 

6. Find the maximum and minimum values of the function 
2 / = a; sin a; + 2 cos a;. 

6. Find maxima, minima, and points of inflection, and sketch the 
curves, for the following: 

(а) y xlogx; 

(б) a;fi - (2/ - x^y = 0. 

49. Maxima and Minima of Functions of Several Variables. 

A function of two variables f{x, y) is said to have a maximum at 
(a, 6), if /(a + A, 6 + fc) — /(«, &) < 0 for sufficiently small 
positive and negative values of h and A;, and a minimum, if 
f{a + h,b + k) ~/(a,6) >0. 

Geometrically, this means that when the point (a, 6, c) on the 
surface z = /(x, y) is higher than all neighboring points, then 
(a, b, c) is a maximum; and when (a, 6, c) is lower than all 
neighboring points, it is a minimum point. At a maximum or a 
minimum point (a, b, c) the curves in which the planes x — a and 
2 / =« b cut the surface have maxima or minima. Therefore, 
fx(a, b) = 0 and /^(a, b) = 0. The conditions /» = 0 and /y = 0 
can be solved simultaneously to give the critical values. 

The testing of the critical values for maxima and minima is 
more difflcult than in the case of functions of one variable. 
However in many applied problems the physical interpretation 
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will determine whether or not the critical values yield maxima 
or minima or neither. An analytical criterion can be established 
for the case of two variables in a manner analogous to the method 
used for one variable. By the use of Taylor’s expansion, it can 
be shown that if /*(a, 6) = 0 and h) = 0, then /(a, 6) is a 
maximum if 

^ ^ b) < 0 

with 

fxx(p>} 5) ^ 0 and b) ^ 0, 

and a minimum if 


D — b) /**(«, b)/yy(a, b) < 0 

with 

fxx(a, b) > 0 and Ay(a, b) > 0. 

In case 

flyicL, b) fxxip^j b) > 0, 


f{a, b) is neither a maximum nor a minimum. If 


fxyi^y b) fxxiP'i b)fyy{ci, b) — 0, 


the test gives no information, just as/"(x) = 0 gives no criterion 
in the case of one variable. 

These considerations can be extended to functions of more than 
two variables. Thus, in the case of a function /(x, y, z) of three 
variables, 





is the necessary condition for a maximum or a minimum. 


Example 1. A long piece of tin 12 in. 
wide is made into a trough by bending up 
the sides to form equal angles with the 
base (Fig. 37). Find the amount to be 
bent up and the angle of inclination of the 
sides that will make the carrying capacity 
a maximum. 

The volume will be a maximum if the area of the trapezoidal cross 
section is a maximum. The area is 


A — 

/2-2x 


Fig. 37 . 


A =* 12x sin 6 — 2x* sin $ + sin 0 cos d\ 
for 12 - 2x is the lower base, 12 - 2x + 2x cos 0 is the upper base, 



162 MATHEMATICS FOR ENGINEERS AND PHYSICISTS §49 
and X sin 6 is the altitude. Then, 


and 


* 12a: cos d — 2a;* cos ^ + a;* cos* ^ — a;* sin* 0 
= a;(12 cos ^ — 2a; cos ^ + a; cos* 6 — x sin* $) 

n A 

* 2 sin ^(6 — 2a; + a; cos 6). 


dA/dx = 0 and dA/dS = 0, if sin d = 0 and a; = 0, which, from physical 
considerations, cannot give a maximum. 

There remain to be satisfied 


and 


6 — 2a; + a; cos ^ = 0 

12 cos 6 — 2x cos 6 + x cos* ^ — a; sin* ^ = 0. 


Solving the first equation for x and substituting in the second yield, 
upon simplification, 

cos ^ or 6 — 60®, and a; = 4. 


Since phjrsical considerations show that a maximum exists, a; == 4 and 
6 - 60® must give the maximum. 

Example 2. Find the maxima and minima of the surface 


?! 

a* 



Now, 

^ — ^ y. 

dx c a** dy c 6*^ 
which vanish when x = y = 0. But 


^ _ JL L 

dx* a^c dy* b^c dx dy 


Hence, D « l/a*6*c* and, consequently, there is no maximum or minimum 
at X « ^ « 0. The surface under consideration is a saddle-shaped 
surface called a hyperbolic paraboloid. The points for which the first 
partial derivatives vanish and D > 0 are called minimax. The reason 
for this odd name appears from a consideration of the shape of the 
hyperbolic paraboloid near the origin of the coordinate system. The 
reader will benefit from sketching it in the vicinity of (0, 0, 0). 


PROBLEMS 

1. Divide a into three parts such that their product is a maximum. 
Test by using the second derivative criterion. 
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2. Find the volume of the largest rectangular parallelepiped that can 
be inscribed in the ellipsoid 


?! 4- _L 5! 

^ c* 


= 1 . 


3. Find the dimensions of the largest rectangular parallelepiped that 
has three faces in the coordinate planes and one vertex in the plane 


a ^ 6 ^ c 


4. A pentagonal frame is composed of a rectangle surmounted by an 
isosceles triangle. What are the dimensions for maximum, area of the 
pentagon if the perimeter is given as P? 

6. A floating anchorage is designed with a body in the form of a right- 
circular cylinder with equal ends that are right-circular cones. If the 
volume is given, find the dimensions giving the minimum surface area. 

6. Given n points Pi whose coordinates are (xi^yi^Zi), (i = 1,2, 
• • • , n). Show that the coordinates of the point P(a:, t/, 2 ), such that 
the sum of the squares of the distances from P to the Pi is a minimum, 
are given by 



60. Constrained Maxima and Minima. In a large number of 
practical and theoretical investigations, it is required that a 
maximum or minimum value of a function be found when the 
variables are connected by some relation. Thus, it may be 
required to find a maximum of u = f{x, y, z), where x, t/, and z 
are connected by the relation y, z) = 0. The resulting 
maximum is called a constrained maximum. 

The method of obtaining maxima and minima described in 
the preceding section can be used to sqlve a problem of con- 
strained maxima and minima, as follows: If the constraining 
relation y,z) =0 can be solved for one of the variables, 
say z, in terms of the remaining two variables, and if the resulting 
expression is substituted for z in u — /(x, y, z), there will be 
obtained a function u = P(x, y). The values of x and y that 
yield maxima and minima of u can be found by the methods of 
Sec. 49. However, the solution of ^(x, y^z) ^ 0 for any one 
of the. variables may be extremely difficult, and it is desirable to 
consider an ingenious device used by Lagrange. 
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To avoid circumlocution the maximum and minimum values 
of a function of any number of variables will be called its extremal 
values. It follows from Sec. 49 that the necessary condition 
for the existence of an extremum of a differentiable function 
/(«!, X 2 , • • • , is the vanishing of the first partial derivatives 
of the function with respect to the independent variables Xi^ 
Xzf • • • f Xn> Inasmuch as the differential of a function is 
defined as 

^ + ■ ■ ■ + 

it is clear that df vanishes for those values of xi, x^, ^ , Xn 

for whiSh the function has extremal values. Conversely, since 
the variables Xi are assumed to be independent, the vanishing 
of the differential is the necessary condition for an extremum. 

It is not difficult to see that, even when some of the variables 
are not independent, the vanishing of the total differential is 
the necessary condition for an extremum. Thus, consider a 
function 

(60-1) M = f(x, y, z), 

where one of the variables, say z, is connected with x and y by 
some constraining relation 

(50-2) <p(x,y,zy=0. 

Regarding x and y as the independent variables, the necessary 
conditions for an extremum give du/dx = 0 and du/dy — 0, or 

dx~ dx'^dz^ ~ * 
du ^ y _ - 

dy dy dzdy 

Then the total differential 




' I '**+1 


and since the expression in the parenthesis is precisely dz, it 
follows that 


(60-3) 
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The total differential of the constraining relation (50-2) is 

(50-4) + + = 

Jjct this equation be multiplied by some undetermined mul- 
tiplier X and then added to (50-3). The result is 

(to + ^ If ) ^ If ) + (i + ^ 

Now, if X is so chosen that 

+ = 0 

dx^^dx 

^-^X^=0 

dy^ dy 

^ -l-X^ = 0 

vix, y, z) = 0 , 

then the necessary condition for an extremum of (50-1) will 
surely be satisfied. 

Thus, in order to determine the extremal values of (50-1), 
all that is necessary is to obtain the solution of the system of 
Eqs. (50-5) for the four unknowns x, y, z, and.X. The multiplier 
X is called a Lagrangian multiplier. 

Example 1. Find the maximum and the minimum distances from 
the origin to the curve 

5x2 Qxy + 5t/2 — 8 = 0. 

The problem here is to determine the extremal values of 

/(x, y) = x2 + 2/2 

subject to the condition 

(p(Xj y) = 5X2 ^ 5y2 _ g - 0, 

Equations (50-5) in this case become 

2x + X(10x + 62/) = 0, 

2y + X(6x + lOy) « 0, 

5x2 Qxy + .52/2 — 8 = 0. 

Multiplying the first of these equations by y and the second by x and 
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then subtracting give 

6X(y* - X*) » 0, 


so that y ^ ±x. Substituting these values of y in the third equation 
gives two equations for the determination of a?, namely, 

2a;* = 1 and a;* = 2. 


The first of these gives / s x* + j/* = 1, and the second gives / = a;* 
+ 1/2 = 4. Obviously, the first value is a minimum, whereas the 
second is a maximum. The curve is an ellipse of semiaxes 2 and 1 
whose major axis makes an angle of 45® with the a;-axis. 

Example 2. Find the dimensions of the rectangular box, without a 
top, of maximum capacity whose surface is 108 sq. in. 

The function to be maximized is 


/(aj, y, = my 

subject to the condition 

(50-6) xy + 2xz + 2yz = 108. 

The first three of Eqs. (50-5) become 

, I 2/2 + Hy + 2z) = 0, 

(50-7) \ xz + \{x + 2z) « 0, 

I xy + X(2a; + 2y) = 0. 


In order to soJLve these equations, multiply the first by a;, the second 
by y, and the last by z, and add. There results 


or 


\{2xy + Axz + Ayz) + = 0, 

\{xy + 2xz + 2yz) + ^xyz = 0. 


Substituting from (50-6) gives 


or 


108X + %xyz - 0, 

. xyz 
““ 72* 


Substituting this value of X in (50-7) and dividing out common factors 
give 

1 - ^ (y + 2«) = 0, 

1 “• ^ (* + 2a) =0, 

1 - ^ (2* + 2y) » p. 
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From the first two of these equations, it is evident that x = y. The 
substitution of a; = 2 / in the third equation gives z = IS/y. Substitut- 
ing for y and z in the first equation yields a; = 6. Thus, a? = 6, = 6, 

and z — S give the desired dimensions. 


PROBLEMS 

1. Work Probs. 1, 2, and 3, Sec. 49, by using Lagrangian multipliers. 

2. Prove that the point of intersection of the medians of a triangle 
possesses the property that the sum of the squares of its distances from 
the vertices is a minimum. 

3. Find the maximum and the minimum of the sum of the angles 
made by a line from the origin with (a) the coordinate axes of a cartesian 
system; (6) the coordinate planes. 

4. Find the maximum distance from the origin to the folium of 
Descartes + y^ — Saxy = 0. 

6. Find the shortest distance from the origin to the plane 

ax + by + cz = d» 

61. Differentiation under the Integral Sign. Integrals whose 
integrands contain a parameter have already occurred in the 
first chapter. Thus, the length of arc of an ellipse is expressi- 
ble as a definite integral containing the eccentricity of the ellipse 
as a parameter.* 

Consider a definite integral 

(51-1) <p(a) = 

mO 

in which the integrand contains a parameter a and where uo and 
ui are constants. As a specific illustration, let 

^(a)=f^.inaxdx. 

In this case the indefinite integral 

F(x, ax dx == — — + C 

is a function of both x and a ; but, upon substitution of the limits, 
there appears a function of a alone, namely 


T 


«>(“) — I sin ax dx = 


T 


cos ax 


2 

0 




* See Sec. 14. 
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Frequently, it becomes necessary to calculate the derivative 
of the function <p{a) when the indefinite integral is complicated 
or even cannot be written down explicitly. Inasmuch as the 
parameter a is independent of x, it appears plausible that in 
some cases it may be permissible to perform the differentiation 
under the integral sign, so that one can use the formula 


• dOL t/uo dot 

This formula turns out to be correct if /(a:, a) and df{Xy a) /da 
are continuous functions in both x and a. Thus, forming the 
difference quotient with the aid of (51-1), 

(51.2) - ipja) ^ a + Aa) - f{x, a) 

Now the limit, as Aa — > 0, of the left-hand member of (51-2) 
is precisely dtp/da^ whereas the limit of the expression under the 
integral sign is df/da. Hence, if it is permissible to interchange 
the order of integration and calculation of the limit, one has 


(51-3) 


da Juo da 


The restrictions imposed on the function /(a;, a) can be shown 
to be sufficient to justify the inversion of the order of these 
operations. 

Suppose next that the limits of integration ui and Uo are func- 
tions of the parameter a, so that 

In this case, one can proceed as follows: Let 
//(-, a)dx = F(x, a) 

BO that 

(51-4) = /(x, a). 

Then, 

( 51 - 5 ) 

=* F(ui, a) — F(wo, a). 
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Assuming the continuity of all the derivatives involved, one 
can write* 


dip _ dF{ui, a) dui ^ dF(uiy a) ^ dF{uoy a) duo _ dFjuo, a) 
da dui da da duo da da ^ 


which, upon making use of (51-4) and (51-5), becomes 


g - /(«., ^ - /(«., 0 )^ + 1 .) - ?(«., »)] 

- /(»., ^ - /(«., «) ^ + £ M «) & 


The partial derivative appearing in this expression means that 
the differentiation is to be performed with respect to a, treating 
Wo and Ui as constants. Hence, making use of (51-3), 


(51-6) 


da 


N dui 


/(mo, a) 


duo 

da 


+ 


Mo(a) 


«) 

da 


dx. 


This formula is known as the formula of Leibnitz, and it 
specializes to (51-3) when ui and uo are independent of a. The 
validity of this formula can be established under somewhat less 
restrictive hypotheses,! but the limitations imposed on the func- 
tion f{Xy a) in the foregoing discussion are usually met in prob- 
lems arising in applied mathematics. 


Example 1. 


Find if <p{a) 


/ 2a 
-a* 


'2a _ 


dx. 


Then 


da J-a* 

/•2a 
"" J-a* 


0^2 

e da; — (—2a) + 6-^2) 

2r‘* — 

<-Ux + 2ae-^^ + 2e-\ 


Example 2. Formula (51-3) is frequently used for evaluating definite 
iiftegrals. Thus, if 

<p(a) = r log (1 + a cos a;) dXy 

* See Sec. 39. 

t See SoKOLNiKOFF, I. S., Advanced Calculus, Sec. 39, p. 121. 
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= X" =■ (l - l+icos J ^ 

1 / . 1 ,, a + COS X \\ir 

« - ( a; H > ) 

a V \/l — 1 + a cos x/\^ 

If 1 /• ^cc — 1 . + 

aL \/l--a*\ 1 — a 1 + a/ J 


Therefore, 


. i / , \ ^ I _ 

a V \/l — av a 0£\/l — 

= IT ^log a + log ^ ^ + ' 


^(a) = IT log (1 + Vl a*) + c. 

But, when a == 0, 

^(0) = log 1 dx = 0. 

Hence, 

0 = TT log 2 + c and c = — tt log 2, 
and 

via) = IT log 

PROBLEMS 

/ ir/2 

Q sin ax dx by using the Leibnitz formula, 
and check your result by direct calculation. 

2. Find dtp/da^ if ^(a) == Jj (1 — a cos x)^ dx. 

J f*a* ^ 

Q tan”*^ ^ dx. 

4. Find dtp/da, if <p{a) = Jq tan (x — a) dx. 

5. Find d^pldx^ if <p{x) = Jq* \/x dx. 


6. Differentiate under the sign and thus evaluate 

by using II ^ > 1* 

^ » JO a — cos X • (a* — 1)^ 


dx 

(a — cos x)® 


if a2 > 1. 


7. Show that 


log (1 — 2a cos X + a*) dx « 0, if a* ^ 1 

« IT log if a* ^ 1. 
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8. Verify that 

1 r* 

y "" Jc Jo /W sin kix - a) da 
is a solution of the differential equation 

g + =/(*). 


where A? is a constant. 
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62. Definition and Evaluation of the Double Integral. The 

double integral is defined and geometrically interpreted in a 
manner entirely analogous to that sketched above for the simple 
integral. Let /(aj, y) be a continuous and single-valued function . 
within a region B (Fig. 38), bounded by a closed curve C, and 
upon the boundary C. Let 
the region R be subdivided ^ 
in any manner into n sub- 
regions ARi, AR 2 , • • • , AJ2„ 
of areas Ail 1 , A A 2 , • • • ,Aiin. 

Let (f», Vi) be any point in 
the subregion ARi, and form 
the sum 

/({<, Vi) AA<. 

The limit of this sum, as n — > 00 and all Ail^ 0, is defined as 
the double integral of f(x, y) over the region B. Thus, 

(52-1) lim V /({,, r,,) Ail, f /(x, y) dA. 

n-* « ^ 

The region B is called the region of integration, corresponding 
to the interval of integration (a, h) in the case of the simple 
integral. The integral (52-1) is sometimes written as 

/ ^ 

In order to evaluate the double integral, it will be simpler to 
consider first the case in which the region R (Fig. 39) is a rec- 
tangle bounded by the lines x = a, x = by y = c, y = d. The 
extension to other types of regions will be indicated later. 
Subdivide R into mn rectangles by drawing the lines x = Xi, 

X == Xz, • • • y X ^ Xn-\y y ^ yiy 2/ = 2/2, * ' * , 2/ = 2/m-l 
Define Aa;, s a;, — Xi^i, where Xq = a and Xn — b, and define 
Ayf s yj — where yo = c and ym = d. Let ABa be the rec- 
tangle bounded by the lines x = XiSiy x — Xiy y — y =*= y^. 
Then, if the area of AR^ is denoted by AAny 

AAii = Axi Ayj. 


i-l 


Fia. 38. 
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Formula (167-3) is known as 
the trapezoidal rule, for it gives 
the value of the sum of the areas 
of the n trapezoids whose bases 
are the ordinates yo, Vu 2 / 2 , • * • , 
JL yn- Figure 157 shows the six 
trapezoids in the case of n = 6. 

If m ?= 2, (167-1) becomes 

£ [yo + XAyo + A*yo] dX 

= 2j/« + 2Ayo + |(|-2^AVo 

= 2yo 4- 2(yi — S/o) + g ( 2/2 - 2yi + yo) 

1 ,4 ,1 

= 3 yo + g + g 3/2, 



or 

(167-4) 


I 3/ d® = 5 (j/o + 43/1 + yt). 

Zq O 


Suppose that there are n + 1 pairs of given values, where n is 
even. If these w + 1 pairs are divided into the groups of three 

pairs with abscissas X 2 if X 2 x+i, X 2 i+ 2 f 

then (167-4) can be applied to each group. Hence, 


f^'ydx = f^'ydx + £ydx + ■ ■ ■ + £>d® 

d d 

= 3 (3/0 + 43/1 + 2/2) + 3 (3/2 + 43/* + yt) 

+ • * • + ^ (Vn-i + 43 /„_i -+• yn) 

= ^ [ 3/0 + 3 /» + Mvi + 3 /* + * ‘ ■ + Vn-l) 

+ 2 (ya + 3/4 -f • • • -I- yn-t)]. 

Formula (167-6) is known as Simpson’s rule with ?» = 2. 
Interpreted geometrically, it gives the value of the sum of the 
areas imder the seeond-degree parabolas that have been passed 
through the points {xk, y 2 <), (® 2 i+i, Vk+i), and (® 2 <+ 2 , 3 ^ 24 + 2 ), [i = 0, 
1, 2, • • • , (n - l)/2]. 



